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PREFACE 


This study is a technical manual describing the operation of 
a quantitative computer model which simulates the Canada Pension Plan. 
The Ontario Government has developed this model, called PENSIM, in 
order to examine in depth the implications of alternative proposals 
for amending the Plan with respect to contributions, benefits and 
fund levels. The model has the capacity to test quickly a wide 
range of alternative amendments to the CPP, including the proposals 
advanced by the federal government in the 1970 white paper, Income 
Security for Canadians. In this way, it can contribute significantly 
to the achievement of a federal-provincial consensus on how the 


Canada Pension Plan should be restructured for the future. 


The development of the PENSIM system and preparation of this 
study was undertaken in the Taxation and Fiscal Policy Branch of the 
Ministry of Treasury, Economics and Intergovernmental Affairs. We 
wish to acknowledge the cooperation of officials in the federal 
Department of Insurance, who made available the actuarial data 


necessary to construct the simulation model. 


H. I. Macdonald, T. M. Russell, 

Deputy Minister, Executive Director, 

Ministry of Treasury, Economics Taxation and Fiscal Policy 
and Intergovernmental Affairs. Division. 


September, 1972. 
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I Introduction 


The 1970 federal white paper, Income Security for Canadians, 


set out proposals for amending the Canada Pension Pia.” These 
proposals involve new contribution levels and revised benefit streams, 
hence would have a major impact both on taxpayers and pension recip- 
ients and on the future size of the CPP fund. It is of little value 
for the provinces to discuss these proposed CPP amendments or alter- 
native amendments with the federal government without having a 
capacity to quantify and assess the implications of changes on the 
future operation of the Plan. Accordingly, the Ontario Government has 
developed a quantitative computer model, called PENSIM, to test the 
implications on contributions, benefits and fund levels of a broad 


range of alternative amendments to the CPP. 


The purpose of this study is to set out in detail the design 
and beenagendal dimensions of this pension simulator system. PENSIM 
uses a methodology which is similar to that employed by the federal 
Department of Tpeusanves but in addition provides great speed and 
flexibility in testing alternatives other than those advanced by 

the federal government. In Ontario's view it is necessary to go 
beyond the scope of testing only the federally proposed amendments. 
Rather, a comprehensive approach should be followed which reviews 
a wide range of alternatives. PENSIM provides a common vehicle which 
can be used by all governments for testing and assessing policy 
alternatives so that the final amendments to the CPP produce the 
best results in terms of benefits, contributions and funding in 


the years ahead. 


aren EEE EE REE ESSERE! 


1. Hon. John Munro, Income Security for Canadians (Ottawa: Department 
of National Health and Welfare, 1970) 


2. Department of Insurance, Canada Pension Plan: Actuarial Report as 
at December 31, 1969. 


7 ae 

=: ve werentnon vis 

Sathe ai7ane ha ae 

mee ft 68: a9 YANe ated 7 iene oaee? nara yet er 
s eA 10g" Hines) © - EDN IN hy Tiel aap sew ed 


? ! 
= a rhnhe ule Vy leeds Sha® 2d’ chpbe! sateen a 
: ‘e ; i ay 7 | Tie : 7 aun! Y, ok ¥ ‘ts 78 voewiige, 
’ ” = 
' Py Tay os 
es o etn.® 0 apeey 
7 7 
——* 1 - 
a § u UV Te. a ow 
yt 
_ 
3 (e [3 See a2 bo 
; ) a7 dl 
a? 1 & 
ee ee an, 
su 62 
: 1 : 
yim oc. 5S lias : 
1 . 
. ; sve rok ly guheqm ask a 
. _ 7 _ 
£ eos? oa 7 
| } “#2 pire . - 
- _ xr 2. 
| . oa 
vy / ‘8 ple é pose oo ly ee ” 
aa 
6 , vt y) f 3 
at 
; ee i tae 
( 
i = — 
x smpatac \ewnry pists 1) S20) sea oe et, a 


* | (ce owas f<s tilaed dace rast Te = 


Sin vey {eed 22.8 Sg taata e tne r ‘idiaicedk>>, sedi 


at ] gli = 


oa J See ee oe 


II METHODOLOGY 


This technical manual serves as a guide for the use of 
PENSIM. It is not intended to be a complete description of the 
methodology of forecasting CPP contribution and benefit levels. 
Rather, its purpose is to provide sufficient information for the 


user to fully utilize the model to test policy variables. 


PENSIM is comprised of ten interactive programs. Each 
program in the system inputs the results of its computations into 
one or more other programs. Due to the nature of the Canada 
Pension Plan, it is more efficient to segment the model into ten 
programs rather than using one large program. This method facil- 


itates a saving of both time and expense. 


As outlined in the PENSIM flow chart,” each program is 
categorized by the level in which it operates. Each level delin- 
eates the amount of interaction with other programs. A Level 1 
program interacts with almost every other program on a higher 
level; a Level 2 program interacts with higher level programs but 


not with lower level programs. 


If options are being tested that involve high level 
programs, it is not necessary to refer to lower level programs in 
the computation. For example, a change in the disability flat 
rate is computed by a Level 4 program. Interaction and consequently, 
cost, are minimized since only two programs must be run to effect 
the change: DISBN.7 and FUND.9. If an option is being tested 
that involves a low level program, a greater number of programs 
must be used. For example, if a new economic assumption is tested, 


then all programs would come into use to test the implications. 


The advantage of this method of construction over a 
single program is that it is not necessary to run all the programs 
in order to test a change in certain variables. The segmented 


programming of PENSIM also gives faster turnaround. The entire PENSIM 


3. See page 5. 
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can be run in less than fifteen minutes and no program uses more than 


24K of octal core. 


The last program in the PENSIM model, CPPBN.10, is ancillary 
to the first nine programs. It does not affect in any way the output 
of the final program FUND.9. It is used to calculate the monthly 
pensions for each type of benefit to correspond with the options being 


tested. 


The methodology incorporated in PENSIM is virtually identi- 
cal to that developed by the Department of Insurance. Where possible, 
advantage was taken of the computer's speed and efficiency to improve 
the model and thus PENSIM works on a more disaggregated framework 
than the Department of Insurance method. Instead of computing esti- 
mates on a five-year basis, they are computed on an annual basis and 
in some cases, age groups are disaggregated. These changes do not 
substantially improve results but they do permit closer inspection of 


intermediate computations. 


One improvement does make some change in the results. When 
calculating modified/modified average earnings, interpolations are 
made from a non-linear equation. PENSIM uses a mathematical tech- 
nique to produce the equation and interpolate the values which results 
in more accurate computations than those employed by the Department 
of Insurance. However, these differences are barely visible in the 


final calculations. 


The PENSIM model is written in standard Fortran IV and it 
is currently running on a PDP-10 computer. The programs are ali tele- 
type interactive and require responses by the user before they will 
operate. If PENSIM is to run on a non-interactive system or a com-~- 
puter other than the PDP-10, a number of changes would need to be made 
for compatibility. For example, the PDP-10 is an octal machine and 
requires no zeroing of arrays before execution. However, the IBM 360 or 
370 is a hexidecimal machine and requires all the null arrays to be 


zeroed before execution. 
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The remainder of this technical manual describes the opera- 
tion of PENSIM. Each program is a self-contained chapter and listings 


of the programs are found in the Appendix. 
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c OF QUTPUT FILES 
POP.1 INDEX. 2 MMAE. 3 RETBN.4 CONT. 5 
PRIM, RSCUR RSCUR RSCUR 
PRIM MICUR MICUR MICUR 
I PRIF RSPRO RSPRO RSPRO 
PRIF MIPRO MIPRO MIPRO 
N TXFPB TXFPB TXFPB 
RS IND RS PNX MMRSC 
P MIIND MIPNX MMMIC 
TFIND TFPNX MMRSP 
U POINT MMRSC MMMIP 
AVE67 MMMIC MMTFP 
Tt ADJFC MMRSP ONE 
MMMIP OPSHN 
MMTFP 
ONE 
SURVM 
ONE TXFPB MMRSC RTRSC CONT 
TFIND MMMLC RTMIC RSCCN 
0 TFPNX MMRSP RTRSP MICCN 
MMMIP RITMIP RSPCN 
U MMTFP RTTFP MIPCN 
DBRSC TFPCN 
T DBMIC 
DBRSP 
Pp DBMIP 
DBTFP 
U DERSC 
DHMIC 
T DHRSP 
DHMIP 
DHTFP 
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DSRSC 
DSMIC 
DSRSP 
DSMIP 
DSTFP 


TFPWF 


RSCCH 
MICCH 
RSPCH 
MIPCH 
TFPCH 


CPPBN. 10 
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IIT POP.1 


A. INTRODUCTION 


Program POP.1 is designed to interpolate population 
Statistics between five-year intervals using the Lagrangean poly- 
nomial technique. The Department of Insurance provided population 
data for the years 1966, 1970, 1975, 1980, ..., 2025. This program 


is specifically designed to interpolate population figures between 


these intervals. 


B. VECTORS AND MATRICES 
W - contains the interpolating weights calculated in the 
program. 


KSUB - is a vector which supplies the appropriate subscript 
for the vector NO. 


FAC - contains the calculated denominator that is used to 
compute the weights. 


XNUM - contains the calculated numerator that is used to 
compute the weights. 


NO - supplies the interpolation reference points for the 
calculation of the values that are stored in FAC and 
XNUM. 


POPM - input vector of male population data. 
POPF - input vector of female population data. 


POP - interpolated population output data. 


C. INPUT FILES AND DATA STATEMENTS 


Input files PRJM1 and PRJF1 contain the population 
statistics for the nineteen age classes for 1966 for males and 


females respectively. 


PRJM and PRJF contain the remaining population statistics 
for each five-year interval from 1970 through to 2025 for each of 


the nineteen age classes for males and females. 
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The nineteen age groups are: 


O- 4 50 - 54 
5- 9 55 - 59 
10 - 14 60 - 64 
15 ~- 19 65 - 69 
20 = 24 70 - 74 
25 = 29 75 - 79 
30 - 34 80 - 84 
35 - 39 85 - 89 
40 - 44 90 + 
45 - 49 


D. METHODOLOGY 


The weights used in calculating the interpolated values 
are derived by the Lagrangean polynomial method. The mi weights 
are computed by 


t, (€-="t,) 
w. (t= : 


—————, i# j; j = 1,2, ..., n=1 
j II (t, ~ t,) 


or in longer form 


eae (t - t,)(t - t4)(t - t,) Oe age & cg. P 
1 (t, - to) (t, - t,)(t, - ty) fin (t, a ae a 
wt (pn (t- t,)(t- t,)(e- t,) ... (e-tiy D 
2 cr 
(t, - t,)(t, ~ t,) (ty > t,) eek (ty ~ ts 4 D 
Woy (e) ‘ (t - t ct - tp) (t - t,) Ree 6 t) 


it. +1 ti) (t, +1 to) ct, +1 3) ii (ty +1 a 


In the program four weights are calculated over a fifteen- 
year interval. There are four of these intervals in the sixty- 
year time horizon of the data. The second, third and fourth 
intervals all have the same weights. The weights in the first 
interval are slightly different due to the fact that this interval 
contains fourteen years rather than fifteen. The program has been 


written to handle this problem. 
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The calculation of the first set of weights is set out 


below. 


= (t - 5)(t = 10) (t - 15) 
i = - 5)(1 - 10)(1 - 15) 


a a A Le: = Oe 5) 
2 agg. - 1)(5 - 10)(5 - 15) 


(t - 1)(t - 5)(t - 15 
3 sok - 1)(10 - 5)(10 - 15) 


t - 1)(t - 5)(t - 10 
4 "ie - 1)(15 - 5)(15 - 10) 


The expression to calculate the population statistic is 


4 
BE ere Wy X, j=1, 2, 3, 4 


4 “/ 30, 35, 40, 45 
45, 50, 55, 60 


where the « refer to 


Xx, 1966 
X. 1970 
Xi kbs be 
Xis 1980 
Xoo 2025 


An interpolated value in year t in the first interval 
would be calculated by 


P. = Wi % = Wo XK, tb We X09 + Wy, X5 


The vector NO(24) contains the reference points for the 
interpolation. Thus, the first six numbers in this vector are iF 
10, 15, 1, 10, 15 which correspond with the values to calculate 


W. and W.. For the first fifteen-year interval there are 12 points, 


i 2 
4 groups of 3 each, because the program uses a four point inter- 
polation and each weight uses three reference years. The second 


Wlock of twelve numbers in NO(24) refers to the second fifteen-year 
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interval. Since the weights are identical for all intervals other 


than the first, it is not necessary to calculate them. 


The vector KSUB(8) contains the appropriate subscripts 


co reference NO(24) in order to make the interpolation work. 


E. OUTPUT FILES 


POP.1 produces one output file named ONE. This file 
contains all the population statistics for Canada less Quebec for 
the years 1966 to 2025. A complete description can be found in 


chapter XIII. 


F. NOTES 


Problems will be encountered in the use of this program 
1f the eens attempting to interpolate population statistics 
different from those that exist in files PRJM1, PRJF1, PRJM, and 
PRJF. Provincial population statistics can be used but care must 


be taken. 


As the program is written, the x8 must correspond tec 
the years as outlined in the methodology. In addition, the pop- 
ulation statistics must correspond to the fifteen-year intervals. 
For example, assume data exist for the years 1966, 1970, 1975, 
1980, 1985, and 1990. Any interpolated values in the first 
interval 1966-1980 will be correct. However, if the user inter- 
polates a value for 1986, the result will be incorrect. The 
problem occurs because there is not enough data to supply ali the 
reference points for the second interval. An interpolated value 
for 1986 not only depends on values for 1980, 1985 and 1990 but 
also 1995. If there is no value for 1995, then results will be 


nonsensical. 


It is possible to interpolate correct values without 


having data to 2025 but it is necessary to have complete intervals. 
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IV INDEX,2 


A. INTRODUCTION 


Program INDEX.2 projects the pension index and YMPE 
levels used in the PENSIM model. The calculation of the pension 
index is based upon section 20 of the Canada Pension Plan and the 


program will generate values up to the year 2025. 


Similarly, the program calculates the earnings index 
for 1975-2025 under the provisions of the Plan and uses this index 
to produce the YMPE's for the time horizon based upon the assump- 


tion being tested. 


The user initializes the program from the terminal and 


inputs the variables he wants to change (see section E). 


The actual data for 1966-1969 is used as the basis 
for future projections and this information is fixed within the 


program. 


Be VECTORS AND MATRICES 


AVINC - this vector lists the annual increase in average 
earnings used in the calculation of the average 
earnings levels each year and thus the earnings 
index for each year 1975-2025. 


AVERN - this vector contains the projected average earnings 
levels for each year 1966-2025 where 1966 = 1.0 
and is escalated by the appropriate AVINC value. 


CPINC ~ this vector lists the annual increase in the 
consumer price index used in the calculation of 
the pension index. 


CPI - this vector contains the projected consumer price 
index for each year 1966-2025 where 1961 = 100.0 
and is subsequently escalated by the appropriate 
CPINC value. 


ERIND - this vector contains the earnings indices for 
each year 1975-2025 based upon the eight-year 
moving average lagged two years as defined in 
section 21 of the Canada Pension Plan. 


ERBAS - this mumber is the sum of the earnings levels for 
1966-1973 used as the base for calculating the 
earnings index. 
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PNIND = this vector lists the pension indices for each 
year 1967-2025 as calculated in the program where 
1967 = 100.0 for PENSIM purposes. The values are 
output as TPPNX.DAT. 

YMPE - this wector stores the YMPE values for each year 
1966-2025 which are input up to 1975. Beyond 
that year they are the 1975 YMPE multiplied by 
the earnings index and rounded down to the nearest 
$100 as required by the Plan. The values are output 
as TXFPB.DAT. 


C. INPUT FILES AND DATA STATEMENTS 


No input files are used but the actual values for 1966- 


1969 of AVINC, CPI, PNIND and YMPE are input in data statements. 


D. METHODOLOGY, 


1. Pension Index 

The Act sets out the method for calculating the pension 
index in section 20. The values for 1967-1969 are input directly 
as data. Briefly, the pension index is the average of the CPI 
levels for the twelve months ending on June 30 of the previous 
year subject to a 2 per cent ceiling on any annual increase. Thus, ) 
it becomes necessary to project CPI levels for each year and find 
the average CPI level centred over January 1 of the previous year. 
This is done by calculating the CPI level for each year based 
upon an annual increase in CPI which the user inputs for each year 
until 1975. The 1975 increase is used for all subsequent years. 
To find the average CPI centred over January 1, the average of 


the CPI levels for the two previous years is calculated. 


Because of the ceiling on the pension index of 2 per 
cent with a floor of 1 per cent, it is necessary to do an annual 
comparison of the previous year's pension index escalated by the 
ceiling and the CPI level fer the 12-month period ending the 
previous June 30. To maintain the PENSIM flexibility, the user 
may specify any ceiling on the pension index e.g. a 3 per cent 


ceiling. 
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For calculation purposes. the pension indices were 


then adjusted so that the 1967 pension index = 100.0. 


2. Earnings Index and YMPE 


The statute explains how the YMPE levels are to be 
calculated in each year. The value in 1966 and 1967 was set at 
$5,000 and escalated until 1975 by the pension index and there~ 
after by the earnings index. The YMPE levels for 1966-1972 are 
known and thus input as data. The user will then specify the 
YMPE levels for 1973-1975 and similarly he will input the annual 
increase in average earnings for 1970-1975 where the 1975 value 


is used in all subsequent years. 


The YMPE calculation for all years beyond 1975 is based 
on the 1975 YMPE level multiplied by the earnings index and rounded 
down to the nearest $100. Thus, the earnings index must be cal- 


culated for each year 1975-2025 where the -1975 value = 1,0. 


The Plan defines the earnings index as the average of 
the annual average earnings levels for the eight years ending two 
years previous ,divided by the average of annual average earnings 
levels for 1966-1973. Thus, it is necessary to calculate the 
level of annual average earnings for each year to 2025 where 


1966 = 1.0. 


The basis for the earnings index is the average of 
these leveis for 1966-1973 and to simplify the calculation, these 
are merely summed and divided into the sum of the earnings levels 
for the eight years ending two years previous as required by law. 


This ratio is the earnings index for gach year 1975-2025. 
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3. Annual Increase in Average Earnings Index 
This index may be used in program MMAE.3 to escalate 


the 1967 level of average earnings taken from the CPP computer 


output of contributor data. 


E. POLICY VARIABLE INPUT 


As described above the user specifies the YMPE levels 
in 1973-75 at the time he is running the program. The program 
asks the question and the user merely types in his response for 


each year. 


Similarly, the annual increase in average earnings and 
in consumer price index is input by the user at the keyboard for 
each year 1970-1975. He may choose to use the economic assumptions 
on page 23 of the Actuarial Report or he may select his own. The 


values input for 1975 are taken for each subsequent year. 


Finally, the user must specify the ceiling he wishes 


to yse on annual increases in the pension index. 


A sample answer to each question is provided at the 


time the question is posed to assist the yser. 


F, TELETYPE QUESTIONS 


The user must respond to the following questions for 


the program to continue execution. 


1) SPECIFY YMPE IN nnnna; 
Where nnnn is.1973, 1974, and 1975. The user 
should respond with the YMPE's that will be used 
in a) the current plan; b) the proposed plan; or 


c) a new plan to be tested for the years indicated. 


2) SPECIFY GROWTH IN AVERAGE EARNINGS IN nnnn; 
Where nunn is 1970, 1971, 1972, 1973, 1974 and 1975. 


Respond with the growth rates outlined on page 23 


of the 1969 Actuarial Report or with new ones. 
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3) SPECIFY GROWTH IN CPI IN nnnn: 
Where nnnn is 1970, 1971, 1972, 1973, 1974, and 
1975. Respond with the consumer price indices 


outlined on page 23 of the 1969 Actuarial Report 


or with new ones. 


4) SPECIFY CELLING ON PENSION INDEX: 
Respond with the desired ceiling for the pension 


index. 


G. OUTPUT FILES 


The pension index projections generated by this program 
are stored in output file TFPNX.DAT. The YMPE projections are 
stored in output file TXFPB.DAT and the index of increase in 
average earnings is stored in output file TFIND.DAT. These files 


may then be used in running the PENSIM model. 


H. NOTES 


Program Index.2 produces a series of index projections 
as specified by the user. However, PENSIM will normally default 
to one of the four schemes tested in the Actuarial Report i.e. 
moderate inflation/reasonable stability and current plan/ proposed 
plan. The user would normally specify one of these plans in which 
case the output files from INDEX.2 will not be used. This program 
was designed, however, for the user who wishes to choose a further 
set of policy alternatives and specifically where the objects of 
examination are the economic assumptions and hence the indices and 


YMPE levels. 
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V MMAE.3 


A. INTRODUCTION 


The program MMAE.3 calculates modified modified average 
earnings for different age groups of CPP contributors through time. 
The results are stored in a data file for use in a number of other 


programs in the PENSIM model. 


The term modified modified average earnings originates 
from the double modification applied to average earnings for the 


Canadian population less Quebec. 


The first modification results from calculating average 
earnings using only those earnings which are equal to or less than 
the year's maximum pensionable earnings (YMPE). Any earnings in 
excess of this amount are ignored. Thus modifed average earnings 
are less than the real average earnings for the population. This 
calculation is one of the more complex in the PENSIM model and is 


discussed in more detail later in this document. 


The second modification takes place by deflating all 
modified average earnings with an index constructed by dividing 
the 1967 base year YMPE by the YMPE in each year. Modified modified 
average earnings then have a 1967 constant dollar value in terms 
of productivity changes and inflation. This allows the programmer 
to test various economic assumptions to blow up modified average 


earnings through time. 


B. VECTORS AND MATRICES 


YMPE - contains the YMPE's of the option being tested for the 
60-year time horizon of the program. 


AVERN - contains calculated average earnings stored in 11 age 
groups, through 60 years for both sexes. 


INDEX - contains #he reasonable stability or moderate -infla- 
tion indices for growth in average wages depending on 
the option being tested. 
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YBE - contains the results of the calculation of the Year's 
Basic Exemption for 60 years, the method depending 
om the assumption being tested. 


VALU - is a utility array which is overwritten several times 
during the course of operation of the program. Its 
final output is the modified modified average earnings 
(MMAE) which is stored in the same manner as the AVERN 
array. 


The eleven age groups are: 


15-19 45-49 
20-24 50-54 
25-29 535-59 
30-34 60-64 
39-39 65-69 
40-44 


These groups are the only ones necessary since con- 
tributions to the CPP cannot begin before age 18 and 
cease after age 69. The sixty years are the years 
from 1966-2025. The order of sexes is males first, 
females second. 


ADJ - is a vector of adjustment factors which is explained 
below. 


C. INPUT FILES AND DATA STATEMENTS 


There are five data files read in by MMAE.3. The first 
reads the annual YMPE's produced by the program INDEX.2. The input 
file read in keys off the word typed in at the beginning of the 
program. A listing of the files can be found in chapter XIII. 

The second file read in, ADJ, contains twenty-two adjustment 
factors that are used to approximate the income distribution of the 


other age groups interpolated from the 40-44 age group. 


The initial contents of the AVERN matrix are read in 
from the data file AVE67. This file contains two groups of eleven 
items (male and female); one item for each of the eleven age 
classifications. The data consists of the 1967 total earnings for 
each age group and sex divided by the total number of contributors 


in each age group and sex. 
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The values read into the vector INDEX are the growth 
indices which are used to escalate average earnings through time. 
They are based on the growth rates found on page 23 of the 1969 
Aetuavial Report. The index which is used is keyed off by the 
teletype response entered by the programmer. These files are 


listed in chapter XIII. 


The final data file is read in during the execution of 
subroutine SPLINE. The file read in is POINT. The data in this 
file is used to set up the interpolating function for the spline 
operation. The first value, J, denotes whether the first (J = 1) 
or second (J = 2) derivatives of the end points are being supplied. 
The two values for D are the first or second derivatives of the 
end points. The next théfty-one values are the X's and the final 
thirty-one are the Y's to give the thirty-one ordered pairs 
which describe the interpolating function. There are four sets 
of points as described above, each of which has two sets of values 


interpolated from it to produce eight sets of interpolated points. 


D. METHODOLOGY 


The basic operation of the program consists first of 
calculating the average earnings matrix. This is a simple operation 
involving nothing more than escalating age-grouped 1967 average 
Canadian (except Quebec) earnings by an appropriate index. which 
is found on page 23 of the 1969 Actuarial Report. If the user 
wishes to use another index it should be entered in the data 


file TFIND, created by program INDEX.2. 


The matrix AVERN is then divided year-wise into the vector 
YMPE for upper-limit values, and then into the vector YBE for lower- 
limit values, by age and sex. The results are stored in the matrix 
VALU. This gives the YMPE and the YBE as a percentage of average 


earnings by age group and sex. 
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Percentage of Contributors’ earnings less 


than the earnings level to total 


contributions . 


- 19 « 


At this point the Department of Insurance derived four 
graphs to calculate modified average earnings. All the graphs 
were computed on the basis of the income distribution of the 40-44 


age group. There are two graphs for males and two for females. 


The first graph in each sex group shows the functional 
relationship between the dollar value of contributors' earnings 
which are less than the YMPE, expressed as a percentage of total 
contributions, and the earnings levels of the contributors' age 


group expressed as a percentage of average earnings. 


The second graph shows the relationship between contrib- 
utors' earnings groups to’ total earnings, and the earnings levels 
of the contributors' age group expressed as a percentage of average 


earnings. The shapes of these graphs are depicted below. 


The percentages which are stored in VALU are then inter- 
polated in turn from these graphs; both the upper and lower-limit 
values are interpolated for each graph which produces eight sets 
of interpolated values. The percentages stored in VALU can be 
considered as X values from which we are interpolating Y values 


on the graphs. 


MALES 40-44 FEMALES 40-44 


Percentage earnings level Percentage earnings level 
to average eornings to average ecrnings 
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earnings. 


earnings to total 


Percentage 
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MALES 40°44 FEMALES 40-44 


to total esrnings. 


earnings 


Percentage 


Percentage earnings level Percentage earnings level 
to averege eorniags. fo averege edriings . 


To transform the interpolated values into income dis- 
tributions for the other age groups a set of adjustment factors was 


calculated by the Department of Insurance. These factors are 


found in the vector ADJ. There are 11 age groups for each sex for a 


total of 22 factors. This method was utilized instead of con- 
structing graphs for all the age groups because of the rather 
gross nature involved in setting up the graphs and using a hand 
interpolation method to draw the functions. It was felt that an 


adjusted interpolation was more accurate than the ocular method. 


The interpolating procedure used in PENSIM is the cubic 
spline method. This calculation is found in subroutine SPLINE of 
the program MMAE.3. The variables used in the explanation of the 


method are the same as those found in the program. 


The method is as follows: 


th 
in the i interval (X, »X, + v define 


= 21 = 


The general cubic then has the form: 


y(x) = Ys -- A, (x - x,) — By (x - x,)° + Cc, (x - x,)° 
a, (a,)2 (a,)3 


" 


' 
Under the conditions that the values of the y, y , y be the same 


t xX 
at X, , , for the intervals (x, x, 4 1) and (X; 44> Xy 4 2) we 
obtain the following equations: 

Replace x with Xe 4d in equation (1), then 


y (x, + D | - A; -- By + Cc, 


or 


Y, f A, + B, _ Cc, 


where the A,s ‘ i 
Differentiating equation (1) we get 
q 
2 
y (x) = Ph RS ~ x,) + 3C, (x - x,) 


——— 


as a, a, a, (a,)2 


Then for the interval (x, + 1°*4 + >) 


q 
9 Aya a eae 
ss Bae aa | 


(1) 


(2) 


B. and C, are the interpolating coefficients. 


(3.1) 


(3.2) 


Equating (3.1) and (3.2) and eliminating A; with equation 


(2) , 


= Y, = B, & 2c, = HA, doer: 


Differentiating equation (1) again gives: 


W 
y (x) = 2B. + 6C , (x - X;) 


a,2 a,2 as 


"t 
Similarly, we evaluate y (x) for the two intervals 


(x, 5X, 4 1) and (x 41? * 4 2)° 
2 


ge Oye Bye a ne ead: 


where i= 1, 2, Pe eet 


If the first derivatives are specified at x, and x» the 


following set of equations is adjoined to equations (2), 


(3) and (4): 


A, = ary 


Rt aoe 


(3) 


Equating the two we get 


(4) 


(5.1) 


(6.1) 
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If the second derivative is specified the following set 
is used instead: 
1 i Ws 


os ak casa 


“4 Sieeiesaas A Biase ice 


To ensure that the point (x, »¥,) lies on the (n - 1) 


of the cubic we use: 


(5.2) 


(6.2) 


(7) 


Finally, the set of equations (2), (3), (4)4 (alg A 6a8) 


or (5.2), (6.2) (depending whether the lst or 2nd deriv- 
atives are specified respectively) and (7) form a tri- 


diagonal system of 3n - 3 equations. 


To solve the tridiagonal system for the 3n - 3 unknown 
coefficients 


A,» B 


4? od 


4? C y ae eS on pk ae 


a step-by-step procedure is used. 


During the operation of the program the user is unaware 
of the interpolation procedure since it is completely automatic. 


It is probably best to regard it as a black box. 


A word of caution is interjected here regarding the use 
of the input file POINT in the SPLINE subroutine. It may be 
necessary to constrain the end points to the derivatives that have 
been specified. This can be accomplished by adding sets of 
ordered pairs to either, or both end points. The fact that the 
spline interpolation method requires the cubic between points 
X and X, agree with the cubic between points X, and X, 7, 


: ee 
at the point xX, in their first and second derivatives makes this 

constraint necessary. One method to ensure correct specification 
is to have the interpolating function plotted on a CALCOMP Incre- 


mental Plotter or some other similar peripheral device to see the 


exact shapes of the function at its end points. 
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Once the interpolation is completed the 5280 values are 
stored in the matrix VALU(11 x 60 x 8). There are four sets of 
interpolated values for males and four for females. The four sets 
are composed of: 

A} proportion of contributors’ earnings less than the YMPE 
B) proportion of contributors’ earnings less than the YBE 
C) proportion of earnings less than the YMPE 


D) proportion of earnings less than the YBE. 


Let the average wage per age group in the year under 
consideration be E and the appropriate year's YMPE be F. The 
final calculation of modified average earnings is given by the 


expression 


E(C - D) + F (1 —- A) 


1-8 


In the program this calculation takes place in one line 
and simultaneously it is adjusted to the proper income group and 
deflated eb 1967 constant dpitars, The deflation is computed by 
dividing the 1967 YMPE ($5,000) by the YMPE of the year in ques- 


tion, 


E. POLICY VARIABLE INPUT 


During the execution of the program the user is asked a 
number of questions. The user must supply an answer before 
execution will continue. The first question relates to the type 
of plan that is to be tested. There are five possibilities each 


represented by a suitable acronym. 


RSCUR - reasonable stability, current plan 
MICUR - moderate inflation, current plan 
RSPRO - reasonable stability, proposed plan 
MIPRO - moderate inflation, proposed plan 


TXFPB - user-specified option. 
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These are the key words which govern the set of YMPE's 
that are read into the program. They are based on the four sets 
of projections as developed by the Department of Insurance, 
Ottawa. The fifth possibility is TXFPB which allows the user to 
specify any set of YMPE's that he wishes. TXFPB is an output 


file produced by program INDEX.2. 


In addition, these codes govern the index that is read 
in to escalate annual average earnings (AVE67). 


RSCUR 


RSPRO - RSIND - reasonable stability index 
MICUR 

MIPRO sd MIIND - moderate inflation index 
TXFPB - TFIND - user-generated index 


If the user has specified the current or proposed plan 
the proper calculation of the YBE is done automatically. If the 
TXFPB option is specified, the user will be asked whether he desires 


the current or proposed plan method of calculating the YBE. 


F. OUTPUT FILES 


MMAE.3 produces one output file of modified modified 
average earnings per run. The output file produced has the 'MM' 
prefix plus the first three letters of the original code word for 
type of plan typed in. For example, a MIPRO run will produce an 
output file named MMMIP. At the end of execution the name of the 
output file generated will appear on the teletype. The output 
file stores the modified modified average earnings matrix by year, 


by age, by sex. (60, 11, 2) 


G. TELETYPE QUESTIONS 


There is one question posed by the teletype that the 
user must respond to. An additional question is asked if the user 


specifies the TXFPB option. 
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1) TYPE IN PLAN TO BE TESTED: 
Respond with one of the five-letter code words outlined 
in section (£). If the user misspelis one of these 


words, the terminal will reply for example,’ RXPRO?, TRY AGAIN. 


2) YBE - CURRENT/PROPOSED? 
Respond with either CURRENT or PROPOSED. This question 
is asked if the user specifies the TXFPB option in 
question 1 above. The computer must know which method of 


calculating the YBE the user desires. 


H. NOTES 


If the user wishes to use a different YBE it must be 


changed in the program. 


To produce an income distribution for a provincial run, 
a new POINT file must be created. The data to create this file 
is available from the CPP Research and Planning Division. Im addition, 
the files AVE67 and ADJFC must be converted to provincial data 


bases. 
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VI RETBN.4 


A. INTRODUCTION 


RETBN.4 is the fourth program in the PENSIM series designed 
to provide benefit costs in the calculation of the CPP fund pro- 
jections. This program specifically calculates the costs of retire- 
ment pensions and the costs of the lump sum death benefit payable to 
the contributor's estate on his death. As well, RETBN.4 generates 
the death benefit factors that are required in WIDOW.6 and DISBN.7 


for the projection of the earnings-related benefit costs. 


The method used in RETBN.4 is adapted from the valuation 
of the federal Department of Insurance for the Canada Pension Plan as 
at December 3lst, 1969. RETBN.4 uses all the assumptions found in 
the Actuarial Report and was designed specifically to test the four 
schemes i.e. moderate inflation, reasonable stability, current and 
proposed plans found in the Report. These in turn are based on the 


federal government's proposals to amend the Canada Pension Plan as 


described in the White Paper on Income Security, 1970. 


In addition to estimating costs based on the four sets of 
assumptions above, RETBN.4 can be used to assess a variety of policy 
alternatives for retirement pensions, and indirectly, all earnings- 
related pensions. In concert with similar changes in the other 
PENSIM programs, fund projections to the year 2025 can be established 


quickly. 


Policy variables handled by RETBN.4 include the retirement 
pension rate, the drop-out rate of lowest annual earnings ratios, 
the retirement age, the earnings test, federal or provincial data, 
the pension index formula, the earnings index formula, and the cost-— 


of-living escalators. 
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VECTORS AND MATRICES 


PART 


PMMAE 


FILE 


SUM 


DROP 


YRSPRT 


EXPOS 


TEM 


DSDRP 


This matrix contains the labour force participation 
rates by sex, year and plan and is input through a 
data statement. 


This matrix contains the modified modified average 
earnings as calculated in program MMAE.3 by age group, 
sex and year. The data is input from the disk and is 
keyed by a switch at the start of execution which asks 
which plan is being tested. 


This matrix is the participating modified modified 
average earnings and is obtained by multiplying each 
element in MMAE by the appropriate participation rate 
in PART. 


This vector contains the input-output file names where 
five options are available for each of YMPE, MMAE, 
death benefit factors, death benefit costs and retire- 
ment pension costs. This option is keyed from the 
switch at the start of execution. 


This matrix is used to sum the PMMAE's for any age group 
over the contribution years of the age group and sex. 


This vector sums the earnings dropped out for males 
because of the drop-out provisions in the Act. The drop- 
out rate is a policy variable in this program and is 
input from the TTY at the start of execution. The drop- 
out calculation occurs in subroutine MIN. 


This vector contains the years of participation for 
each age group depending on the year of death or retire- 
ment (LASTYR) and the first year of participation in 
the CPP (FRSTYR). 


This vector contains the years of participation that 
are dropped out from the calculation of total pension- 
able earnings in subroutine MIN. Each element is 
obtained by multiplying YRSPRT by the drop-out rate 
input from the TTY. 


This counter is used in a do-loop where each iteration 
is death at a different age. 


This matrix is total pensionable earnings by age group 
and sex and is the difference between SUM and DROP. 


This vector is a temporary one in which is stored the 
product of YRSPRT* non-dropped-out percentage *(1/ 
retirement rate) (input from TTY at start of execution). 
Because of the maturing aspects of the Plan, if the 
years of participation are less than 10, there is 

no drop-out provision and TEM1 is merely YRSPRT * 1/ 
retirement rate. 


This matrix is a special case used only for retirement 
at age 65 and is used temporarily to find an average 
TEM. for that iteration only. 


This matrix is used to increase the factor by a disa- 
bility drop-out factor. This data is input in a data 
statement and is used in the program based on the NR 
iteration. 


sew): aoe gi othe : 


- ioe 


eta? is ae py ' 


“= 


pe iaeaees Se 


“aie 


el, 
oe. 


a ce? q°-Seani 28? ew $3 ay 4) abet afdt 
sos has 200us oat) ets Ve, ere SOL30g i Seoos BhY (=5 
at. 253 Jom Vague sh Miner say odd. aeun 262665 wi cF 
su afi at @& 2 a0G Ph reds Sijes MPa IS ig po ay pine 
o:.% pres wide WS Ase toa yall ee wel. iat ee 
ul Lis LS ga? r SUate ons 3a * ant mast 2cias 
Ae cao) GCE USS UA EIS 44s Ss 

; 
; : ee “ay “ati iT 
; G ee pes Lie & y : SOL Se fyosg 25] Foe 
: aaa rid « is. CPRTEGS)| Tad 
ST Sy) Su Ast) 
i : i pRtnInUS : Feo ett h 
7 2s |” : 0°57 a > ae 
; a) is way te es ig rei 

: ao [Fx ec OIF ee 
r Lat ae ; 
7 ph ara iaées’ &. irish) no 44 | 
@& ¢ Sy hes) & see] G3 
a é ieee q Le Ue 7 Ze il 
al 7 ie e$p7acs > = sa: «) ond ' se -3 
oe eet sd (et mess « sy au 1 
» 66 IGS in ae ne Ol RT Ce cep 
Tawwomwe VEC pf. Je 7 ae iis $} lazen Fite Ses 
ss a 227 ea Bhi Ssem Js 3 Gee 
} ea i" iat Sie Me es 2 ae ee 
“i pan” Clave iMG) Sik £e ose Greco? On 
$243 MASs 104 

- « ¢ 
> oh 7 TO Re 2370 wee i267 eae a 4 < 2h) aly: 
e ea 48 Ges: - Leiceena é iS tor a 2 
gui So rGp2as) 12543 Sf 2a 
“00 6 > 92 m=. 6e¢scees a. Ge) «2 CO eaee psec 
¢ ae). FOG: @. £fat = tT 34661 tice mise 
. sti Ged © ous =. Gaig? 62 tee Seas 
= Ser? 

ra) : ; 
j - { 


Lae 


aii Se 


ae ay 


RTFAC - This matrix is for temporary storage of the retirement/ 
death benefit factors and the column of factors is then 
converted to row form in matrix FACTOR. 


FACTOR - This matrix contains the death benefit factors for the 
calculation of death benefit costs. After this calcu- 
lation, the last six rows are overwritten by minor 
changes to compute the retirement benefit costs. 


ESC1 - This vector contains the first escalation factors which 
are used to escalate factors for age groups 65+. The 
calculations of ESCl elements are based on Department 
of Insurance methods and are made in the program. 


ESC2 - This vector contains the second escalation factors which 
are used in the calculation of death benefit costs, 


SURVRT - This matrix contains the survival rates by age group, 
sex and year and is input from disk file SURVM.DAT. 


POP - This matrix contains the population by age group, sex 
and year and is input from disk file ONE.DAT. 


NMDTH - This matrix contains the number of deaths i.e. (1 - 
survival rate)* population* the appropriate death benefit 
factor for the calculation of the death benefit costs. 


DTHCS1 - This matrix is the summation of death benefit costs for 
ail age groups up to age 65 by sex and is then over- 
written by multiplying the first nine age groups by the 
second escalation factor and later overwritten by the 
total death benefit costs. 


DTHCS2 - This matrix is the summation of death benefit costs for 
all age groups from age 65+ by sex and is then over- 
written by multiplying the last six age groups by the 
pension index. 


PNI.ND - This vector contains the pension index for each year and 
is input from disk depending on the plan being tested. 


’WTCST - This matrix contains retirement pension costs by sex 
and year and can be output in table form or in data form 
for use in FUND.9. 

DREDC - This vector contains factors that reduce the death bene- 


fit costs because of the ceiling of six months imputed 
retirement pension or 10 per ceat of the year's YMPE. 


C. INPUT FILES AND DATA STATEMENTS 


The most important source of data required in program RETBN.4 
is the appropriate matrix of modified modified average earnings pro- 


duced by program MMAE.3. 


Because the retirement and death benefit factors are totally 
dependent on earnings levels, this file must be read into matrix MMAE 


and used as the basis of calculation. 
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Similarly, it is necessary to read in the population 
figures as generated by POP.1 in file ONE. This source of data is 


read into vector POP. 


The appropriate levels of YMPE and the pension index are 
read in from the disk and the selection of the proper file is keyed 


from the terminal by the user in his choice of plan being tested. 


For the plan under consideration, the user will select one 


of the code words listed in chapter V. 


Two other files which input from the disk are sources of 
data on survival rates and reduction factors. The former is SURVRT 
used in the calculation of death benefit costs and the latter DREDC 
is used to account for the ceiling on the death benefit of 10 per cent 


of the year's YMPE. 


Also used as data are the labour force participation rates 
of vector PART; disability drop-out rates (DSDRP) and a series of 
numbers representing the ages of contributors in 1966 (FRSTYR) and the 
ages of contributors in the year of their death (LASTYR), necessary 


to calculate the death benefit rate. 


D. METHODOLOGY 


For this description, program RETBN.4 can be split into 
two distinct parts. The first is the calculation of the retirement 
benefit factors and death benefit factors, and the second part is 
the calculation of the retirement benefit costs and death benefit 


costs based on the earlier calculations. 


Because the factors generated in the first section of this 
program are earnings-related, they are also used in WIDOW.6 and 
DISBN.7 to calculate the earnings-related portion of those supple- 


mentary benefits. 
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The program is initialized by the user from the teletype 
when he selects the policy variables he wishes to use. By typing in 
any of the plans mentioned in chapter V, the user invokes a series 
of switches which are used to call in the other data e.g. modified 


modified average earnings, YMPE levels, pension index values, etc. 


The calculation of the benefit factor for any age group 
involves summing the modified modified average earnings (MMAE) amounts 
for each year of participation. Thus for someone who participates 
in the CPP from age 18 to age 65 this sum will be the total of 47 


earnings amounts. 


To ensure that the projections reflect only those for par- 
ticipants in the Canada Pension Plan, the MMAE's are multiplied by the 
appropriate labour force participation rate depending on the year, 
sex, and policy options being assessed. The result of this operation 
produces a matrix of "participating" MMAE's (PMMAE) for each year, 


sex, and age group. 


To facilitate the summation of these PMMAE's for each age 
group, it is first necessary to do a transformation of this matrix. 
In so doing, the method departs slightly from that of the Department 
of Insurance, though the results are the same. The federal method 
involves a diagonal step-wise addition of these PMMAE's for every 
five-year interval and for each five-year age group. However, be- 
cause PENSIM generates annual data, it is easier and more accurate 
to transform the matrix from a diagonal series through time for each 
age group to a matrix of rows, each of which contains up to forty- 


seven PMMAE's for that particular age group. 


To retain harmony with the federal method, RETBN.4 uses 
age groups as they stood in 1965. Because those people who will not 
be born until many years in the future must be considered, the range 
of these age groups is the group aged 65-69 in 1965 back to the group 
aged -46 to -50 in 1965. Those age groups 10-14 in 1965 to -46 to -50 


stored in the non-positive partitions of the matrices so that these 
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age groups would not be involved in any partial calculations. The 
remaining age groups i.e. 15-19 to 65-69 in 1965 must be treated 
differently because they would not participate in the CPP for the 


full 47 years. 


The transformed matrix of PMMAE's is given the name PINDIV 
and a similar transformation is done on matrix MMAE so that it be- 
comes matrix INDIV. This latter operation is necessary for the sub- 


sequent calculation of earnings dropped out. 


To find the sum of pensionable earnings for any age group 
up to the retirement age, all the values on the row are summed up. 
For those age groups that do not have a positive storage location, 
there will be 47 values. The other age groups will have fewer values 
e.g. age group 25-29 in 1965 will have 40 values to sum i.e. 65 minus 
25. The calculation of death benefit factors is similar except that 
these factors must be calculated not only for age of death or retire- 
ment at 65 but also for death at ages 25, 30, 35, 40, 45, 50, 55, 60, 
and 65. The calculation of the factors for the different ages of 
death are done iteratively in a do-loop (NR). Thus when NR is set 
at value 8, this is the calculation of the death benefit factor at 


age 60; when NR = 7, the calculation is for death at age 55, etc. 


The calculation of the factor is done by summing the years 
of participation in the CPP, (SUM), as described above, and then 


adapting this amount by various policy choices. 


The first such variation is the drop-out provision of the 
Canada Pension Plan. Under the legislation, whenever the participa- 
tion exceeds 10 years, the lowest 15 per cent of annual earnings 


levels are dropped from the sum in calculating retirement benefits. 


This calculation takes place in subroutine MIN. Matrix 
INDIV is duplicated in the subroutine (DUP) and an exposure factor 
is determined for each age group given their years of participation 


(YRSPRT). Thus with a drop-out rate of 15 per cent and a maximum 
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participation of 47 years, the exposure factor would be 7.05 for that 
age group. Fewer years of participation yields a smaller exposure 
factor. For the retirement benefit calculation, any participation 


of less than ten years does not invoke the drop-out feature. 


The calculation of earnings dropped out involves finding 
the lowest annual MMAE and summing these for each year up to the ex- 
posure factor. For example, with an exposure factor of 7.05, find the 
seven lowest MMAE's for the age group, sum them and then add .05 


times the eighth lowest MMAE. 


The subroutine does this in a series of do-loops where the 
lowest MMAE for the age group is determined and accumulated. This 
value is then increased by 5,000 so it will not be selected again and 
then the process is repeated to find the next lowest MMAE. This 
method continues until the lowest MMAE's are summed and the decimal 
part of the exposure factor is added. This amount is cumulated in 


matrix DROP. 


It should be pointed out that no earnings are dropped out 
for females unless the drop-out rate exceeds 25 per cent. This is 
necessary because of the pattern of years of nil earnings. Because 
females will have more years of nil earnings than males, the drop-out 
provision is less important. However, should the drop-out rate in- 
crease substantially -- PENSIM uses all drop-out rates exceeding 25 


per cent -- the drop-out factor becomes more significant and thus 


female earnings will be dropped out in the calculation as well. 


Because subroutine MIN has altered the table of MMAE's in 
finding the lowest values, it is necessary to duplicate the INDIV 


matrix each time the subroutine is called. 


After the drop-out values (DROP) are returned to the main 
program, they are multiplied by the lowest participation rate because 
the calculation of earnings dropped out did not use participating 


MMAE's. Similarly, only .5 of this value is used,to account for the 
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years of nil earnings that would exist during the years that the pen- 


sioner could have contributed, but did not. 


The resulting amount of total pensionable earnings is de- 
termined by subtracting the earnings dropped out, where applicable, 


from the total pensionable earnings calculated previously. 


Another factor that affects the retirement pension is the 
retirement rate. In the statutory formula, this rate is 25 per cent 
and this is the figure usually entered by the user when he initializes 
the program. To account for this factor, or use any alternative re- 


tirement rate, a vector calle TEMl is calculated. 


Each element in this vector is determined by multiplying 
the years of participation for each age group by the percentage of 
years not dropped out (the residual of the exposure factors) and then 


multiplying this product by the inverse of the retirement rate. 


The inverse of this rate is used because the next step 
involves dividing the TEM1L value into the total pensionable earnings. 
In effect, total pensionable earnings are multiplied by the retire- 
ment rate and divided by the years of participation that remain after 
allowing for the drop-out factor. The amount so calculated is then 
multiplied by a disability factor (DSDRP) which allows for the years 
when contributions may not have been made because of disability. 

Once this calculation has been made, the result is the death benefit 


factor for that age group for death at the given age. 


A slight variation occurs in the calculation of retirement 
benefit factors in that the factors for age groups 60-64, and 65-69, 
are averaged and then the disability drop-out rate is applied. This 


is done by creating a temporary vector called TEM. 


The factors are then transformed from columns to rows so 
that further calculations can take place. The resulting matrix 
FACTOR contains the retirement/death benefit factor for each age group 


from 25-29 to 65-69 for every five years from 1970 to 2025. 
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In further calculations, the age group is collapsed from a 
five-year span to the median age. Thus age group 25-29 in the future 
will be represented by age 27%, etc. To achieve these values each 
pair of successive death benefit factors is averaged e.g. the factor 
for age 27% is one-half the sum of factors for age group 25-29 and 
age group 30-34. The factor for age 22% is determined by subtracting 
the .factor for age 25 from age 27% and deducting this difference from 
age 25 to get the resulting factor for age 22%. The matrix FACTOR 
now contains the factors for ages 22% to 67% in five-year intervals 


of age and time. 


The method for the short-run calculation of these factors 
differs slightly and thus some changes are made at this point. 
Similarly, a linear interpolation is used to find the factors for 
every year based on the quinquennial value calculated previously. 
Thus the matrix now contains factors for ages 22% to 67% for every 


year from 1970 to 2025. 


To determine the factors beyond age 65, it is necessary to 
use the value for 67% and escalate these by an amount that allows for 
pension index ratios, the retirement rate, and YMPE levels. This is 
done in the vector ESCl which is the first escalation factor. This 
file is calculated by finding the initial retirement pension for each 
year and dividing it by a ratio of pension indices where the pension 


index in 1967 for the purposes of calculation is 100.0. 


Once these factors are calculated for each year, they are 
multiplied by the retirement/death benefit factor for age 67% in that 
year and the escalated factor is then carried forward to older age 
groups by moving it forward in time. Thus the retirement/death 
benefit factor for age 67% in 1980 say, is multiplied by the first 
escalation factor in 1980 and this escalated value is assigned to 
age group 67% in 1980, age group 72% in 1985, age group 77% in 1990, 
etc. This process continues for every year and expands the matrix 


FACTOR out to the age group 90+. 
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The matrix FACTOR is then output onto the disk for use in 


programs WIDOW.6 and DISBN.7. 


The next step in the method is the calculation of the death 
benefit costs. To find these costs, the number of deaths in each age 
group for every year is determined and multiplied by the appropriate 
death benefit factor. This is done in matrix NMDTH where for each 
year the appropriate mortality rate (1 - survival rate) is applied to 
each age group and sex population figure, and this multiplied by the 


proper death benefit factor. 


This amount is summed for all ages up to age 62%, and this 
sum is then multiplied by the second escalation factor (ESC2) which 
reflects the amount of the initial retirement pension for each year, 
but does not include a pension index ratio because the death benefit 


factors up to age 62% were not escalated. 


The cost amounts for ages 67% to 90+ are also summed and 
this time the total is multiplied by the pension index for the year 
because it was divided out in the calculation of the first escalation 


factor. 


The two resulting sums for all age groups are then combined 
and multiplied by .5 to account for years of nil earnings. This 
amount is further reduced by a death reduction factor (DREDC) which 
allows for the ceiling on the death benefit of 10 per cent of the 
YMPE for the year or one-half of the annual imputed retirement pension. 
This amount is then totalled for both sexes and output onto the disk 


for use in program FUND.9 or may be produced in tabular form. 


The calculation of retirement benefit costs requires over- 
writing that part of the FACTOR matrix for ages 65+ because of the 
slight differences in calculation of the retirement benefit factor. 
The method of escalation for ages greater than 65 is the same as for 
the death benefit factors, and again short-run values for these 


factors must be provided by matrix RFAC. 
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To account for the earnings test in age group 65-69, the 
retirement benefit factor for age 67% is reduced by one-half for 
males and one-quarter for females. These are the values used by the 
Department of Insurance in its valuation. If the user wishes to 
omit the earnings test for his calculation, this section will not be 
called. Similarly, if the user wishes to alter the terms of existing 
earnings test, the programming can be changed to accommodate the new 


actuarial assumptions. 


The calculation of the retirement benefit costs is done by 
multiplying the population for each age group, sex, and year by the 
appropriate retirement benefit factor and then multiplying this amount 
by the pension index to replace this value that was divided out in the 
calculation of the first escalation factor. These amounts are then 
summed for all age groups and for both sexes in each year. A short 
term factor is applied to account for the maturing aspects of the 
Canada Pension Plan until 1975 and the resulting cost amounts are out- 
put onto the disk as data for use in program FUND.9 or in tabular 


form where it is disaggregated by year and sex. 


E. POLICY VARIABLE INPUT 


During the execution of the program, the user is asked 
several policy questions. The first relates to the type of plan being 
tested. The response is one of the five mentioned above in chapter V 


i.e. RSCUR, MICUR, RSPRO, MIPRO, or TXFPB. 


As mentioned earlier, this response automatically keys other 
data sources to the proper sets of economic assumptions i.e. reason- 
able stability or moderate inflation. The data required includes 
the pension index values, the values of reduction factors, and the 


YMPE levels. 


Other policy variables that must be input include the re- 
tirement pension rate, the drop-out rate and the choice of whether 
to include the earnings test in the calculation of retirement benefit 


costs. 
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The final option to the user is the form of his output 


e.g. in data form for use in FUND.9 or in table form by year and sex. 


F, OUTPUT FILES 


Each execution of RETBN.4 produces three output files. The 
retirement benefit costs are prefixed "RT"; the death benefit costs 
are prefixed "DB"; and the death benefit factors are prefixed "DH". 
The suffix to each output file name is composed of three letters which 
will depend upon the plan being tested e.g. for moderate inflation 
with the proposed plan, the retirement benefit costs are output in 


file RIMIP. 


Each file contains values for each year from 1970 to 2025 
and represents total figures when in data form, or a breakdown by sex 


when in tabular form. 


G. TELETYPE QUESTIONS 


1) TYPE IN PLAN TO BE TESTED: 
Respond with one of the five-letter code words outlined in 
chapter V. 

2) TYPE IN DROP-OUT RATE (0-1): 
Respond with a decimal value in the proper range. The 
statutory drop-out rate is .15. 

3) TYPE IN RETIREMENT PENSION RATE (0-1): 
Respond with a decimal value in the proper range. The 
statutory retirement pension rate is .25. 

4) DO YOU WANT EARNINGS TEST? (YES/NO): 
Respond with YES or NO. The earnings test is included 
in the Canada Pension Plan but may be omitted as a policy 
variable. 

5) OUTPUT? (TABLE/DATA): 
Respond with TABLE or DATA. DATA is used when a full 


PENSIM run is being done to accommodate program FUND.9. 
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TABLE is used when a policy change is being tested on the 
retirement or death benefit alone. 

6) ECONOMIC ASSUMPTION? (RS/MI): 
When the user specifies a TXFPB plan, he must tell the 
program whether he wants reasonable stability or moderate 


inflation assumptions. 


This program has a built-in clock function that records 
the time used for execution (generally 15 seconds of CPU time) and 
this is output on the terminal when execution is completed. The 
program also types out the names of the disk files that are generated 


by the program. 


Because of the number of vectors and matrices that are used 
in the program, most of which must be retained throughout, RETBN.4 
requires 24K memory capacity when run on a PDP-10 and usually requires 


approximately five minutes of actual clock time to run. 


ay 


=e : wn 
ticbiee Yeas wis shud sua xce-al ee 


ua m3 “pix BL a poss oy ee ae” on 


leans i Eousiad hips lanes a ead a beanee 


ae 7 
\ a Sages A ; 


he >« ora) oul ot sein en Se eee old SH eng: Si agree 5 


ys | | | " \siiagens sauia Cpe 


i] 
& Pi +. i 
H ® j c 7% 3663) ae be f 
oS [eh Ale. a — 7 Sas So 
; 
MS ste 4 (Tene: fy. S57 ae 
eh Ea 
a 
i i 
’ 7! Ti P = a) @ 
La 
F 
j 
ane -_ 
eae . 
oa 0 
‘gee 
1 : a 
7 O¢ 
- o 
- - 
\ ot 
7 a= 7 
} a 
{ 


VII CONT.5 


A. INTRODUCTION 


Program CONT.5 calculates the total contributions paid 


into the Canada Pension Plan from 1966 to 2025. This program can 


also be used to determine provincial contributions to the CPP for 


the same time span. 


The program has the facility to test changes in the 


contribution rate and different methods of calculating the year's 


basic exemption. The output can be specified in a formatted file 


or a data file for use in program FUND.9. 


B. VECTORS AND MATRICES 


POP > 
MMAE - 
ETADJ - 
PART = 
ADJ > 
LEMP. oo 
CONRTE - 


is the population data matrix read in from data 
file ONE. Parts of this array are overwritten 
several times in the program. 


contains the modified modified average earnings 
matrix read in from the appropriate data file that 
matches the option to be tested. 


contains the earnings test adjustment factors. 


stores the labour force participation rates for 

the population. There are two sets of participation 
rates for males; one for the years 1966 to 1985 

and one for 1986 to 2025. There are two groups 

of four sets of participation rates for females; 

the first group for moderate inflation assumptions 

and the second for reasonable stability assumptions. 
The four sets are used for the year groupings 1966- 
1975, 1976-1985, 1986-2005, and 2006-2025 respectively. 


contains adjustment factors to account for the CPP 
overpayment that employers cannot claim. 


is a utility array. 


contains the contribution rate for every year from 
1966-2025. 


C. INPUT FILES AND DATA STATEMENTS 


The participation rates and adjustment factors stored 


in the matrices PART and ADJ were developed by the Department of 


Insurance and are entered in data statements. 
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The first two data files that are read in are keyed 
off the code word that the user types in response to the first 
question. The modified modified average earnings are read in 
from the output file produced by MMAE.3. The appropriate YMPE's 


are read in from the output file produced by INDEX.2. 


The remmining data are read in from files ONE and OPSHN. 
ONE contains all the population statistics produced by POP.1l. 
OPSHN is an option file that enables the user to test various 
methods of calculating the YBE and contribution rates. The first 
three values read in are associated with the variable names X1YBE, 
KX1, and X2YBE. The first entry, X1YBE, is the value of the flat 
rate YBE in dollars. If, as in the proposed plan, a flat rate YBE 
is to be tested, it would be set at 600. The third entry, X2YBE, 
is the value of the ratio used in calculating the YBE as a constant 
proportion of the YMPE. The second entry, KXl, is a switch that 
allows the user to test either method of calculating the YBE with 
any of the plans. The table below outlines the switch values for 


various YBE options. 


Code Proportional Flat Rate 
RSCUR KXl = 1 KXl = 3 
MICUR KXl = 2 KX1l = 4 
RSPRO KXl = 3 KXl = 5 
MIPRO RX1 = 4 KXl = 1 
TXFPB KXl = 5 KXl = 2 


If the user wishes to understand the mechanics of this 
switch he should refer to statement 999 and the vector FILE in 


the program listing. 


If a flat rate option is being tested, the value of 
X2YBE, the proportional rate must be set at .12 and if a propor- 
tional rate is being tested, the value of X1YBE must be set at 600. 
This step is necessary if the user requires detailed contribution 
tables for output because the correct titles are keyed off these 
values. If the output is a data file which is read in for the fund 


calculation, this step can be ignored. 
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The last sixty values (ten per line) stored in OPSHN 
are the contribution rates. They are stored as decimal fractions 
e.g. 036. This method of entering the contribution rate allows 
the user to sgecify any contribution rate that is desired in any 
year. For a MIPRO run say, the contribution rates would all be 
.036. If the user wishes to test the effects of changing the 
contribution rate to .046 in 1986 and beyond, then the twenty-first 


to sixtieth values would be .046. 


D. METHODOLOGY 


Since the volume of contributions paid into the CPP 
are based on a percentage of contributory earnings for the employed 
population, the basic goal of the program is to estimate contrib- 
utory earnings on an annual basis. The first step towards this 
goal is to calculate the earnings that are exempted from contri- 
butions. This is accomplished by computing the YBE for every 


year by the method that has been specified by the user. 


The labour force participation rates are then multiplied 
by the population statistics to determine the employed labour force 
through time. The working population is subsequently multiplied by 
modified modified average earnings to obtain age grouped earnings 
through time. Since the first age group includes individuals aged 
15-19, it must be scaled to exclude those who are 15, 16 and 17. Two 
scale factors - .39 and..4 - are used, the first for years before 


2000, and the second for the remainder. 


At this point, the age grouped earnings are totalled 
by year. These annual total earnings are still in 1966 constant 
dollar terms because the modified modified average earnings were 
deflated in the program MMAE.3. It is here that the economic 
assumptions enter the calculation of total contributions. The 
grouped earnings are multiplied by a growth index derived by 


dividing the YMPE in 1966 (5,000) into the YMPE of the appropriate 
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year. In other words, this growth index is similar to the earnings 
index which is specified by the user in the program INDEX.2. Thus, 
whatever the assumptions made concerning the earnings index, they 
are reflected in the YMPE's and consequently, the earnings totals 


for each year, 


Annual contributory earnings are then calculated by 
multiplying the YBE by the participating population, totalling 
the values for each year and subtracting them from the escalated 
earnings. An adjustment factor is included to account for CPP 


overpayments. 


Finally, contributory earnings are multiplied by the 
contribution rate vector, CONRTE, to calculate annual contributions 
to the CPP. The remainder of the program sets up the format for 


the output files. 


Ee POLICY VARIABLE INPUT 


The YMPE's used in the program are keyed from the code 
word that is entered in response to the first question in the 


program. 


As explained above in section C, the method of calcu- 
lating the YBE and the contribution rates are specified in the 


file OPSHN. 


CONT.5 produces a maximum of two output files. Each 
output file is optional. The program will ask the user if he 
desires table output which contains the annual contribution rate, 
contributions for males and females, total contributions and the 
growth in total contributions. In addition, the user can specify 
data output for the fund program. The table output will be placed 
in a file called CONT. The data output is placed in a file that 


matches the code word entered at the beginning of the program. 
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RSCUR = RSCCN 
MICUR r MICCN 
RS PRO - RSPCN 
MIPRO = MIPCN 
TXFPB : TFPCN 


The output data file stores total contributions by 


year. 


G. TELETYPE QUESTIONS 


The user must respond to all teletype questions for 


execution to continue. 


1) TYPE IN PLAN TO BE TESTED: 
Respond with one of the five-letter code words 
outlined in chapter V. If the user misspells 
one of these words, the terminal will reply, for 


example, RXPRO?, TRY AGAIN. 


2) OUTPUT? (TABLE/DATA): 
If the tabular output is desired, respond with 
TABLE; if a data file is desired, respond with 


DATA; and if both are desired, respond with BOTH. 


H. NOTES 


CONT.5 requires one change if provincial contributions 
are to be forecast. It is necessary to supply a set of labour 
force participation rates that are appropriate for the province 


in which the contributions are being made. 


VIII WIDOW.6 


A. INTRODUCTION 


The purpose of this program is to generate values for 
widows' benefits paid out under the Canada Pension Plan for the 
period 1968-2025. The method used in computing values is identi- 
cal to the method used by the federal Department of Insurance with 
a few exceptions. These exceptions are due mainly to the fact that 
WIDOW.6 performs the calculations in one-year intervals whereas the 


Department of Insurance uses five-year intervals. 


The program was designed for maximum flexibility so that 
the user may test numerous schemes, including the four schemes out- 
lined in the Actuarial Report (1969). When running the program 
from a terminal, the program generates certain questions which the 
user must answer. The answers supplied by the user instruct the 


program as to which scheme is to be tested. 


B. VECTORS AND MATRICES 


AGEDIS - this matrix contains coefficients for the distribution 
of wives by age for each age group of males. 


COEFF - this matrix stores the number of surviving widows for 
each year of widowhood by age group. 


DEATHS - this vector stores the number of male deaths for any 
given year by age group. 


DTHFC - this matrix contains the male death benefit factors 
created in the program RETBN.4 (see Chapter VI). 


ERBO65 - this vector stores the earnings-related benefits for 
widows 65 years of age and over. 


ERBU65 - this vector stores the earnings-related benefits for 
widows under 65 years of age. 


ESC1 - this vector contains the escalation coefficients for 
each year from 1968 to 2025. 


ESCDBF - this vector contains the escalated male death benefit 
factors for a given year for each age group. 


FAC68 - this vector contains the male death benefit factors for 
1968 for the first 10 age groups. The factors 
are not calculated in RETBN. 4 for 1968-69. Hence 
the federal Department of Insurance figures are inserted. 
1969 factors are interpolated from 1968 and 1970 factors. 
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FFAC - this vector contains the female death benefit factors 
for each year from 1968 to 2025. 


FRBU65 - this vector stores the flat rate benefits for widows 
under 65 years of age. 


MAXPEN - this vector stores the maximum pension receivable in 
each year from 1968 to 2025. 


PNIND - this vector contains the pension indices for each year 
from 1967 to 2025 (1967 = 100.0). 


POP - this matrix contains population figures and age dis~ 
tribution -for males for each year 1966-2025. 


PRPMRD - this matrix contains coefficients indicating the pro- 
portion of males that are married by age group 
for different periods in time. 


PWEWB - this matrix contains coefficients indicating the pro- 


portion of widows eligible for Widows benefits by 
age group for each year 1968-2025. 


PWNRP1 - these vectors contain coefficients indicating the pro- 

PWNRP2 portion of widow beneficiaries not receiving retirement 
pensions by age group for years of widowhood 1968-69 
(PWNRP1) and 1970 + (PWNRP2). 


PMARRY - this matrix contains probabilites that a widow will 
remarry by age group and for different periods in time. 


STOTAL - this vector stores the averages by widow's age group, 
of the number of male deaths weighted by the es- 
calated male death benefit factor for each male age group. 


SURVIV - this matrix contains the quinquennial survival proba- 
bilities of widows by age group for different periods in 
time. 


SURVM - this matrix contains the one-year survival probabilities 


for males by age group for different periods in time. 


YMPE - this vector contains the year's maximum pensionabie 
earnings for each year 1966-2025. 


OAS - this is a 3-element vector containing 
1) the annual flat rate payment for the current 
plan 


2) the ennual fiat rate payment for the pro- 
posed plan beginning in 1973 

3) the user-specified proposed flat rate pay- 
ment beginning in the year specified by the 
user. 


ERPCTU - these vectors are each 3-element vectors containing the 
ERPCTO percentage of the benefit that is earnings-related for 
widows under 65 (ERPCTU) and for widows 65 and over. 
The order is as follows: 
1) the percentage under the current plan 
2) the percentage under the proposed plan to 
begin in 1973 
3) the user-specified percentage beginning in 
the year specified by the user. 
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C. INPUT FILES AND DATA STATEMENTS 


This program reads in numerous data files, some of which 


are dependent on the scheme being tested and some of which are not. 


Those input files which are independent of the scheme 


being tested are: 


The names associated with each of the schemes that the 


program can test are: 


1) RSCUR reasonable stability current plan 
2) MICUR moderate inflation current plan 

3) RSPRO reasonable stability proposed plan 
4) MIPRO moderate inflation proposed plan 
5) TXFPB all other plans 


For each scheme, three additional files are read in. 


For RSCUR: RSCUR 
RSPNX 
DHRSC 


For MICUR: MICUR 
MIPNX 
DHMIC 


For RSPRO: RSPRO 
RSPNX 
DHRSP 
For MIPRO: MIPRO 
MIPNX 
DHMIP 
For TXFPB: TXFPB 


TFPNX 
DHTFP 


Those variables which are initialized in data statements 


in the program are: 

ACEDIS 

FAC68 

FFAC 

PRPMRD 

The user is reminded that before he can run the WIDOW.6 pro- 
gram, he must run the programs POP.1, INDEX.2, MMAE.3, and RETBN.4 


(see chapters III - VI for instructions). 
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D. METHODOLOGY 


As mentioned earlier, the method used in this program is 
identical to that used by the federal Department of Insurance with 
a few differences. These differences are outlined in later para- 


graphs of this section. 


The program uses the year of widowhood as the basis for 
all calculations, performing all calculations separately for each 
year of widowhood. The first task of the program is to determine 
the number of women who become widows in any given year. This re- 
quires the calculation of the number of male deaths in that year 
(calculated as the difference between the size of the male popula- 
tion in that year and the number of males who survive until the 
next year). The number of widows is simply the proportion of 
those males who were married at the time of their death. The next 
step is to compute for this year of widowhood the number of widows 
that are still alive in each of the succeeding years until 2025. This 
is accomplished by applying survival probabilities to the number of 


widows for various age groups and for various periods in time. 


To this point in the program, the methods of this program 
and of the Department of Insurance are identical with the exception 
that the program performs calculations in one-year intervals rather 
than five-year intervals. However, since many of the coefficients of 
probability for survival and remarriage etc. used in the program 
are available only in five-year intervals, values for these coefficients 


for the missing years are interpolated linearly within the program. 


Having computed the number of surviving widows, the actual 
benefits can be calculated. It is at this point that the most strik- 
ing difference between the two methods (WIDOW.6 versus the Department 
of Insurance) occurs. This will become apparent in the following para- 


graphs. 
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As mentioned before, the basis for all calculations is 
the year of widowhood. The Department of Insurance performs cal- 
culations in five-year intervals (1970, 1975, 1980, etc.) with one 
exception (1968). The Department of Insurance divides the ages 
of the widows in any year of widowhood into 16 age groups (under 20, 
20-24, 25-59, ...., 80-84, 85-89, 90 and over). With the exception 
of the first and last, these age groups are five-year groups. 
WIDOW.6 can perform calculations for every year of widowhood from 
1968 to 2025 (one-year intervals). However, the five-year age groups 
for widows defined above are retained in WIDOW.6 since the interpo- 
lated population projections used by WIDOW.6 were prepared using 
these five-year age groups. The use of five-year intervals and five- 
year age groups is an obvious computational convenience. The only 
problem case'is for 1968. WIDOW.6, with its one-year intervals, must 


make certain assumptions before proceeding with the calculations. 


As an example, consider the calculation of the flat rate 
portion of the widows’ benefit. This payment is made only to widows 
under the age of 65 (at age 65, this is replaced by the Old Age Se- 
curity). Hence, the flat rate benefit is concerned only with the 
first ten of the sixteen age groups defined above (from under 20 up 
to 60-64). The flat rate benefit for any age group in a given year 
is the product of the number of surviving widows in that age group 
for that year, the percentage of widows in that age group who are 
eligible for widows' benefits, and the user-specified value of the 
flat rate payment. The total flat rate benefit for that year is the 


sum of the flat rate benefits for each age group. 


For the year of widowhood 1970, the Department of Insurance 
computes flat rate benefits for the years 1970, 1975, 1980, etc. (five- 
year intervals). For 1970, the calculation will include the age group 
60-64. But, in 1975, the surviving widows of this age group will be 
five years older (i.e. in the 65-69 age group) and will receive OAS. 


However, the surviving widows of the 55-59 age group in 1970 will be 
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60-64 in 1975. Similarly, the 50-54 age group in 1970 will be 
60-64 in 1980, and so on. Hence, for the year of widowhood 197%, 
the flat rate benefit calculation will include all age groups up to 
and including 60-64 in 1970, 55-59 in 1975, 50-54 in 1980, and so on. 
Similar statements can be made for the years of widowhood after 

1970 as computed by the Department of Insurance. The year of 
widowhood 1968 is the only problem year. The Department of 
Insurance calculates benefits for 1968, 1970, 1975, 1980, etc. 
Consider those widows in the 60-64 age group in 1968. In 1970, 

the surviving widows of this group will be two years older. This 
would place them in a 62-66 age group. Since no 62-66 age group 

has been defined, the Department of Insurance assumes that 40 per 
cent (those aged 65 and 66) of the widows are receiving OAS in 

1970 and that 60 per cent (those aged 62, 63, 64) of the widows 

are still in the 60-64 age group. Similarly, in 1975, the surviving 
widows for the 55-59 age group of 1968 are 62-66 years old. The 


same sixty-forty treatment is applied in this case. 


WIDOW.6 uses similar assumptions, but on a much more 
widespread basis. For the year of widowhood 1968, WIDOW.6 cal- 
culates benefits for each year from 1968 on. Consider those 
widows in the 60-64 age group in 1968. In the years from 1969 to 
1975, the surviving widows of this group will be in the age groups 
61-65, 62-66, 63-67, 64-68, 65-69, 66-70, and 67-71, respectively. 
Hence, in 1969 only 80 per cent of the surviving widows will 
receive flat rate benefits, in 1970 60 per cent, in 1971 40 per 
cent, and in 1972 20 per cent. By 1973, all of the surviving 
widows of the 60-64 age group will be 65 or over and will receive 
no further flat rate benefits. But, the survivors of the 55-59 
age group of 1968 make up the 60-64 age group in 1973. Thus it 
starts all over again. In 1974, only 80 per cent of the surviving 
widows (of the 55-59 age group of 1968) will receive flat rate 


benefits and so on. 
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This is the method used by WIDOW.6. It has the advantage 
that it treats all years of widowhood in the same manner. There 


are no special cases for special years. 


A further comparison of the two methods can be obtained 
by examining one of the worksheets prepared by the Department of 
Insurance. Consider the worksheet for the year of widowhood 1968. 
On that section of the page dealing exclusively with the calculation 
of the benefit, there is a "diagonal" line running from the top 
right-hand corner to the bottom. left-hand corner. This:*line moves 
downwards in steps of five years for each five-year age group when 
moving from right to left on the page. If the "diagonal" line 
associated with the WIDOW.6 method were superimposed on the sheet, 
it would move downwards in five steps of one year each for each 


five-year age group when moving from right to left on the page. 


The earnings-related portion of the widows’ Benefits has 
two components. First, the payment to widows under the age of 65 
for any given year is dependent upon the number of widows surviving 
in that year and upon the user-supplied percentage of benefits that 
are earnings-related. Second, the payment to widows 65 years of 
age and over is calculated in a similar manner to the under 65 
payment. However, since a known proportion of widows 65.and over 
are receiving other forms of retirement pensions, their earnings- 
related benefit is reduced by a certain factor. This factor, called 
the reduction factor, is calculated in WIDOW.6 using the same 


method as the Department of Insurance. 


The discussion in earlier paragraphs pertaining to the 
differences in the two methods of calculating the flat rate benefit 
applies equally to the earnings-related benefits. Widows under the 
age of 65 receive the flat rate payment and the earnings-related 
payment. Widows 65 years of age and over receive OAS and an earn- 


ings-related payment (which may be reduced by a reduction factor). 
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The last step in the program is the actual "adding up" 
of the benefits calculated for each year of widowhood to arrive 
at the total benefits (both flat rate and earnings-trelated) paid 
out for each year in the period 1968 to 2025. The "adding up" 
procedure of the Department of Insurance requires the use of a 
high order interpolation formula with pre-determined weights. 
WIDOW.6, due to its design, requires only a simple addition of 


numbers. 


E. POLICY VARIABLE INPUT 


Most of the policy variables which influence the cal- 
culation of the widows' benefits are initially specified when 
running INDEX.2, MMAE.3, and RETBN.4. However, certain policy 
variables are unique to WIDOW.6. First of all, the user specifies 
a new proposed annual flat rate payment and the year he wishes that 
it begin. Also, the user specifies new percentages of benefits 
that are earnings-related for widows under 65 years of age and for 
widows 65 and over. Under the current plans the program assumes 
that the annual flat rate is 300 dollars and that it does not 
change. Under the proposed plan, the program assumes that the 
annual flat rate becomes 960 dollars in 1973. Under the current 
plan, the program assumes that the percentages of benefits that 
are earnings-related for widows under 65 and for widows 65 and 
over are 37.5 and 60 per cent respectively, and that they do not 
change. Under the proposed plan, these percentages both become 


re 


F. OUTPUT FILES 


Each run of this program produces one Output file. The 


name of the output file depends on the plan being tested. 


RSCUR = RSCWD 
MICUR ~ MICWD 
RSPRO - RSPWD 
MIPRO - MIPWD 


TXFPB ~ TFPWD 
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In addition, when the run is complete, a message will 
be printed on the terminal giving the "END OF RUN" statement, the 


CPU time required for the run, and tte name of the output file. 


G. TELETYPE QUESTIONS 


1) TYPE IN PLAN TO BE TESTED - 
The user types in one of the 5 possible 5-digit codes 
mentioned in chapter V depending on the plan to be 


tested. 


la) TYPE IN GROWTH RATE FOR AVERAGE EARNINGS AFTER 1975 
(e.g., FOR 5.5%, TYPE 1.055) --- 
This question appears only if TXFPB is specified above, 
and is self-explanatory. The program assumes that for 
reasonable stability, the value is 1.035, and that for 


moderate inflation the value is 1.055. 


2) TYPE IN RETIREMENT PENSION RATE (0-1) --- 
Generally, the answer to this question is 25 per cent. 


Hence, the user must type in .25 (i.e. in decimal form). 


3) TYPE IN STARTING YEAR (DEFAULT IS 1968) --- 
The user types in the year of widowhood for which cal- 
culations are to begin. The year must be in the range 


1968 to 2025 inclusive. 


4) TYPE IN PROPOSED ANNUAL FLAT RATE PAYMENT, AND THE YEAR IT 
BEGINS --- 
See chapter VIII, section E for the program's basic 


assumptions about flat rate payments. 
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The following is a list of suggested user responses 

depending on the scheme being tested: 

i) If the user is running a current plan (RSCUR or 
MICUR), he should type in 300, 1968 in response to 
this question. 

ii) If the user is running a proposed plan (RSPRO or 
MIPRO), he shouid type in 960, 1973 in response to 
this question. 

iii) If the user is running a different plan, he can 
specify other proposals. For example, under RSPRO, 
MIPRO, or TXFPB, if the user types in x,y where y> 
1973, then the program will generate flat rate pay- 
ments of 300 dollars beginning in 1968, 960 dollars 
beginning in 1973, and x dollars beginning in the 


year y. 


5) TYPE IN PERCENTAGE (BETWEEN 0 AND 1) OF BENEFIT THAT IS EARNINGS 
RELATED 
(1) FOR WIDOWS UNDER 65 --- 
(2) FOR WIDOWS 65 AND OVER --- 


BEGINNING IN THE YEAR 


See chapter VIII, section E for the programs's assumptions 


about earnings-related benefits. 


The following is a list of suggested user responses 
depending on the scheme being tested: 
i) If the user is running a current plan (RSCUR or 
MICUR), he should type in 
3/5 
6 
1968 
in response to this question. 
ii) If the user is running a proposed plan (RSPRO or 
MIPRO), he should type in 
ep 
75 
973 


in response to this question. 
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iii) If the user is running a different plan, he can 
specify other proposals. For example, under RSPRO, 
MIPRO or TXFPB, if the user types 


x 


N 


tf 
where y>1973 and O<x,z<l, then the program will generate 
earnings related payments according to the following 


percentage (expressed as decimals): 


Widows under 65 Widows 65 and Over 
Beginning in 1968 0375 .6 
1973 ee Pe de 
tf x z 


H. NOTES 


The program consists of a main program and one sub- 
routine, entitled REDFAC. The purpose of this subroutine is to 
compute the reduction factors for the calculation of the earnings- 


related benefit for widows 65 years of age and over. 


The program requires 12K of core and generally uses 
about 80 seconds of CPU time for complete execution. When 
running the program, the user, after answering all questions, may 
have to wait as long as five minutes (clock time) before receiving 
the 'END OF RUN' message. This delay is normal and its length 


depends on the workload in the computer. 


When preparing to run this program, it is recommended 


that the user make use of the Run Sheet for the program. 


The WIDOW.6 program does not produce results that are 
identical to those produced by the federal Department of Insurance. 
The reason for these differences is related directly to the different 
time intervals for which the calculations have been made (WIDOW.6 
covers a one-year interval while the Department of Insurance covers 


a five-year interval). 
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IX DISBN.7 


A. INTRODUCTION 


Program DISBN.7 calculates the disability pensions in the 
PENSIM model. It also calculates the pension paid to the husband or 
wife of a disabled contributor if the proposed scheme is being tested 
or if the user desires it as an option in some other scheme. The 
program has the facility to test different flat rate payments and 
different ratios of the imputed retirement pension for the earnings- 


related portion of the disability pension. 


B. VECTORS AND MATRICES 


POP contains the population data read in from ONE. 

DTHFAC - contains the death factors beginning in age group 20-24 
and incrementing in five-year intervals to the 90+ age 
group. These factors are a result of output produced 
by RETBN.4. 


PRVRTE - stores two sets of disability prevalence rates. The 
first set is used for the earnings-related calculation 
and was developed by the Department of Insurance. The 
second set is used for the calculation of the flat rate 
disability benefits and is found in Schedules 2 and 3 
of the 1964 Actuarial Report. 4 In each set there are 
six quinquennial year groupings and nine age groupings. 
The values for 1975 are used from 1970 to 1979. 


DISINS - contains four proportions indicating the ratio of con- 
tributions made by persons insured for disability. 
The first is for males and the remaining three are for 
females and pertain to the years 1970-1979, 1980-1989 
and 1990+ respectively. 


PROP - contains ten sets of proportions of total population 
insured for disability benefits. These are used in the 
flat rate caiculations. The first two sets are for 
maies for the years 1970-1989 and 1990+ using reason- 
able stability or moderate inflation assumptions. For 
females there are two groups of four sets each, one 
group for moderate inflation and one group for reason- 
eble stability. The four sets pertain to the year 
intervals 1970-1979, 1980-1989, 1990-2009, and 2010+ 
respectively. The other dimension of the array is for 
the nine age classes which are: 


20 - 24 45 - 49 
yy +) 50 - 54 
30 - 34 je OE 
be a 60 - 64 
40 - 44 


“Department of National Health and Welfare, The Canada Pension Plan: 
Actuarial Report, November 1964, p. 81. 
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DISFR - stores the disability flate rate calculations. 


YMPE -~ contains the appropriate YMPE's for the scheme being 
tested. 
PNIND - contains the pension index that corresponds to the 


scheme being tested. 
DISERN - stores the disability earnings-related calculations. 
PHWHWC - is an acronym for proportion of husbands and wives at 
home with children. The vector is used in the com- 
putation of wives' benefits and was developed by the 


Department of Insurance. 


FACTOR - contains the proportion of wives entitled to wives' 
benefits. 


FMTH65 - contains the number of females per 1000 married to 
husbands over age 65 who were disabled at age 65 and 
continue to meet the disability test. This data was 
developed by the Department of Insurance. 

WIFEBN - is used in the calculation of wives' benefits. 

OAS - is a vector which contains the current plan flat rate 


payment in 1968, the proposed plan flat rate payment 
in 1973 and a user-specified flat rate. 


C. INPUT FILES AND DATA STATEMENTS 


There are four data files read into DISBN.7. As in other 
programs in the PENSIM model, three of these files are keyed from the 
code word entered at the beginning of the program. These three files 
contain the death benefit factors, YMPE's and the pension indices re- 


spectively. The code words govern the following files: 


DEATH BENEFIT PENSION 
YMPE FACTORS INDEX 
RSCUR DHRSC RS PNX 
RSPRO DHRSP 
MICUR DHMIC MI PNX 
MI PRO DHMI P 
TXF PB DHTFP TFPNX 


In addition the population data is read in from file ONE. 


D. METHODOLOGY 


The initial operation of the program is common to both the 
computation of earnings-related and flat rate pensions. The first 
age group is scaled to account for those participants who are in the 


22-24 group. Subsequently, the population is multiplied by the 


- ui therauls 


; - 7 : any ae ieee _ Lae te r : i ' ; mn! 
. Ane ae - 
“UNS aie : 7 - - 7 
Dire. Ww cl Rd cag hiss Pea aac i" ; Pe pape gt i” j Kis : 
: mi 


i ke 7 

> tuts Woe! Paid wet eae ar cay 
Son 46 A oH bal Cheat. inane oO : 
jae. bcosk gall: 4 Piar ac ae bol eth bn cae | 
epee Te aaa ade ed 


mW “nara 3o daidaiueat a a2 si oui “- ie 


si : bi >i :oN ae eo ben wr SC cat wS . io ew oatcunp | a 44 _ : | Saf) - + 
due 7 ohne at? asls, Cerne, ray was ety ae digdangarn 
vt bint 4Saf*..bo99; ies sane @ ‘hoe, EVER TG 


a i eA ' Me's 16 6 wei-@s 23464: is B dial Wat ¢88 2an- 
> oe ‘o ‘24 a 69 o653 3 j 2fca = is on mi ee os a3 ace [2o34 
sesdt ugedl  .aetyrs4 o€2 MH gasantget of? Sa GoRelge Dire aes 
¥o : ~ ya “Ms SD Di , 4 » 2 | ae ee ; 26 4. i Zi p?aus. 
7 a s* » 
Ktiga SS é la soe 48 (iit) L. aeee 
' 6. SIA 
uM de De ab 7 
s = — 
AG yi Os8 
ert] Scan 
oe 
as ve vty e.3 or) 
ee i [ee 
=~ 
nau 3 ecs 


4 A ‘ P aan 
, ay | ryt. 2 2) cb ak 4 = 7 ws ii. = e¢ ’ : Tied : “o- ag ‘ 
Bs 7 ; 
= 4 a a) 
p= a 
Oe 
- ee Ad hy | - 
a 
- i _ 
vo y= 
923 = end = i ar be i - a at a via a ae 
f oy 
3:5 ad eet, 2/87 205% Gah iio. Cas) SO ELIE, on oats 
: yey 
o@3 a2! och aay é ing Fat SOG / a4 SS 59) GSS BI Geet? ae € ‘ 


= net uh Opus ait lacs ipiselioges al} @tbiee soured nee 


ey 2 


appropriate prevalence rates to determine the numbers of persons to 
whom the disability pension will apply. At this point the two calcu- 
lations diverge. The earnings-related portion of the disability benefit 
is multiplied by the death benefit factor to incorporate a lifetime 
earnings relation for the imputed retirement pension. The flat rate 
portion is multiplied by the proportion insured for disability to 
determine total recipients of the flat rate portion in every year. 

These values are then copied into another matrix for use in the wives' 


benefit calculation. 


Both benefits are then totalled for all the age groups for 
each year. The costs are determined by multiplication by the annual 
flat rate kar be the flat rate portion and by the ratio of the 
imputed retirement pension for the earnings-related portion. The 
latter is escalated by the growth rate of the three-year moving aver- 


age of the YMPE's and the former by the pension index. 


If the user desires a wives' pension, he indicates this to 
the computer in response to a teletype question. Since the number of 
wives under 65 whose husbands have been disabled has been calculated 
it is only necessary to compute the number over 65. To accomplish 
this the number of females over 65 is multiplied by the ratio in- 
dicating the number of females married to husbands over 65 who were 
disabled at age 65 and continue to meet the disability test. It is 
now necessary to determine the number of wives with disabled husbands 
who have children. By multiplying the total number of wives with 
disabled husbands by the proportion of husbands and wives at home with 
children, the number of wives eligible for the benefit is determined. 
The annual flat rate pension is multiplied by these values to cost the 


benefit. 


E. POLICY VARIABLE INPUT 


Economic assumptions are entered in the model through the 
input of the YMPE's and the pension indices. If the user is running 


a TXFPB option, the program inquires if moderate inflation or reasonable 
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stability assumptions are desired. With this information the program 
can use the appropriate proportion of the total population insured 


for disability. 


The choice of the proper YMPE and pension index keys off 


the initial code word entered. 


YMPE PENSION INDEX 
MIPRO MI PNX 
MICUR 

RSPRO RSPNX 
RSCUR 

TXFPB TFPNX 


During the courses of execution of the program the user must 
indicate the flat rate value of disability benefits and the ratio of 


the retirement pension for the earnings-related portion. 


F. OUTPUT FILES 


The program produces two output files; the first for dis- 
ability benefits and the second for wives’ benefits. The output file 
for disability benefits stores the earnings-related portion first and 
then the flat rate portion. If no wives' benefits are desired, as in 
the current plan, then this output file is filled with zeros. The out- 


put files produced are: 


YMPE CODE DISABILITY PENSION WLVES' PENSION 
RSCUR DSRSC RSCWF 
MILCUR DSMIC MI CWF 
RSPRO DSRSP RSPWF 
MI PRO DSMIP MIPWF 
TXF PB DSTFP TEPWF 


The values stored in these files are for the years 1970 


to 2025. 


G. TELETYPE QUESTIONS 


Before execution will continue the user must respond to 
the following questions: 
1) TYPE IN PLAN TO BE TESTED: 


Respond with one of the five-letter code words outlined in 
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chapter V. If the user misspells one of these words, the 
terminal will reply for example, RXPRO? TRY AGAIN, 
2) ECONOMIC ASSUMPTION - (RS/ML): 
Respond with reasonable stability (RS) or moderate inflation 
(ML). This question is asked if the TXFPB option is specified. 
3) ANNUAL FLAT RATE BENEFIT AND YEAR IT BEGINS: 
Respond with a value for the flat rate benefit and the year 
this benefit is to commence. For example, if the annual 
flat rate is to be $2000 in 1990, type in 2000. 1990. If the 
user is running a current plan option the flat rate is set at 
$300 until the year the user wishes to change it. If it is 
not to be changed, 300. 1970 should be typed in. If the user 
is running a proposed plan option the flat rate is set at 
$300 up to 1973, $960 from 1973 to the year the user wishes 
to change it. If 2000. 1980 was entered, the costing would 


use the following values in the years indicated: 


$ 300 1970 - 1972 
$ 960 1973 - 1979 
$2000 1980 - 2025 


If no change is to be made in the proposed plan, the user 
would enter 960. 1973. When a TXFPB option is specified, 
the following additional questions result. 

3a) FLAT RATE - (CURRENT/PROPOSED): 
If current is entered the flat rate payment option will 
operate as the current plan in question 3. If proposed is 
entered, the option will operate as the proposed plan. 

4) WIVES' BENEFITS (YES/NO): 
Self-explanatory response. 

5) RATIO OF IMPUTED RETIREMENT PENSION AND YEAR IT BEGINS: 
Follow the same procedure as for question 3. 

5a) EARNINGS RELATED - (CURRENT/ PROPOSED) 


Follow the same procedure as for question 3a. 
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H. NOTES 


In section G there is no question pertaining to the amount 
of the wives' benefits. They are always assumed to be identical to 
the flat rate disability pension except in the case where there is no 
wives' benefits at all. If the user wishes to alter wives' benefits 
to some other value, a programming change is necessary to DISBN.7. 

As it currently operates, DISBN.7 calculates wives' benefits on the 
basis of $960 in 1973 and escalates it to the year that a change takes 
place. If a current plan is specified, there should be no wives' 
benefits. If these benefits are specified and current option is in 


effect, the wives’ pensions will be set at $960 until it is changed. 


a ipa ie i - marie hone Bog ise 5 
ie AD cla ee mae nih an ici Ue e x aa 


as fiat a 4 ru 


ig e 


7 an 


a 


X CHILD.8 


A. INTRODUCTION 


CHILD.&8 is the eighth program in the PENSIM model. It 
calculates the benefit costs for children who have been orphaned 
by the death of their father (regardless of whether their mother 


survives) and for children whose fathers have been totally disabled. 


The program calculates these benefits under various 
user specified options. Benefits are projected for the years 
1968 to 2025 for orphans and for the years 1970 to 2025 for children 


of disabled beneficiaries. 


B. VECTORS AND MATRICES 


DUMMY - is a dummy vector which is used to skip over the 
population data that is not needed in the program. 


POP - contains the population data read in from the data 
file ONE. Only the first four age groups of each 
sex are used. They are: 


QO - «4 
5.- 9 
10 = 14 
15 - 19 
PROB - contains the probabilities of being an orphan for 


the four age groups. The probabilities are for the 
years 1968, 1970, 1975, 1980 ete. to 2025 in five- 
year intervals. These probabilities were developed 
by the Department of Insurance, Ottawa. 


PNIND - contains the appropriate pension index. 

ADJ - contains adjustment factors to scale the first five 
years of orphans benefits so that they will be in 
line with actual trends. 

SDFIFD - contains the percentage of surviving disabled fathers 
insured for disability for the first four age groups 
and for the years 1975, 1980, 1990, 1995 and 2000 plus. 
This data was developed by the Department of Insurance. 


PROBL = contains the linearly interpolated values of PROB 
for the years 1968 to 2025. 


PROB2 - contains the linearly interpolated values of SDIFID 
for the years 1970 to 2025. 


ORPHBN - is used for the calculation of orphans' benefits. 
CHLDBN - is used for the calculation of children's benefits. 
YEAR - contains the years 1968 to 2025, 


FILE - contains the input key words and output data file 
titles . 
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C. INPUT FILES AND DATA STATEMENTS 


The contents of the arrays SDFIFD and PROB are entered 


in data statements at the beginning of the program. 


The first data read in are the population statistics. 
The first four age groups for males are read in; the remaining 
statistics for males are skipped and then the first four age groups 
for females are read in. Subsequent to reading the population 


statistics, the program reads in the user specified pension index. 


D. METHODOLOGY 


The program contains two logically separate sections; 
the first calculates orphans' benefits and the second calculates 
shuttle benefits. One of the first operations common to both 
sections is a scaling of the age group 15-19. Since the CPP 
makes payments to children who are less than 18, it is necessary 


to scale this age group. The scale factor used is. .6. 


After the scaling ‘the pvogram begins to calculate 
orphans' benefits, Since the probabilities of being an orphan 
are given in five-year intervals, it is necessary to interpolate 
values for the years between the intervals. The program uses a 
linear method to interpolate these values. To calculate the 
number of orphans in any year, these probabilities are multiplied 
by the corresponding population statistics. The number of orphans 
are then totalled for each age group. To cost the benefits, the 
annual payment made to orphans is multiplied by the number of 


orphans. 


The calculation of children's benefits is very similar. 
Instead of interpolating probabilities of being an orphan, the 
percentage of surviving disabled fathers insured for disability 


are interpolated using the same method as before. 
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Economic assumptions enter the computation via the 
pension index which escalates the annual benefits paid to orphans 


and children. 


E. POLICY VARIABLE INPUT 


As in previous programs in the PENSIM model, CHILD.& 
expects one of the five code words before execution will continue. 
Unlike other programs, CHILD.8 does not read in the YMPE's. The 
only policy variable the program uses is the pension index. Two 
assumptions are implied in this index. If the options being 
tested assume moderate inflation, then the pension index will grow 
at its current two per cent maximum. In addition, the maximum 
can be altered to suit the user's options. These assumptions are 
specified in the program INDEX.2. The code words call the following 


files for the pension index: 


MIPRO 
RS PRO 
RSCUR RS PNX 
TXFPB - TFPNX 


During the course of execution of the program, the user 
must indicate the flat rate values of the orphans' and children's 


benefits so that the costs of these benefits can be calculated. 


F. OUTPUT FILES 


The program produces one output file that contains 
the caleulated values for both orphans' and children's benefits. 
In the output file, the fifty-eight values for orphans' benefits 
are output first, followed by the fifty-eight values for children's 
benefits. The program calculates values for orphans’ benefits 
from 1968 to 2025 while children's benefits are calculated from 
1970 to 2025. As a result, the first two entries of the children's 


benefits are zero. 
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The name of the output file produced is keyed off the 
original code word. 
RSCUR = RSCCH 
MICUR - MICCH 
RSPRO - RSPCH 


MIPRO =- MIPCH 
TXFPB - TFPCH 


G. TELETYPE QUESTIONS 


For execution to continue the following questions 


must be answered. 


1) TYPE IN PLAN TO BE TESTED: 
Respond with one of the five-letter code words 
outlined in chapter V. If the user misspells 
one of these words, the terminal will reply for example, 
RXPRO? , TRY AGAIN. 

2) ENTER ANNUAL VALUE OF ORPHANS BENEFIT AND YEAR IT BEGINS 
Respond with a value for orphans' benefit and 
the year this benefit is to begin. For example, 
if the user wishes the orphans’ benefit to be 
changed to $600 a year in 1980, type in 600. 1980. 
If the current level of benefits are to be main- 
tained, type in the unescalated amount in 1968 
i.e. 300. 1968. If the user enters some value for 
1980, the program automatically escalates $300 
until 1979. 

3) ENTER ANNUAL VALUE OF CHILDRENS BENEFIT AND YEAR IT BEGINS 


Pollow the same procedure as for question 2. 
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XI FUND.9 


A. INTRODUCTION 


The purpose of this program is to generate summary 
tables for the Canada Pension Plan Fund for the period 1970-2025. 
The program takes as its input, the output filés of the five 
programs described in chapters VI-X and calculates total benefits, 
expenses of administration and an annual value for the fund. The 
method used in making all cakeulations in this program is identical 


to the method used by the federal Department of Insurance. 


When running the program from a terminal, the program 
generates certain questions which the user must answer. The 
answers supplied by the user instruct the program as to which plan 
is being tested and provide useful information for the output 


tables. 


Be VECTORS AND MATRICES 


AIESPY - is a 56-element vector containing annual figures 
(1970«2025) for the interest earned since the 
preceding year on outstanding investments. 


ANNINT 


is a 60-element vector containing annual figures 
(1966-2025) for the total amount of interest 
collected since the CPP fund began in 1966. 


CHILD - is a 56-element vector containing annual figures 
(1970-2025) for children's benefits (see chapter X). 


CPPCN - is a 56-element vector containing annual figures 
(1970-2025) for contributions to the fund 
(see chapter VII). 


DHBEN < is a 56-element vector containing annual figures 
(1970-2025) for death benefits (see chapter VI). 


ERDSB - is a 56-element vector containing annual figures 
(1970-2025) for earnings-related disability 
benefits (see chapter IX). 


ERWID - is a 56-element vector containing annual figures 
(1970-2025) for earnings-related widows' benefits 
(see chapter VIII). 


FRDSB <= is a 56-element vector containing annual figures 
(1970-2025) for flat rate disability benefits 
(see chapter IX). 
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FRWID - is a 56-element vector containing annual figures 
(1970-2025) for flat rate widows' benefits 


(see chapter VIII). 


MATURE - is a 56-element vector containing annual figures 
(1970-2025) for the value of loans that are paid 
back in each year. 


ORPHAN - is a 56-element vector containing annual figures 
(1970-2025) for orphans' benefits (see chapter X). 


RATE - is a 56-element vector containing annual figures 
(1970-2025) for the interest rate charged on 
investments made in each year. 


RTBEN - is a 56-element vector containing annual figures 
(1970-2025) for retirement benefits (see chapter VI). 


TOTAL - is a 56-element vector containing annual figures 
(1970-2025) for total benefits and expenses of 
administration. 


TOTINV - is a 60-element vector containing annual figures 
(1966-2025) for the total amount available for 
investment. 

WFDSB - is a 56-element vector containing annual figures 
(1970-2025) for wives' disability benefits 
(see chapter IX). 


CNRATE - is a 60-element vector containing annual contribution 
rates (1966-2025). 


C. INPUT FILES 


The input files to FUND.9 are all dependent on the plan 
being tested. For the five types of plans available, namely 
Reasonable stability current plan - RSCUR 
Moderate inflation current plan -  MICUR 
Reasonable stability proposed plan - RSPRO 


Moderate inflation proposed plan -  MIPRO 
All other plans -  TXFPB 


The input files are: 


Contributions RSCCN MICCN RSPCN MIPCN TFPCN 
Retirement benefits RTRSC RTMIC RTRSP RIMIP RITFP 
Death benefits DBRSC DBMIC DBRSP DBMIP DBTFP 
Widows' benefits RSCWD MICWD RSPWD MIPWD TFPWD 
Disability benefits DSRSC DSMIC DSRSP DSMIP DSTFP 
Wives' benefits RSCWF MICWF RSPWF MIPWF - TFPWF 


Children's & orphans’ benefits RSCCH MICCH RSPCH MIPCH TFPCH 
Interest rate on investments RSINT MIINT RSINT MIINT TFINT 


where all files have the extension ".DAT". 
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The first four elements of the vectors ANNINT and 
TOTINV are initialized with actual values for the years 1966-1969. 


See section B for the description of these variables. 


The user is reminded that before he can run the FUND.9 
program, he must run the eight programs of chapters III - X, all 


for the same plan. 


The last data source that must be read in is file OPSHN. 
This data contains the annual total contribution rate for each year 


1966-2025. 


D. METHODOLOGY 


The first table generated by FUND.9 is a summary table 
of total benefits and expenses of administration for each year 
(1970-2025). Annual figures for each type of benefit are read in 
from input files. Expenses of administration, according to the 
Actuarial Report (1969), are defined to be one-tenth of one per 
cent of contributory earnings for each year. Contributory earnings 
are calculated by dividing total contributions by the contribution 


rate. 


The second table generated by FUND.9 contains three 
figures for each year (1970-2025), namely, total contributions, 
total benefits and expenses of administration, and the value of 
the fund. The most important calculation in this section is the 


computation of the value of the fund. 


The value of the fund in any given year is calculated 
as follows. First, the value of interest earned since the preceding 
year on outstanding investments is calculated. Next, the total 
amount available for investment in that year is computed. It is 


equal to the sum of: 


1) total contributions 


2) maturing loans 
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3) interest earned in the preceding year 
on outstanding investments 


4) less total benefits and expenses of 
administration. 


The program assumes that loans mature after twenty years. 


If the total amount available for investment is negative, 


it may be necessary to recall certain loans (see section E) in 


order to prevent the fund from going bankrupt. 


The remaining steps in the calculation of the value 


of the fund require the definition of the following program 


variables: 


TOT6 - 


COL12 - 


COL13 - 


COL14 - 


COL15 - 


is the cumulative total value of loans that have 
been paid back. This variable has a different value 
for every year after the first loan has been repaid. 


is the cumulative total value of investments that 
have been made. This variable has a different value 
for every year (1970-2025). 


is the total value of outstanding investments 
(i.e. TOT7 - TOT6). This variable has a different 
value for every year. 


is equal to COL12 plus interest earned in the first 
quarter of the next year on outstanding investments. 
This variable has a different value for every year 
and allows for the three-month time lag on investment 
interest. 


is the difference between contributions and total 
expenses (i.e. benefits and expenses of administration) 
in 1969 multiplied by a factor equivalent to the 
interest rate for 1969 raised to the power (YEAR - 1969) 
where YEAR ranges from 1970 to 2025. This variable 

has a different value for every year. 


is one-eighth of the following quantity: COL14 less 
the difference between contributions and total expenses 
in the current year. This variable has a different 
value for every year. 


is the sum of COL13 and COL15. This variable has a 
different value for every year. 


The third table generated by FUND.9 is a summary table 


designed to show the behaviour of the fund in more detail. The 


figures appearing in this table were calculated in the earlier 


steps of the program. "Total Revenue" is simply the sum of 


‘Contributions", "Interest", and "Loan Repayments", "Net Revenue” 
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is equal to "Total Revenue" less "Benefits and Expenses". Note 
that "Net Revenue" is the same as the total amount available for 


investment. 


E. POLICY VARIABLE INPUT 


Most of the information which the user must provide 
via the teletype when executing FUND.9 is, in fact, for inforn- 
ation purposes only. It is used to print a summary page describing 


in detail the plan being tested. 


It is possible that under certain plans the Canada Pension 
Plan will become exhausted at some future point in time. In order 
to help prevent this, the Canada Pension Plan empowers the Minister 
of Finance to recall loans. In FUND.9, the user has the option of 
i) calling back the most recently made loans, or 
ii) calling back no loans at all, or 


iii) calling back the loans with the shortest 
time to maturity. 


The method used by the Department of Insurance was the first of 

the above. Once the loan has been recalled, interest is no longer 
payable. The program assumes that the full loan made in any year 
will be recalled and will not redeem only a partial issue. Thus 
any excess recalled and not needed to cover expenses will be loaned 


out again as a new investment. 


Should the user indicate that he is executing the program 
for Ontario only, the fund data for 1966-1969 is taken as 55 per 
cent of the total for those years. Similar ratios could be used 


for the other provinces. 
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OUTPUT FILES 


Each run of this program produces one output file. The 


name of the output file depends on the plan being tested. 


RSCUR - #£RSCFD 
MICUR - MICFD 
RSPRO - £RSPFD 
MIPRO = £MIPFD 
TXFPB - . TFPFD 


In addition, when the run is complete, a message will 


be printed on the terminal giving the"END OF RUN" statement, the 


CPU time required for the run, and the name of the output file. 


G. TELETYPE QUESTIONS 


1) TYPE IN PLAN TO BE TESTED 


The user types in one of the appropriate 5-digit codes. 
Should the user type in some other combination, an 
error message will come up and the user can retype 


the correct version. 


la) IS TXFPB PLAN MODERATE INFLATION? (Y OR N) 


2) 


3) 


4) 


(DEFAULT IS REASONABLE STABILITY) --- 
This question will only be typed for the TXFPB options 
and the user's response will generate the proper set of 
economic assumptions for the interest rate of the 
fund. 

TYPE IN POPULATION SEGMENT USING THESE CODES: 

1 = NFLD; 2 - PEI; 3 - NS; 4 - NB; 5 = QUE; 6 = ONT; 

7 = MAN; 8 - SASK; 9 - ALTA; 10 - BC; 11 - CANADA; 

12 = CANADA LESS QUE (DEFAULT IS 12) --- 


Self-explanatory. 


TF EXPENSES EXCEED REVENUE, WHICH LOANS WILL BE RECALLED? 
(TYPE -1 FOR MOST RECENT; 0 FOR NONE; 1 FOR EARLIEST); 


The user can specify here the recall provisions to 
be usede 
TYPE IN YEAR'S BASIC EXEMPTION - 
The user can specify the YBE to be used. The standard 


response is 600. 
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RETIREMENT BENEFITS 

TYPE IN RETIREMENT BENEFIT RATIO AND PERCENTAGE OF 

LOWEST AER'S TO DROP OUT --- 
The user will specify which retirement ratio and 
drop out ratio were used in the PENSIM run and 
this information will appear on the summary page 
of the output file. 

IS THERE AN EARNINGS TEST (Y OR N) - 
The user responds with the appropriate letter 


depending on which assumptions were used. 


DISABILITY BENEFITS 
TYPE IN ANNUAL FLAT RATE (IN DOLLARS) AND YEAR IT BEGINS - 


The user types in the value that was used in DISBN.7 
and the year in which the changes were made. 


TYPE IN PERCENTAGE OF BENEFIT THAT IS EARNINGS RELATED 
AND YEAR IT BEGINS - 


The user types in the value that was used in DISBN.7 
and the year in which the changes were made. 


CHILDREN'S BENEFITS 
TYPE IN ANNUAL FLAT RATE (IN DOLLARS) AND YEAR IT BEGINS - 


The user types in the value that was used in CHILD.8 
and the year any change was made. If no changes were 
made, the initial amount is 300 in 1967. 


WIVES' BENEFITS 
TYPE IN ANNUAL FLAT RATE (IN DOLLARS) AND YEAR IT BEGINS - 


The user will type in the amount to be paid to the 
wife of a disabled pensioner if any, and the year 
in which such payments will begin. 


WIDOWS ' BENEFITS 
TYPE IN ANNUAL FLAT RATE (IN DOLLARS) AND YEAR IT BEGINS - 


The user types in the values used in WIDOW.6 and 
the year in which any changes take place. 


TYPE IN PERCENTAGE OF BENEFIT THAT IS EARNINGS RELATED AND 
YEAR IT BEGINS --- 


The user types in the values used in WIDOW.6 and 


the year in which any changes take place. 
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13) ORPHANS' BENEFITS 
TYPE IN ANNUAL FLAT RATE (IN DOLLARS) AND YEAR IT BEGINS 


The user types in the value that was used in CHILD.8 
and the year any change was made. If no changes were 


made, the initial amount is 300 in 1967. 
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XII CPPBN.10 


A. INTRODUCTION 


CPPBN.10 is the tenth and final program in the PENSIM 
model. This program projects maximum monthly benefits under the 
Canada Pension Plan for each type of pension. In addition, maximum 
monthly OAS and GIS payments are projected for single pensioners 


so that comparisons can easily be made. 


The projections are made for each year from 1973 to 
2025 and the amounts will vary given the alternatives under con- 
sideration. Thus, for any set of policy variables being tested 
actuarially, PENSIM can also generate the levels of benefits 
actually payable to recipients on an individual basis. Thus, 
the macroscopic view of the CPP fund can be shown together with 


the microscopic benefit levels. 


B. VECTORS AND MATRICES 


YMPE - this vector contains the YMPE values for 1966-2025 
and they are read in from the disk depending on 
which plan is being tested. 


RETBN - this vector contains the maximum retirement benefits 
available in each year and is calculated in the 
program. 

WIDOW - this vector contains the maximum widow's benefit 


available in each year and is calculated in the 
program by summing WIDFLT and WIDERN. 


WIDFLT - this vector contains the widow's flat rate benefit 
for each year 1973-2025. 


WIDERN - this vector contains the widow's earnings-related 
benefit in each year and is calculated as a percent- 
age of the imputed retirement pension. 


LMPSM - this vector contains the lump sum death benefit in 
any year and is set at 10 per cent of the YMPE for 
each year. 


DISABL - this vector contains the maximum disability benefit 
for each year and is the sum of DISFLT and DISERN. 


DISFLT - contains the flat rate benefit for disability for 
each year 1973-2025. 


DISERN - contains the disability earnings-related benefit in 
each year and is calculated as a percentage of the 
imputed retirement pension. 
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DISWIF - this vector contains the monthly benefit payable 
to a wife of a disabled pensioner for each 
year 1973-2025. 


CHLDBN - contains the monthly benefit paid to orphans and 
children of the disabled. 


OAS - this vector contains the maximum monthly OAS payment 
to anyone aged 65+ for every year 1973-2025. 


GIS - this vector contains the maximum monthly GIS payment 
to a single pensioner for every year 1973-2025. 


CPT - this vector contains the two sets of CPI growth 
assumptions for calculating the OAS and GIS benefits. 


MATUR =- this vector contains the factors used to show the 


maturing aspects of the plan and is found in a 
data statement. 


C. INPUT FILES AND DATA STATEMENTS 


The YMPE values are read in from the appropriate data file 
depending on the plan being tested. The vector MATUR is input through a 


data statement to allow for the maturing aspects of the Canada Pension Plan. 


Similarly, the vector CPI is read in from a data statement 
and this contains the annual increases in the consumer price index 
for each year to 2025, first for the moderate inflation assumptions 


and then for the reasonable stability assumptions. 


De. METHODOLOGY 


When execution of the program begins, the user is asked 
several questions (see section G) that specify the plan, the pension 
rates and monthly amounts. Once this information is provided, the 


program projects the benefit levels by type for each year to 2025. 


The first calculation is that of the retirement pension. 
This calculation multiplies the retirement pension rate by the 
average earnings ratio and then by the three-year moving average 
of YMPE's up to the year of retirement. This formula accounts for 


the maturing aspects of the Plan and calculates the monthly benefit. 
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Once the retirement pension has been calculated, the 
program uses the widow's ratio to calculate the widow's earnings- 
related portion and similarly for the disability earnings-related 
portion. The lump sum death benefit is calculated at this point as 


10 per cent of the year's YMPE. 


The next step is the calculation of the flat rate 
benefits. This process begins with the 1973 monthly flat rate 
benefits as specified by the user and escalates the amount annually 


by a pension index which is also input by the user. 


The flat rate benefits calculated are the widow's and 
disability flat rates, the benefits to orphans and children of the 
disabled, and the payment to the wife of a disabled pensioner, if 


this is applicable. 


Once these calculations have been done, the earnings- 
related and flat rate portions of supplementary benefits are added 


to get the total maximum widow's and disability pension. 


The final calculation is the OAS-GIS pension which 
begins with the 1972 maximum monthly amounts for single pensioners, 
and escalates these benefits annually by the increase in the 


consumer price index given the economic assumptions being tested. 


The program then writes out the annual benefit levels 


by type on disk file BENFT. 


E. POLICY VARIABLE INPUT 


The user must specify at the start of execution the 
plan he is using, the average earnings rates to be used and the 
ceiling on the pension index which is to be applied to the flat 


rate benefits. 
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Similarly, the user must input to the program the 
widow and disability earnings-related ratios and the monthly flat 
rate amounts to be used in 1973. It is at this point that the 
user can specify the benefit, if any, to the wife of a disabled 


pensioner. 


F. OUTPUT FILES 


The output from program CPPBN.10 is in file BENFT. 
This data requires a full 132 width printer because of the number 


of benefits that are being calculated. 


G. TELETYPE QUESTIONS 


1) TYPE IN PLAN TO BE TESTED: 
The user will respond with one of the five 


possibilities listed in chapter V. 


This response automatically inputs the proper data source 
of YMPE values for calculating the retirement benefit and it selects 


the proper set of economic assumptions. 


2) TYPE IN PENSIONER'S LIFETIME AVERAGE EARNINGS RATIO (0 = 1): 
Respond with a decimal value in the proper range. 

3) TYPE IN CEILING ON PENSION INDEX in 1973 ET SEQ: 
Respond with a value such as 1.02. This is the 
statutory limit of the current plan. 

4) TYPE IN RETIREMENT RATE (0 - 1): 
Respond with a decimal value in the proper range. 
The statutory value is .25. 

5) TYPE IN WIDOW'S MONTHLY FLAT RATE IN 1973: 
Respond with a monthly dollar value. For example, 
the current plan amount is 28.15 and the proposed 


value is 80. 
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6) 


8) 


9) 
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TYPE IN WIDOW'S EARNINGS-RELATED RATIO (0 = 1): 
Respond with a decimal value in the proper range. 
Currently the rate is .375 and the proposed rate 
is .75. 

TYPE IN ORPHAN'S MONTHLY FLAT RATE IN 1973: 
Respond with a dollar value. The amount is 


currently 28.15. 


TYPE IN DISABILITY MONTHLY FLAT RATE IN 1973: 
Respond with a monthly dollar value. For example, 
the current plan is 28.15 and the proposed amount 
is 80. 

TYPE IN DISABILITY EARNINGS-RELATED RATIO IN 1973 (0 = 1): 
Respond with a decimal value in the proper range. 
Currently, the rate is .75 and the proposed rate 


is 1605 
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XIII DATA SOURCES 


The tables in this chapter were used in the PENSIM 
model and encompass most of the actuarial assumptions incorporated 


in the programs. 


In all cases, the source of the data was with the 
federal Department of Insurance. This allowed compatability of 
methodology and assured a common basis for comparisons. While 
some of the data used was published in the Actuarial Report, a 
significant amount of information was taken directly from Depart- 


ment of Insurance worksheets. 


Where necessary, original data was used as the basis 
for interpolation and is so identified on the appropriate table. 
Each table is identified as to the nature of the data it contains, 
the vector/array or input file name and the program (s) in which 


the data is used. 
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I POPULATION - CANADA LESS QUEBEC, MALES - PRJM1, PRJM - POP.1 


All Figures in Thousands 


AGE/YEAR 1966 1970 1975 1980 
0-4 803.1 658.3 644.4 733.9 
5-9 S2223 836.5 672.5 660.0 

10 - 14 749.5 833.4 857.6 696.0 
15 - 19 644.9 734.9 852.7 879.93 
20 - 24 492.9 632.0 749.3 868.7 
2. = 25 430.6 487.8 646.4 765.0 
30 - 34 444.1 447.6 505.4 665.0 
35. - 39 463.1 458.5 466.1 be Ps 
40 - 44 449.1 473.9 468.2 477.3 
45 - 49 Rh de es 437.9 47 hez 466.5 
50 - 54 365%5 378.0 428.3 461.7 
55 - 59 303.9 339 64 362.0 410.7 
60 - 64 245.6 21463 Jha 29 338.5 
65 - 69 19202 | 243.7 281.9 
70 - 74 Lac 153.24 tie 19945 
75 - 79 LIQ eZ 110.9 141.8 127 
80 - 84 63.43 68.6 68.8 69.7 
85 - 89 Pe A 30.7 3245 32.6 
90 and over 8.3 oY Liv LL. 


The data files PRJM1 and PRJF1 are the 1966 columns in the male and 
female population data respectively. PRJM and PRJF contain the re- 


maining population statistics. 
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I POPULATION - CANADA LESS QUEBEC, MALES - PRJM1, PRJM - POP.1 


All Figures in Thousands 


AGE/YEAR 1985 1990 1995 2000 
0-4 843.2 875.7 670.0 875.7 
5-9 421.5 861.8 895.8 891.9 

10 - 14 686.1 799.5 891.9 928.3 

15 - 19 720.1 712.3 827.5 922.8 

20 - 24 896.7 740.2 734.4 850.9 

25:=— 29 885.3 915.7 762.0 758.3 

30 - 34 784.6 906.9 939.3 788.7 

33a> 459 686.3 806.5 930.0 964.5 

40 - 44 536.7 698.2 818.2 942.1 

45 - 49 476.5 ey is Py 695.5 814.1 

50 - 54 457.8 468.4 525, 6 682.6 

55 - 59 443.5 441.0. 451.7 507.2 

60 - 64 384.7 416.2 415.3 425.9 

65 - 69 302.9 344.8 373.5 373.3 

70 - 74 232.3 251.0 286.6 S1L1L.Z 

75 - 79 148.1 T7343 188.0 7 8 

80 - 84 80.1 92.9 109.1 118.4 

85 - 89 33.0 38.0 US 51.8 

90 and over 11.9 12.1 13.6 15.9 
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I POPULATION - CANADA LESS QUEBEC, MALES - PRJM1, PRJM - POP. 1 


Ali Figures in Thousands 


AGE/YEAR 2005 2010 2015 
0-4 917.4 975.8 1,027.1 
5-9 898.9 941.8 1,001.3 

10 - 14 926.2 $35,0 979.8 

15 = 15 960.1 959.6 970.3 

20 - 24 946.3 985.6 986.8 

Rats Me ae 875.5 9/1.9 1,012.6 

30 - 34 786.7 904.8 1,002.4 

35 - 39 816.8 816.6 935.6 

40 - 44 977.3 832.4 833.4 

45 - 49 936.3 971.4 829.4 

50 - 54 798,35 YLEG 952.0 

55 - 59 658.0 768.8 883.0 

60 - 64 478.3 618.5 721.8 

65 - 69 383.5 430.3 555.4 

Wim 7h S1l3 SLY 28 358.9 

7TH: = 79 234.0 234.1 240.5 

80 - 84 135.9 147.7 147.7 

85 - 89 56e2 64.5 pil Pes 

90 and over 18.5 20.3 ye | 


2020 
1,064.6 
1,054.0 
1,041.1 
1,016.9 
999.0 
1,015.4 
1,044.6 
1,034.3 
952.5 
831.0 
814.2 
916.3 
828.2 
647.4 
463.2 
269.9 
151.7 
70.1 


25.3 


2025 
1,098.4 
1,093.0 
1,095.9 
1,080.1 
1,047.2 
1,029.3 
1,049.4 
1,078.3 
1,051.5 
948.6 
816.1 
785.0 
859.5 
742.4 
539.9 
348.4 
170.3 


72.0 


continued .... 


i. 060 sf 


~~ « » * 
a) 4 Peg +L 


¢.ai¢o, 
‘ 
at) 
C028 
4 
| 


; 6 anoae 


ee chal un O. Ts 
i ve 


BBE! 


ata 2 
e620, D 


JSUT. og 


By aan 
ie 
b. R00; 5 
i 
a,eL? 


ad: 


'hPse 


a” 
i guatt 
dee 
buat 


ea ae 


ALLE 
y eee 


o80t 


‘ee i bs 
yal F 


mh a 


5 
} 


wey 
AJ 
ee 


co; 


oe 
ABH -o- 


idly’ 


CA 


e dee 


paneer 


I POPULATION - CANADA LESS QUEBEC, FEMALES - PRJF1, PRJF - POP.1 


All Figures in Thousands 


AGE/YEAR 1966 1970 1975 1980 
0-4 761.8 626.3 613.5 717.8 
5-9 792.7 793.6 639.3 628.3 

10 - 14 715.8 799.7 813.6 661.3 
15 - 19 626.5 702.8 820.9 836.3 
20 - 24 494.3 625.2 i265 842.8 
25 - 29 430.5 491.7 645.9 745.6 
30 - 34 425.4 446.9 Jee 668.2 
35 - 39 446.8 439.0 465.8 532.5 
40 - 44 453.3 464.0 450.0 478.2 
45 - 49 396.3 443.9 464.8 451.9 
50 - 54 3952 384.2 441.1 462.6 
55 - 59 292.0 336.0 379.4 435.4 
60 - 64 242.7 aiose 325.9 368.4 
65 - 69 205.7 221.8 258.2 306.2 
70 - 74 174.8 181.5 197.3 230.6 
75 - 79 126.0 L377 148.7 160.6 
80 - 84 77.3 86.7 96.8 105.0 
85 - 89 35.4 40.7 Nive 53.2 
90 and over 13,2 15.0 17.4 20.3 
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I POPULATION - CANADA LESS OUEBEC, FEMALES-PRJF1, PRJF - POP.1 


AGE/YEAR 
0-4 
5-9 

10 - 14 
15 - 19 
20 - 24 
Zone 29 
30 - 34 
i bs a 
40 - 44 
45 - 49 
50 - 54 
55 - 59 
60 - 64 
65 - 69 
70 - 74 
fous 19 
80 - 84 
85 - 89 
90 and over 


1985 
802.3 
733.6 
652.8 
686.9 
860.8 
867.5 
769.6 
689.9 
545.7 
480.1 
450.2 
457.1 
423.3 
347.1 
273.7 
189.1 
113.7 

58.1 


Lave 


All Figures in Thousands 


1990 
833.0 
819.2 
759.9 
680.6 
713.9 
887.4 
893.1 
792.7 
703.1 
548.2 
478.8 
445.9 
444.7 
399.3 
310.7 


225.1 


1995 
825.0 
851.2 
847.5 
789.7 
709.7 
743.0 
915.1 
917.7 
806.5 
704.3 
546.3 
474.2 
434.9 
419.8 
358.1 
255.9 
160.2 

74.9 


28.2 
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I POPULATION - CANADA LESS QUEBEC, FEMALES - PRJF1, PRJF - POP.1 


All Figures in Thousands 


AGE/YEAR 2005 2010 2015 2020 2025 
0-4 872.1 927.5 976.5 1,012.0 1,044.2 
5-9 853.4 894.2 950.7 1,000.9 1,037.8 

10 - 14 877.2 887.3 929.7 988.1 1,040.1 

15 - 19 915.4 912.6 924.6 968.8 1,029.2 

20 - 24 912.3 949.9 949.0 962.9 1,009.0 

25 - 29 853.9 946.8 986.2 987.2 1,003.1 

30 - 34 iiane 887.5 982.1 1,023.3 1,026.4 

35 - 39 802.6 803.9 919.4 1,015.4 1,058.3 

40 - 44 957.2 819.9 822.5 938.3 1,034.9 

45 - 49 932.3 957.5 822.1 825.2 940.4 

50 - 54 803.1 926.5 951.8 818.8 822.3 

55 - 59 692.9 792.7 913.7 938.9 809.3 

60 - 64 527.8 674.4 770.9 887.8 912.6 

65 - 69 437.8 498.7 636.7 727.4 837.3 

70 - 74 368.8 392.7 447.3 5712 652.4 

75 - 79 310.8 304.3 324.0 369.1 471.2 

80 - 84 210.1 221.3 216.7 230.8 262.8 

85 - 89 101.9 117.4 123.7 tet.1 128.9 


90 and over ie eA Eas 51.8 55.5 =o PP 
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AGE GROUP IN 1965 


AGE GROUP 
18 - 19 
20 - 24 
25 - 29 
30 - 34 
35 = 39 
40 - 44 
45 - 49 
50 - 54 
55 = 59 
60 - 64 
65 - 69 


to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 


= "oy = 


II ADJUSTMENT FACTORS - ADJFC - MMAE.3 


MALES FEMALES 
Le L359L 1.04537 
91515 1.00941 
99412 1.04689 
1.00023 1.03419 
1.00247 1.01398 
1.00067 1.00034 
99521 99267 
. 98426 .98501 
96821 94721 
95424 Phe BF de 
- 93388 . 98098 


PEL DISABILITY DROP-OUT FACTORS - DSDRP - RETBN.4 


DEATH OR RETIREMENT AT AGE: 


25 30 354 40 
-46 1.002 1.003 1.004 1.005 
-41 1.002 pa eleks) 1.004 1.005 
-36 1.002 1.003 1.004 1.005 
-31 1.002 1.003 1.004 1.005 
-26 1.002 1.003 1.004 1.005 
-21 1.002 1.003 1.004 1.005 
-16 1.002 jp WAL 1.004 1.005 
-11 1.002 1.003 1.004 1.005 
-6 1.002 1.003 1.004 1.005 
«1 1.002 1.003 1.004 1.005 
4 1.002 1.003 1.004 1.005 

9 1.002 1.003 1.004 1.005 
14 1.002 1.003 1.004 1.005 
19 1.002 1.003 1.004 1.005 
24 1.000 1.002 1.003 1.004 
29 1.000 1.000 1.003 1.004 
34 1.000 1.000 1.000 1.003 
39 1.000 1.000 1.000 1.000 
44 1.000 1.000 1.000 1.000 
49 1.000 1.000 1.000 1.000 
54 1.000 1.000 1.000 1.000 
59 1.000 1.000 1.000 1.000 
64 1.000 1.000 1.000 1.000 
69 1.000 1.000 1.000 1.000 


* Interpolated from adjoining columns 
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III DISABILITY DROP-OUT FACTORS - DSDRP - RETBEN.4 


AGE GROUP IN 1965 DEATH OR RETIREMENT AT AGE: 
45% 50 558 60 65 
-50 ta -46 1.006 1.007 1.0085 1.010 1.015 
-45 to -41 1.006 1.007 1.0085 1.010 pat 
-40 to -36 1.006 1.007 1.0085 1.010 1; Ua 
-35 to -31 1.006 1.007 1.0085 1.010 L635 
-30 to -26 1.006 1.007 1.0085 1.010 1.015 
-25 to -21 1.006 1.007 1.0085 1.010 1.015 
-20 to -16 1.006 1.007 1.0085 1.010 Nae ee) 
-15 to -1l 1.006 1.007 1.0085 1.010 gH i PO 
-10 to -6 1.006 1.007 1.0085 1.010 1.015 
-5 to -l 1.006 1.007 1.0085 1.010 +6915 
0 to 4 1.006 1.007 1.0085 1.010 LL 
5 to 9 1.006 1.007 1.0085 1.010 0G 
10 to 14 1.006 1.007 1.0085 1.010 7,013 
Pte LF 1.006 1.007 1.0085 1.010 1.013 
20 to 24 1.0055 1.007 1.0085 1.010 BiG 
25 to 29 1.0055 1.007 1.0085 1,019 1.015 
30 to 34 1.0045 1.006 1.0085 LOLs 5 Is Be 
a5 a SF 1.003 1.006 1.0085 Fes 1.016 
40 to 44 1.000 1.004 1.0075 1.011 1.017 
45 to 49 1.000 1.000 1.0050 1.010 1.017 
50 to 54 1.000 1.000 1.0000 1.007 1.012 
55 to 59 1.000 1.000 1.0000 1.000 1.012 
60 to 64 1.000 1.000 1.0000 1.000 2012 
65 to 69 1.000 1.000 1.0000 1.000 1.000 


* Interpolated from adjoining columns 


Iv SHORT RANGE RETIREMENT BENEFIT FACTORS —- RFAC - RETBN.4 


YEAR DEATH AT AGE 60 DEATH AT AGE 65 
Males Females Males Females 
1970 61 21 167 ok 
1971 71 Eu 148 55 
1972 86 30 232 45 
1973 114 20 189 73 


1974 itd 42 309 61 
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V DEATH BENEFIT COST REDUCTION FACTORS - RSDRD - RETBN.4 


(REASONABLE STABILITY) 


YEAR* 


1970 
1975 
1980 
1990 
2000 
2001+ 


i 


.098 
082 
.070 
-055 
.048 
045 


* Values for intermediate years are interpolated linearly 


VI DEATH BENEFIT COST REDUCTION FACTORS - MIDRD - RETBN.4 


(MODERATE INFLATION) 


YEAR* 


1970 
Ay pe, 
1980 
1990 
2000 
2001+ 


098 
-082 
.073 
.059 
[052 
052 


* Values for intermediate years are interpolated linearly 


AGE 


AGE GROUP 1966 - 1969 
M F 
- 24 .9984 .9995 
= 29 .9985 .9993 
~- 34 .9984 .9991 
= 39 .9979 .9987 
= &4 .9970 .9981 
~ oe .9942 .9969 
= 54 .9906 .9951 
- 59 .9848 .9928 
- 64 .9768 .9881 
- 69 .9648 .9805 
ey f° .9475 .9673 
- 79 9219 .9437 
- 84 .8813 .9059 
- 89 .8232 .8484 
+ .7274 .7486 


VII ONE YEAR SURVIVAL PROBABILITIES - SURVM - RETBN.4; WIDOW.6 


.9984 
.9985 
9984 
2979 
.9970 
.9943 
9907 
. 9850 
.9770 
. 9650 
9479 
Gated 
.8815 
.8232 
lake 


M 


1970 - 1974 


FE 


7990 
«9993 
9991 
.9987 
9981 
.9970 
9952 
.9930 
. 9883 
.9809 
9676 
9443 
.9067 
.8490 
.7482 


1975 - 1979 


M 


9985 
.9986 
9985 
.9980 
«9971 
~9945 
.9912 
.9856 
.9778 
9658 
9490 
9239 
.8825 
8231 
.7267 


F 


9995 
.9994 
9991 
.9987 
9981 
9971 
.9954 
~9932 
9889 
.9816 
. 9680 
~9455 
.9079 
.8502 
7477 
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ONE YEAR SURVIVAL PROBABILITIES - SURVM - RETBN.4; WIDOW.6 


GROUP 


AGE 


1985 - 1989 


M 


9987 
.9987 
.9986 
9981 
9912 
29951 
~9921 
.9870 
9794 
9682 
es 
S921 
8841 
~8229 
OF 


Ry Bier 


GROUP 1980 - 1984 
M OF 

wn Dh .9986 .9996 
- 29 .9986 .9994 
ak .9985 .9992 
- 39 .9980 .9988 
- 44 .9972 .9982 
- 49 .9948 .9972 
- 54 .9917 .9956 
- 59 .9863 .9935 
- 64 .9786 .9893 
- 69 .9669 .9821 
= Fh .9503 .9686 
= 79 9251 .9464 
- 84 8834 .9089 
- 89 8230 .8514 
+ .7260 .7472 
a9 
2 
OA 
- 39 
- 44 
- 49 
- 54 
- 59 
- 64 
- 69 
- 74 
- 7 
- 84 
- 89 
os 


TES 


.9988 
.9988 
19987 
~99S2 
29973 
EPP 
,9926 
. 9878 
.9805 
9700 
» 9538 
“92I2 
8848 
One 
.7243 


1999 


| 


59996 
Ra de bes 
9992 
.9988 
990. 
.9974 
. 9960 
<9939 
.9900 
.9829 
~9702 
9477 
.9107 
8538 
7461 


F 


9996 
9994 
9992 
.9988 
9982 
.9973 
9958 
9937 
9896 
9824 
9693 
9471 
9097 
8525 
7467 


1990 - 1994 


M 


9987 
9987 


.9986 
9981 
-9973 
~9953 
.9924 
9875 
9801 
9693 
-9530 
.9268 
. 8846 
. 8228 
~7247 


= 


9988 
9988 
.9987 
«9982 
.9974 
99565 
-9927 
. 9879 
.9807 
. 9703 
9542 
9276 
.8850 
8227 
7241 


2000 - 


2025 


| 


.9996 
“9995 
-999e 
.9988 
.9983 
9974 
9961 
9940 
9901 
. 9830 
9704 
9478 
~9109 
.8540 
. 7460 


VIII EARNINGS TEST ADJUSTMENT FACTORS - ETADJ - CONT. 5 


YEAR 


1973 
1974 
1975 
1976 
1977 
1978 
1979 
1980+ 


MALES AND FEMALES 


6933 
.867 
. 800 
. 680 
.560 
440 
. 320 
. 200 


F 


.9996 
«9995 
9992 
.9988 
9982 
9973 
.9959 
.9938 
.9898 
9827 
9698 
9475 
«9103 
I32 
. 7463 
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IX LABOUR FORCE PARTICIPATION RATES —- PART - RETBN.4, CONT.5 


MALES 
AGE GROUP 1966 - 1984 1985+ 
18-19 ~65 -68 
20-24 1.00 1.00 
25-34 1.00 1.00 
35-49 93 94 
50-59 87 .88 
60-64 78 «/9 
65-69 oe 54 


FEMALES 
(Moderate Inflation) 


AGE GROUP 1966 - 1974 1975 - 1984 1985 - 2004 2005+ 


i 


18-19 46 50 55 58 

20-24 . 66 69 73 75 

Pes ar 44 50 oH i 58 

30-59 - 36 42 48 49 

60-64 me de 29 34 te 

65-69 ~13 13 14 14 
FEMALES 


(Reasonable Stability) 


AGE GROUP 1966 - 1974 1975 - 1984 1985 - 2004 2005+ 


ee eee eee iO 


18-19 - 46 49 53 ae 
20-24 66 68 71 73 
25-29 44 50 56 57 
30-59 36 Al 47 48 
60-64 ie) .29 33 34 
65-69 13 ~i3 14 14 


X OVERPAYMENT ADJUSTMENT FACTORS - ADJ - CONT.5 


YEAR MALES AND FEMALES 
1966-1970 . 100 
1971 .096 
1972 092 
1973 .083 
1974 074 
1975 .065 
1976 .056 
1977 052 
1978 .048 
1979 044 


1980-2025 .040 
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XI PROJECTED SURVIVAL RATIOS FOR QUINQUENNIAL AGE GROUPS - SURVP - WIDOW.6 


YEAR 
1968-70 


1970-75 
1975-80 
1980-85 
1985-90 
1990-95 


1995-2000 


YEAR 
1965-70 


1970-75 
1975-80 
1980-85 
1985-90 
1990-95 


1995-2000 


* These Age Groups differ 


15-19 
to 
20-24 


997 
997 
997 
998 
998 
998 


998 


SURVIVAL AGE GROUPS 


22-26 27-31 
to to 
27-31 32-36 
«999 998 
997 995 
997 .996 
998 996 
998 996 
998 996 
998 996 


SURVIVAL AGE GROUPS 


20-24 25-29 
to to 
25-29 30-34 
997 996 
297 996 
997 997 
998 997 
998 Fh i 
998 997 
998 997 


slightly for 1968-70 


continued .... 


32-36 
to 
37-41 


998 
«993 
994 
994 
~994 
~994 


994 


30-34 
to 
35~39 


994 
994 
994 
2995 
995 
~995 


«995 
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XI PROJECTED SURVIVAL RATIOS FOR QUINQUENNIAL AGE GROUPS - SURVP - WIDOW.6 


SURVIVAL AGE GROUPS 


37-41 42-46 47-51 52-56 
to to to to 

_YEAR 42-46 47-51 52-56 57-61 
1968-70 .997 996 992 989 
1970-75 -990 986 .979 968 
1975-80 991 987 980 969 
1980-85 991 987 980 -969 
1985-90 992 988 981 970 
1990-95 992 988 981 971 
1995-2000 992 988 982 [Siz 


SURVIVAL AGE GROUPS 


35-39 40-44 45-49 50-54 
to to to to 

YEAR 40-44 45-49 50-54 55-59 
1965-70 992 989 981 eis 
1970-75 - 992 - 988 982 974 
1975-80 2993 989 983 mi 
1980-85 -993 ~ 989 983 0975 
1985-90 993 990 . 984 976 
1990-95... 993 -990 984 976 
1995-2000 993 -990 985 977 
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XI PROJECTED SURVIVAL RATIOS FOR _QUINQUENNIAL AGE GROUPS - SURVP - WIDOW.6 


SURVIVAL AGE GROUPS 


caccesleaeeicisiiiiininithinstctt catalina tances EI EEL 


57-61 62-66 - 67-71 72-76 

to to ta to 
1968-70 982 971 952 .920 
1970-75 947 .915 . 860 FFL 
1975-80 ; 949 «917 . 863 774 
1980-85 «951 919 . 864 . 776 
1985-90 952 920 . 866 mee 
1990-95 » 953 «922 867 .779 
1995-2000 954 » 922 . 868 0/19 


SURVIVAL AGE GROUPS 


calc cce i n  e eCN  ne 


55-59 60-64 65-69 70-74 
me) to to to 

YEAR 60-64 65-69 10-74 75-79 
1965-70 956 928 885 »8izZ 
1970-75 958 931 890 .816 
1975-80 “960 oo 893 818 
1980-85 961 . 936 894 . 820 
1985-90 - 962 937 895 822 
1990-95 963 938 897 ~823 
1995-2000 964 939 897 824 
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XI PROJECTED SURVIVAL RATIOS FOR QUINQUENNIAL AGE GROUPS - SURVP - WIDOW.6 


SURVIVAL AGE GROUPS 


77-81 82-86 87-91 92+ 
to to to to 
_YEAR _ 82-86 87-91 92-96 97+ 
1968-70 866 782 668 503 
1970-75 -641 ~475 -000 -000 
1975-80 645 478 .000 .000 
1980-85 -647 480 299 .000 
1985-90 .648 -482 299 ‘hae 
1990-95 -650 -483 299 Pe le 
1995-2000 -650 484 ~ 299 bce 


SURVIVAL AGE GROUPS 


75-79 80-84 85-89 90+ 
to to to to 

mie 80-84 85-89 90-94 95+ 
1965-70 698 541 365 179 
1970-75 704 547 367 6477 
1975-80 707 ‘55. 368 .176 
1980-85 709 554 369 .176 
1985-90 710 556 2371 “4177 
1990-95 712 “557 <3i2 177 


1995-2000 712 Be Po «373 177 
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XII PERCENTAGE OF WIDOW BENEFICIARIES NOT IN RECEIPT OF AGE RETIREMENT PENSIONS 


BASED ON AGE IN 1968 - PWNRP1 - WIDOW.6 


AGE GROUP PERCENTAGE 
Less than 20 55 

20 - 24 83 

zo 29 78 

30 - 34 73 

35 - 39 67 

40 - 44 61 

45 - 49 54 . 
50 - 54 47 

55 - 59 40 

60 - 64 Hs 

65 - 69 15 


BASED ON AGE IN 1965 -— PWNRP2 - WIDOW.6 


AGE GROUP PERCENTAGE 
Less than 20 85 
20 - 24 50 
25 - 29 a5 
30 - 34 69 
35 - 39 63 
40 - 44 57 
45 - 49 50 
50 - 54 43 
55 - 59 36 
60 - 64 PA 


ae 


XIII AGE DISTRIBUTION OF WIVES - AGEDIS - WIDOW.6 


HUSBAND'S WIFE'S % DIST. HUSBAND'S WIFE'S % DIST. 
AGE GROUP AGE GROUP OF WIVES AGE GROUP AGE GROUP OF WIVES 
20 - 24 - 20 22 50 - 54 30 - 34 2 
20 - 24 69 35. = 39 5 
25 - 29 9 40 - 44 17 
30 - 34 0 45 - 49 38 
35 - 39 0 50 - 54 32 
40 - 44 0 55.- 99 6 
25 = 29 - 20 3 ob et 35 - 39 3 
20 - 24 40 40 - 44 6 
25 - 29 49 45 - 49 18 
30 - 34 8 50 - 54 37 
35 - 39 0 55 - 59 30 
40 - 44 0 60 - 64 6 
30 - 34 - 20 0 60 - 64 40 - 44 3 
20 = 24 9 45 - 49 7 
25 = 29 39 50 - 54 19 
30 - 34 43 55 = 59 35 
35 - 39 9 60 - 64 29 
40 - 44 0 65 - 69 7 
a= 3 20 - 24 2 65 - 69 45 - 49 5 
25 - 29 10 50 - 54 7 
30 - 34 38 55 - 59 Ly 
35 - 39 42 60 - 64 34 
40 - 44 8 65 - 69 28 
45 - 49 0 70 - 74 7 
40 - 44 25 - 29 3 70 - 74 50 - 54 6 
30 - 34 13 55 = 59 7 
a5 = 3? 40 60 - 64 20 
40 - 44 36 65 - 69 34 
45 - 49 8 70 - 74 27 
50 - 54 0 75 - 79 6 
45 - 49 25 - 29 1 
30 - 34 4 
35 - 39 16 
40 - 44 38 
45 - 49 34 
50 - 54 7 
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73 = 79 


80 - 84 


85 - 89 
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XIII AGE DISTRIBUTION OF WIVES - AGEDIS - WIDOW.6 


WIFE'S AGE 
__ GROUP 
55 - 59 
60 - 64 
65 - 69 
70 - 74 
75 - 79 
80 - 84 
60 - 64 
65 - 69 
70 - 74 
75 = 79 
80 - 84 
85 - 89 
65 - 69 
70 - 74 
75 Hd 
80 - 84 
85 - 89 
90 - 

70 - 74 
75 = 79 
80 - 84 
85 - 89 
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XIV PROPORTIONS MARRIED - PRPMRD - WIDOW.6 


20 = 24 30 30 30 30 30 
25 = 29 70 70 70 70 70 
30 -°34 82 82 82 82 82 
ee ae 86 86 86 86 86 
40 - 44 88 88 88 88 88 
45 - 49 88 88 88 88 89 
50 - 54 87 87 88 88 88 
33 = 59 84 85 85 86 86 
60 - 64 82 83 84 85 85 
65 - 69 80 81 82 83 84 
70 - 74 74 76 78 80 81 
75 - 79 67 70 72 74 76 
80 - 84 me 58 60 62 64 
85 - 90 43 46 48 50 52 
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XV PROBABILITIES FOR RE-MARRIAGE OF WIDOWS - REMAR.- WIDOW. 6 


For years 6f widowhood 1968-69 


Under 20° 20-24 25-29 30-34 35-39 40-44 45-49 50-54 55-59 60-64 


1968 . 688 . 789 ~869 ~923 957 e915 986 993 997 .995 


1970 -270 - 430 613 .765 . 864 924 961 981 992 1.0 
1975 541 676 795 . 880 934 967 984 6993 71.0 1.0 
1980 - 761 - 838 «902 943 970 986 6994 1.0 1.0 1.0 
1985 841 .903 945 ~970 986 994 1.0 1.0 1.0 1.0 


For years of widowhood after 1970 


Under 20 20-24 25-29 30-34 35-39 40-44 45-49 50-54 55-59 60-64 


a_————__ men ————— TT rE TT lr — Cl —— TT lC(C  Lhr———— ee ——E EES 


1970 . 288 451 627 MES, . 868 924 . 960 .980 4991 =..996 
1975 417 569 aha 837 .908 «952 975 989 3997 Tie 
1980 .700 794 872 924 960 981 oI92 -997 1.0 1.0 
1985 .813 . 880 pete 961 981 992 997 1.0 1.0 1.0 


1990 . 880 ~931 961 981 e992 .997 2,0 1.0 1.0 1.0 
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XVI PERCENTAGE OF WIDOWS ENTITLED TO WIDOW BENEFITS - PWEWB - WIDOW. 6 


YEAR OF UNDER 
WIDOWHOOD 50 50-54 55-59 60-64 65-69 70-74 75-79 80-84 85-89 90+ 


wt. Chine Lontn ae lle leer EE ES SU 


1968 90 85 75 70 60 55 0 0 0 0 
1970 20 90 80 60 45 30 20 0 0 0 
Ato 95 90 85 %5 60 45 30 20 0 0 
1980 95 I2 90 85 75 60 45 30 20 0 
1990 95 95 » fe 90 85 80 75 60 45 30 
2000 95 95 be 95 95 90 85 80 75 60 
2015 95 95 95 95 a0 95 90 90 85 80 
2030 95 Je 95 95 72 be 95 99 90 90 


Annual values for this variable were obtained by linear interpolation. 
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XVII UNESCALATED FEMALE DEATH BENEFIT FACTORS - FFAC - WIDOW. 6 


1968 - 72 0 
1973 23 
1974 oe 
1975 41 
1976 64 
1977 88 
1978 pa! 
1979 135 
1980 158 
1981 181 
1982 204 
1983 228 
1984 ve BN 


1985 - 2025 dependent upon plan being tested 
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XVIII UNESCALATED MALE DEATH BENEFIT FACTORS IN 1968 - FAC68 -— WIDOW.6 


AGE GROUP OF HUSBAND 


20 - 24 600 
25 - 29 793 
30 - 34 881 
35 - 39 864 
40 - 44 855 
45 - 49 855 
50 - 54 832 
55 - 59 785 
60 - 64 699 


65 - 69 405 


wy 
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XIX FLAT RATE PREVALENCE RATES - PRVRTE - DISBN.7 


AGE/YEAR 1975 1980 1985 1990 1995 2000+ 
22-24 0.04 0.05 0.05 0.05 0.05 0.05 
25-29 0.10 0.13 0.13 0.13 0.13 0.13 
30-34 0.24 0.35 0.37 0.36 0.36 0.37 
35~39 0.40 0.60 0.67 0.69 0.68 0.68 
40-44 0.66 0.96 1.07 1.12 1.13 1.12 
45-49 0.96 1539 1.55 1.63 1.67 1.68 
50-54 1.78 2.61 2693 3.07 x ae 3.16 
35-59 2.88 4.27 4.78 4.97 5.08 oe 4 
60-64 5.00 7.64 8.66 9.01 9.21 9.31 
XX EARNINGS RELATED PREVALENCE RATES —- PRVRTE - DISBN.7 
AGE/YEAR 979 1980 1985 1990 1995 2000+ 
22-24 0.05 0.05 0.05 0.05 0.05 0.05 
25-29 0.12 0.14 0.14 0.14 0.14 0.14 
30-34 0.30 0.36 0.38 0.38 0.38 0.38 
35-39 0.50 0.62 0.68 0.70 0.70 0.70 
40-44 0.82 1.00 1.08 Lake pS 1 ae Be 
45-49 Lea? 1.45 eat 1.64 1.6/ 1.69 
50-54 z.19 2.3 2.96 3.08 S615 3.18 
55-59 3.54 4.46 4.82 4.98 5.08 oe 
60-64 6.14 8.09 8.74 9.02 Fak 9.20 


XXI NUMBER OF FEMALES PER 1000 MARRIED TO HUSBANDS OVER 65 
WHO WERE DISABLED AT 65 AND CONTINUE TO MEET THE 
DISABILITY TEST - FMTH65 - DISBN.7 


AGE/YEAR 1975 1980 1985 1990 1995 2000 2005 2010 2015 2020+ 
45-49 0.490 0.731 .963 1.040 1.087 1.109 1.120 1.120 1.120 1.120 
50-54 1.365 2.421 3.257 3.655 3.844 3.942 3.991 4.000 4.000 4.000 


55-59 3.605 6.532 9.124 10.393 11.064 11.387 11.552 11.593 11.600 11.600 
60-64 7.350 14.112 19.846 23.189 24.881 25.790 26.223 26.360 26.392 26.400 


XXII PROPORTION INSURED FOR DISABILITY - PROP - DISBN.7 


MALES 
AGE/YEAR 1970 — 1985 1990+ 
22-24 40 40 
25-39 95 95 
40-54 90 91 
55-64 85 86 


continued ... 
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XXII PROPORTION INSURED FOR DISABILITY - PROP — DISBN.7 


AGE/YEAR 


22-24 
25~29 
30-34 
35-59 
60-64 


AGE/YEAR 


22-24 
25-29 
30-34 
35-59 
60-64 


1970 - 1975 


. 20 
»45 


1970 = 1975 


FEMALES 
(Moderate Inflation) 


1980 - 1985 1990 - 2005 2010+ 
s2l sZe .24 
ae F . 60 ~61 
47 54 Ae fo 
42 49 - 50 
- 36 43 44 

FEMALES 


(Reasonable Stability) 


1980 - 1985 1990 - 2005 2010+ 
21 bd 23 
be F oo? 60 
47 re 54 
41 .48 49 
36 42 43 


XXIII PROPORTION OF HUSBAND AND WIVES AT HOME WITH CHILDREN - PHWHWC 


AGE GROUP 


under 25 
25-29 
30-34 
35-39 
40-44 
45-49 
50-54 
55-59 
60-64 
65 + 


576 
. 738 
942 
947 
894 
843 
741 
558 
-356 
114 


XXIV PROPORTION OF HUSBANDS TO TOTAL POPULATION - FACTOR - DISBN.7 


AGE GROUP 


22-24 
25-29 
30-34 
35-39 
40-44 
45-49 
50-54 
35-59 
60-64 
65. + 


576 
. 738 
942 
947 
894 
. 843 
741 
558 
- 356 
.114 


- DISBN.7 
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XXV PROPORTION OF CONTRIBUTIONS MADE BY PERSONS 
INSURED FOR DISABILITY - DISINS - DISBN.7 


MALES FEMALES 
1970 = 2025 1970 - 1975 1980 - 1985 1990+ 
°95 .70 /5 . 80 


XXVI PROBABILITIES OF BEING AN ORPHAN — PROB - CHILD.8 


YEAR AGE GROUP 
0-4 2 = 9 10 - 14 be Yj 

1968 - 105 ~ L535 234 .324 
1970 ~425 -670 984 1.383 
Bi 328 1.710 2.589 3.577 
1980 soli 1.773 3.563 4,947 
1985 - 486 } ay 8 | Seat h 5/17 
1990 - 468 1.648 3.427 3.322 
Lg9D aor 1.574 3+296 Jaav 
2000 423 1,502 3.156 5.097 
2005 - 403 1.424 3.023 4.878 
2010 - 400 1.380 2.888 4.672 
2015 BR 1.378 2.196 4.466 
2020 he fe hears 2.805 4.376 
2025 401 nA FP 2.801 4.367 


XXVII SURVIVING DISABLED FATHERS INSURED FOR DISABILITY - SDFIFD - CHILD.8 


YEAR AGE GROUP 

eet Mai iae ee ee 0 oa ATR TA ARs ee ke 
1975 399 ~655 - 986 1.301 
1980 .488 eoLL Lene 1.636 
1985 Ba fe sore Lesse 1.748 
1990 oa. .909 1.392 1.853 
1995 - 546 .923 1.420 1.895 
2000+ .548 929 1.432 1.916 
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XIV SAMPLE OUTPUT DATA 


As described in chapter XI, PENSIM projects the fund 
levels for each year to 2025. Because of some slight differences 
in contribution and benefit results, described in chapters VI - X, 
the fund projections do not exactly duplicate the federal est- 


imates. 


The data shown in this chapter is the output for the 
Canada Pension Plan operating under the moderate inflation 
economic assumptions with the white paper proposals. However, the 
main tables in the Actuarial Report are based on a non-escalating 
O.A.S. and G.1I.S. while PENSIM assumes these payments will increase 
with the cost of living. Thus for comparison purposes, the PENSIM 
results in this chapter must be read with the capsule federal pro- 
jections on page 52 of the Actuarial Report which does allow for 


such escalation. 


The PENSIM output format is similar to that of the 
Actuarial Report where Table 1 in this chapter duplicates the 
style of Table 1B in the Report and Table 2 is similar to Table 
3 in the Report. As explained in chapter XI, the third table in 
this chapter gives disaggregated information for revenues and 


costs of the plan. 


The following table compares the PENSIM results with 
those of the Department of Insurance under the escalating 0.A.S.- 


G.1I.S. assumptions for the moderate inflation proposed plan. 


PROJECTED INVESTMENT FUND LEVELS 
($ Billions) 


YEAR PENS IM FEDERAL? DIFFERENCE 


1975 

1980 1 
1985 1 
1990 2 
1995 2 
2000 1 


5. See page 52 of the Actuarial Report. 
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APPENDIX 


The ten source programs used in the PENSIM model are 
listed in this appendix for those who wish to duplicate the 
results of the Ontario model. All programs are written in 
standard FORTRAN - IV and were designed to operate on an inter- 


active basis using a PDP-10 computer. 


Other users are cautioned that this computer is an 
octal machine, differing from the IBM 360 or 370 series. Thus 
programming changes will be necessary to accomodate the PENSIM 


model on other computers. 


While the program versions in this appendix are the 
most recent ones, continuous updates and revisions are taking 
place. However, the programs listed here were used to generate 
the output in chapter xiv and any output can be compared on that 


basis. 
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VERSION AS OF 12:00 AUGUST 11,1972 
PROGRAMME POP. 1 
LAGRANGEAN POLYNOMIAL INTERPOLATION PROGRAM USING 
FOUR INTERPOLATION POINTS HOVING THROUGH TIME OM A FIFTER! 
YEAR INTERVAL. THE PROGRAM IS USED TO INTERPOLATE POPULATION 
STATISTICS, BETWEEN FIVE YEAR INTERVALS. 
PROGRAIi WRITTEN BY T.S.LETT ON JAN. 5 1972 FOR TI'E 
TAXATION AND FISCAL POLICY SRAMCH OF THE ONTARIO GOV'T. 


DIMENSION W(4,60),KSUB(8),FAC(4), XNUM(15),NO(24), POPII(1S, 13) 
+, POPF (19,13), POP(15,60, 2) 

DATA NO. /5,10,15,1,10,15,1,5,15,1,5,10, 20, 25, 30, 15, 25, 30,15, 
¥*20, 30,15, 20, 25/ 

DATA. KSUB 17577; 10, 13,.16,13, 227 

DATA W,FAC, XNUM, POP /240*0.,4*0.,15%*0.,2280*0./ 

DATA LM,NN,NNN,L,LL,LK /0,1,15,1,1,0/ 

CALL IFILE(20,'PRuUM1') 
CALL IFILE(21,'PRUJF1') 
CALL IFILE(22, 'PRUM') 
CALL IFILE(23,'PRUJF') 
READ(20,100)(POPM(1,1) 
READ(21,100)(POPF(I,1) 
READ(22,100)((POPM( I,J 
READ(23,100)((POPF(I,J 
FORMAT(12G) 

GO TO & 

L=2 

LK=15 

DO 10 J=1,4 

H=ELiN+l 

DO 10 i=1,15 
KK=KSUB(LM) 

IFCJ.NE.1) GO TO l 
fPELSEG. 2) GOs TO 1 

Ahi Sd a i ae AS 

GO TO 

FAC(J)=C(CLM-L) *5<NOCKK) ) #C CLMAL) #5-NOCKK4+1) )*CCLI-L) #5 
+-NO(KK+2)) 

YNUMCT )=CL4LK-NOCKK) )*® CLF LK-NOCKK+1) )*# CL FL K-NOCKKE2 )) 
WCJd, 1) =XNUMCL)/FACCJ) 

CONTINUE 

IFCLM.EQ.8) GO TO 5 

DO 11 I=NN,NNN 

LFCPCEO. 16.0R. 1.EQ.31.0R.1.EQ.46) LL=LL+3 

DO 11 J=1,19 

POP(J,1,1)=W(1,1)*POPM(JU,LL)+W(2, 1 )*POPHCJ,LL+1)4+(3, 1) 
+*POPM(JU,LL+2)+W(4, 1) *POPM(J, LL+3) 
POP(J,1,2)=WC1,1)*POPF(J,LL)+#W(2,1)*POPF(U,LL+1) +h (3, 1) 
+*POPF (J,LL+2)+W(4, 1 )*POPF(JU,LL+3) 
CONTINUE 

IFCNNN.EQ.60) GO TO 14 

GO TO 3 

DO 12 M=1,4 

DO 12.N81,15 

WOM, N+#15)=W(h, ND 

WOM, N+30)=W(M,N) 

W(M,N+45 ) =W(N,N) 

NN=16 

NNN=60 

GO TO 13 

CALL OFILE(24, 'ONE') 
WRITE(24,15)¢((POPC(I,J,K),J=1,60), 1=1, 
FORNAT(1H ,10(F6.1,1X)) 

STOP 

END 
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C VERSION AS OF 12:00 AUGUST 11,1972 
C PROGRAMME IMDEX. 2 
REAL YMPE(60), PNIND(60),CPI(60),ERIND(60),AVINC(G0),CPIHC(CO) 
REAL AVERN(60) 
DATA AVINC/1.0,1.07,1.066,1.065, 56*0, / 
DATA YMPE/5000.,5000.,5100.,5200.,5300.,5400.,5500.,53*0,/ 
DATA PNIND/109.4,109.4,111.6,113.8,116.1,118.4, 540, / 
DATA CPI/111.4,115.4,120.1,125.5,129.7,133.4, 54*0./ 
DO 100 J=1,3 
K=1972+J 
TYPE 110,K 
110 FORMAT(' SPECIFY YMPE IN ',14,': ',$) 
ACCEPT 120, YMPE(7#J) 
120 FORMAT(G) 
100 CONTINUE 
DO 130 J=1,6 
K=1969+4J 
TYPE 140, K 
140 ees SPECIFY GROWTH IN AVERAGE EARNINGS IN ',14,'CEG 1.035) 
as rp 
ACCEPT 120, AVINC(4tdJ) 
130 CONTINUE 
DO 210 J=1,6 
K=1969+J 
TYPE 220, K 
220 FORMAT(' SPECIFY CHANGE IN CPI IN ', 14," (EG 1.01): ',35) 
ACCEPT 120, CPINC(4+J) 
210 CONTINUE 
TYPE 142 
142 FORMAT(' SPECIFY CEILING ON PENSION INDEX (EG 1.02): ',$) 
ACCEPT 120, CEILNG 
DO 150 J=11,60 
150 AVINC(J)=AVINC(10) 
AVERN(1)=1. 
DO 160 J=2,60 
160 AVERN(J)=AVERN(J-1) *AVINC(J) 
DO 170 J=1,8 
170 ERBAS=ERBAS+AVERN(J) 
DO 180 J=11,60 
SUM=0. 
DO 190 K=2,9 
190 SUM=AVERN(J-K)+SUM 
180 ERIND(JU)=SUN/ERBAS 
DO 200 J=11,60 
200 YMPECU)=CIFIX(CYMPE(10)*ERIND(J))/100. ))*100. 
DO 230 JHll,60 
230 CPINC(JU)=CPINC(10) 
DO 240 J=7,60 
240 CPI(JU)=CPI (J-1)*CPINC(J) 
DO 250 J=2,60 
250 CPI(J)=((CPI(U)+CPI (J-1))/2. ) 
C THIS CALCULATION OF CPI GIVES US THE VALUE OF AVERAGE CPI For 
C THE 12 MONTHS ENDING ON JUNE 30 OF YR "yg" 
DO 260 J=7,60 
PNIND(J)=PNIND(J-1)*CEILNG 
IFCCPI(J-1).LT. PNIND(J))PNIND(J)=CPI (J-1) 
IFCPNIND(J). LT. (1. 01*PNIND(JU-1))) PINE (CU) =PMNIND(J-1) 
260 CONTINUE 
DO 270 J=1,60 
AVERN(J)=AVERN(J)/1.07 
270 PNIND(J)=PNIND(J)/1. 094 
CALL OFILE(25, 'TXFPB') 
CALL OFILE(26, 'TFPNX') 
CALL OFILE(27,'TFIND') 


WRITE(25,300) CYMPE(J),J=1,60) 
WRITEC(26,301) (PNIND(J),J=1,60) 
WRITE(27,302) (CAVERN(J),J=1, 60) 


301 FORMAT(10(F6.2,1X) 
302 FORMAT(10(F6.3,1X) 
STOP 
END 


( 
( 
300 FORMAT(10(F7.0,1X)) 
) 
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VERSION As OF 12:00 AUGUST 11,1972 

PROGRAMME MHMAE.3. THIS IS THE THIRD PROGRAMME IN TRE CPP 
SET OF PROGRAMMES.IT IS USED TO CALCULATE MODIFIED MO0D- 
IFIED AVERAGE EARNINGS USING THE METHODS OF THRE DEPT. OF 
INSURANCE IN OTTAWA.FOR ADDITIONAL INFORMATION CONSULT 
THE CPP DOCUMENTATION FILE.PROGRAMME WRITTEN BY G.JACOSS 
AND T.LETT FOR THE TAXATION AND FISCAL POLICY BRANCH, 
JAN-FEB 1972. 


3.09 © C2 Cl Co C3: CoG 


REAL YMPE(60),AVERN(11,60, 2), INDEX(60),YPE(60),ADU(22),FILE(10) 
COMMON /IODATA/ VALU(11, 60,8) 
DATA FILE/'RSCUR', 'MICUR", "RSPRO', 'MIPRO', 'TXFPB', "MMRSC', "MMIC 
+" "MMRSP!', 'MMMIP!, "MMTFP!/ 


C 

C THIS SECTION IS USED TO INPUT VALUES OF YMPE THAT ARE TO BE 
C USED IN THE PROGRAMME CPP. THE FOUR POSSIBILITIES ARE: 

C RSCUR WHICH IS REASONABLE STABILITY UNDER THE CURRENT PLAN 
C RSPRO WHICH IS REASONABLE STABILITY UNDER THE PROPOSED PLAN 
C MICUR WHICH IS MODERATE INFLATION UNDER THRE CURRENT PLAN 

C MIPRO WHICH IS MODERATE INFLATION UNDER THE PROPOSED PLAIN 

C THE VALUES ARE STORED IN THE DATA FILE WITH TRE SAME NANE 

C WITH TEN VALUES PER LINE. THE VALUES OF YMPE ARE ODTAINED CY 
C USING THE ASSUMPTIONS ON PG. 23 OF THE ACTUARIAL REPORT FOR 
C THE INFLATION FACTORS AND BY APPLYING THE EARNINGS INDEX FOR 
C THE YEARS 1966-2025 

C 

& TYPE 20 

1 FORMAT(10G) 


20 FORMAT(' TYPE IN PLAN TO BE TESTED: ',$) 
ACCEPT25, FNAME 

25 FORMAT(AS) 
DO 26 IK=1,5 

26 |FCFNAME.EQ.FILEC(IK))GO TO 28 
TYPE27, FNAME 

27 FORMAT(1X,A5,' ?, TRY AGAIN'/) 


GO TO 2 
28 CALLAOETEECZS,ETLECLK+5)) 
CALL IFILE(22, FNAME) 
READ(22,1) “(YMPEC 7, 1=1,60) 
CALL IFILE(28,"ADJFC') 
READ(28,427) (ADUJ(1),1=1, 22) 
427 FORMAT(9G) 


THE NEXT SECTION INPUTS THE AVERACE EARNINGS OF MALES AND FEMALE 
IN 1967 AND MULTIPLIES THIS FIGURE BY AN INDEX TO REFLECT 

THE INCREASE ANNUALLY AND THUS CALCULATE AVERAGE EARNINGS BY AGE 
AND SEX IN ALL YEARS 1966-2025. THE INDEX FACTORS 

ARE OBTAINED FROM THE ASSUMPTIONS OF THE ANNUAL INCREASE 

IN AVERAGE EARNINGS FOUND ON PAGE 23 OF THE ACTUARIAL 

REPORT. THE REASONABLE STABILITY INDEX IS STORED IN RSIND.DAT 

AND THE MODERATE INFLATION INDEX IS STORED IN MIIND.DAT 


hehe vel se hele be ie Bs, 


CALL IFILE(23, 'AVE67') 
READ(23,1) (CAVERN(I,2,K), 1=1,11),K=1, 2) 
GO TO (30,32, 30, 32,31),I1K 
32 CALL IFILE(24, 'MIIND') 
GO TO 35 
30 CALL IFILE(24, 'RSIND') 
GO TO 35 
31 CALL IFILE(24, 'TFIND') 
35 READ(24,1) (CINDEX(1! ),1=1,60) 
DO 4O K=1,2 
DO 40 J=1,60 
DO 40 1=1,11 
4O AVERNC(I,J,K)=AVERN( 1, 2,K)*INDEX(J) 


C THE NEXT SECTION IS USED TO CALCULATE YBE WHICH IS 12% OF Y?'PE 
C ROUNDED DOWN TO THE NEAREST 100 IN THE CURRENT PLANS AND A 

C FLAT AMOUNT OF 600 IN THE PROPOSED PLANKS 
C 


GO TO (60,60,61,61,62),I1K 


! 


ee tay 
ns = La “a a0 aay 


" fi ivy 
a On, va au 
Pay «Yo 9h INR uy ‘ng 


Pave? Sik ane | 

| | | . 7 

fe) ata Aa era ee, coal a ) 
A eee Bad ee WAL SUT sae te Sh a Ba ain | 

. hy} Ly ’ j 1 tigi . iy Nite an Eph ay rea mie i bs hp is ! iat vs ~ 

as es tu th) VR ‘b PA Yt ea: i iu ore iin %, 

Peein “Ys 8 AO TAD ail 

bg Seay we eS our ROP sabe ; 

tay | pet ae AVAL war ae tea | z 

ne 7. ay’ wage H SAY a. 1b vane a4 fh # ig a “a7: 4 


(id 7 Ast ; 
"s 7 


JM rt ah oe Abr aa MoT 
A . i . : ‘Sc wY7 2 
. Pe CVT 


_) 


t oer a “Beh 


o 
pes 
= 


[ay 
Leavy Ane a 
t fen) a i 
of. Oana + Sh Sua 7 
Mudie . ESV : 


a Vat ue. et age es 


xAVe 9? i 4) PH a ge 


MAS LSPS eee. 
. VSS OP GSeSRae 


. 
ay) 
q 


&) 42 


4 
4 
Pee 
ea 


yp at TrAGaAg. oT ares ia or 
Ny eds CO? Ree Ga Gir “aA 
hs ‘a cv OMS A): Gd hav Te 4S 
2B LB FOB a Aa 
wy 
ve ; : . 


» & 


ap Seecs ta & 


62 TYPE63 
G3 FORMATC1X'YBE CALCULATION CURRENT/PROPOSED ?'/) 
ACCEPT25,ANS 
IFCANS. EQ. 'CURRE') GO TO 60 
61 DO 55 1=1,60 
55 YBE(1)=600 
GO TO 70 
60 DO GS 1=1,60 
6S YBECI)=(IFIX(.12*YMPEC1 )+1)/100)*100 
Cc 
C WE INSERT THE +1 ABOVE TO ACCOUNT FOR INTEGER DIVISION 
C THE NEXT SECTION CALCULATES THE X VALUES REQUIRED FOR THE SPLINE 
C SUBROUTINE TO CALCULATE THE CORRESPONDING Y VALUES USED IN 
C THE CALCULATION OF MODIFIED NODIFIED AVERAGE EARNINGS. 
C 
70 JJ=4 
DO 81 K=1,2 
DO 80 J=1,60 
DO 80 121,11 
VALU(I,JU,K+dd )=100*#YMPE(J)/AVERN(I,J,K) 
80 VALU(I,J,K+1+JJ)=100*YBE(U)/AVERN(1,J,K) 


81 JJ=5 

CALL SPLINE 
C 
C THE FOLLOWING SECTION CALCULATES MODIFIED MODIFIED AVERACE 
C EARNINGS IN ONE STEP.THE INTERPOLATED Y VALUES FROM SPLINE 
C ARE USED TO CALCULATE THE TEMPORARY MODIFIED AVERAGE EARN- 
C INGS WHICH ARE IN TURN USED TO CALCULATE THE MODIFIED AVE- 
C RAGE EARNINGS FROM WHICH MODIFIED MODIFIED AVERAGE EARNINGS 
C ARE DERIVED.SEE COLUMN 10 WORKSHEET DO-5 FOR METHOD. 
C 

JJ=0 

DO 90 K=1,2 


po 91 I=#1,11 
DO 91 J=1,60 
91 VALUC(I, Ju, K)=(CCAVERNC 1, Ju, K) *(VALUCI, J, Jd+3)-VALUCI,J, 
&Juth)))+CYMPE(J) *(100-VALU(1,JU,Jut1))))/C100-VALUCI,d, 
&IU+2))) ¥ADI C1 + (Ud*3 -K41) #CYMPEC 1) /YMPECJ)) 
$0 JJ=4 


THIS CALCULATION IS USED TO OBTAIN TRE MODIFIED 

MODIFIED AVERAGE EARNINGS LEVELS BY AGE GROUP AND SEX AND 
YEAR (1966-2025). THE METHOD USED IS ADAPTED FROM TRAT OF 
THE FEDERAL DEPT. OF INSURANCE THAT THEY USED IN THE 1969 
CPP ACTUARIAL REPORT. THE MATRIX AVERN IS EXPLAINED ABOVE. 
THE RESULTANT 'Y' FACTORS FROM SPLINE ARE STORED IN VALU, 
THE VECTOR ADJ CONTAINS TWO GROUPS OF 11 ADJUSTMENT FACTORS 
(MALE AND FEMALE) WHICH SCALE THE RESULTS FROM SPLINE. 

THIS OPERATION IS NECESSARY BECAUSE THE DATA USED IN 

SPLINE RELIED ON THE 40-44 AGE GROUP. THE ADJUSTMENT FACTCR 
SCALES THE RESULTS INTO 11 AGE GROUPS INSTEAD OF ONE. 

THE VALUES IN THIS VECTOR WERE OBTAINED FROM WORKSHEETS 
DO-S AND DO-8 FROM THE DEPT. OF INSURANCE. THE FINAL FACTOR 
WHICH IS THE RATIO OF YMPE'S IS USED TO BRING THE MODIFIED 
AVERAGE EARNING FIGURE BACK TO 1966 CONSTANT VALUES AND 
THUS OBTAIN THE MODIFIED MODIFIED AVERAGE EARNINGS 

MATRIX FOR ALL FURTHER CALCULATIONS IN TRIS MODEL. 

TYPE 94,FILECIK+5) 

94 FORMATC' OUTPUT FILE IS AS, RAT Z2 


OOO t3..) Oo) Ot Ooo oo C2 ot) 2D 


WRITE(29,95) (((VALUC I,J, K),J=1,60), 1=1,11),K=1,2) 
95 FORMAT(5(F7.0,1X)) 


THE OUTPUT IS STORED IN FILES WITH AN "HM" PREFIX(MODIFIED 
MODIFIED) WITH A CORRESPONDING SUFFIX FROM THE OPTION 
TESTED. THESE FILES ARE STORED ON DECTAPE D&23 FOR FUTURE 
USE IN PROGRAMME RETBH.4. 


RPOAOOMNMNNINNIO 


00 FORNAT(' ',9(F7.3,1X)) 
STOP 
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END 
SUBROUTINE SPLINE 


THIS SUBROUTINE INTERPOLATES Y VALUES FROM GIVEN 

X VALUES USING THE CUBIC SPLINE METHOD.A SET OF ORDERED PAIRS 
CREATES AN INTERPOLATING EQUATION. THE ORDERED PAIRS ARE READ 

IN FROM THE FILE POINT.DAT ALONG WITH THE FIRST DERIVATIVES CF 
THE END POINTS(D(1),D(2)).THE INTERPOLATING X VALUES ARE CARRIED 
THROUGH COMMON IN THE ARRAY VALU. THE INTERPOLATED Y VALUES 

ARE RETURNED THROUGH THE SAME ARRAY TO THE MAIN PROGRAMME, 


oe AaNIAaAANAANAYN 


REAL X(31),Y(31),C0(100),W(100) 
INTEGER KSUB(8),D(2) 
COMMON /IODATA/ VALU(11, 60,8) 
DATA C,W,K1,N,KSUB/200*0.,0,31,1,2,3,4,5,6,7, 8/ 
CALL IFILE(20,'POINT') 
IF(K2.EQ.8) GO TO 1 
IFCK1.EQ.0.0R.K1.EQ.2) MI 
IF(K1.EQ.4.0R.K1.EQ.6) MI 
READ(20,111)uU, (D(1) 
READ(20,111)(X(1),1 
READ(20,111)(Y(1),1 
11 FORMAT (10G/9G/9G/3G 


1 

Cc 

C THIS SCALING IS USED ONLY FOR THE CPP DATA 

C AND SHOULD BE REMOVED FOR ANY OTHER DATA SOURCES. 
C 


D 
X 
35 Ms 
Q 
Y #°¥(02) 


100 C(2)=0.0 
W(2)=0.5*Q*Q*D(1) 
200 M=N=-2 
IF(M.LE.0) GO TO 350 
DO 1000 !=1,M 
Al=Q 
Q=X(1+2)- X(1+1) 
H=A1/Q 
C(3*1 )=-H/(2.0-C(3*!I-1)) 
WC3¥1 )=(-YI-WC3*1-1))/(2.0 - C(3*I-1)) 
C(3%1+1)=-H*H/(H-C(31 )) 
WO3*1+1)=CYI-WC3*1))/(H-CC3*1 )) 
V¥heYClez)~ YCR+1) 
C(3*1+2)=1.0/(1. 0-C(3*1+1)) 
1000 WC 3*1+2)=(YI-W(3*14+1))/(1.0-C(3*1+1)) 
350 IF(J.EQ.1) GO TO 400 
C(3*N=-3)=(Q4Q*D(2)/2. O-W(3*N-4))/(3.0- CC3*N-4)) 
GO TO 500 
400 C(3*N-3)=(Q*D(2 )-YI-WC3*N-4) )/(02. 0-CC3*K-4)) 
500 M=3*N-6 
[F(i4.LE.0) GO TO 700 
DO 600 I1=1,M 
=[.i- | [+3 
600 CC 1 )=WCL)-CC1)*CC1+1) 
700 IF(J.EQ.1) GO T9 800 
C(1}=¥(2) = ¥C1)-WU2)-C(3) 
2) 


t 

800 C(2) 
8 K1=K1+1 
K2=KSUB(K1) 
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LATEST UPDATE 12:00 AUGUST 11,1972 

PROGRAM RETBH.4 IS USED TO CALCULATE THE COSTS OF THE 
CANADA PENSION PLAN RETIREMENT BENEFITS. THE PROGRA!S 
REQUIRES A DATA FILE OF MODIFIEC MODIFIED AVERAGE 

EARNINGS (SEE PROGRAM MMAE.3) AND A DESICNATION OF 

BEING TESTED (INPUT FROM TTY). 
MMAE(11,60,2), PART(11,10), PMMAE(11,60, 2),SUM(-12/11, 2) 


FILE(25),P0P(19,60, 2), TEM1(-12/11),FACTOR(15, 60, 2) 
TPNER(-12/11, 2), TEM(-12/11, 2),DSDRP(-12/11,9) 
YMPE(60), PNIND(60),ESC1(60),RTFAC(-12/11,2),DREDC(60) 
PINDIV(-12/11,47, 2), RTCST(60, 3), ESC2(60),RFAC(2,5,2) 


SURVRT(15, 8, 2), NMDTH(15, 60, 2),DTHCS1(60, 3),DTHCS2(60, 2) 


INTEGER JJ(22),NN(13), YR(60),FRSTYR(-12/10), LASTYR(10) 
INTEGER YRSPRT(-12/10) 
COMION/ I ODATA/INDIV(-12/11,47, 2), EXPOS(-12/11),DROP(-12/11,2), 


se? 


DATA 


Re Re L wa 


DATA 


NR,DRPRT 

PAPT/.65, 3*1. 0, 3%.93, 2*. 87,.78, .53,.68, 3*1.0, 3%. 94, 2*. 88, 
679,.54, .46, .66,.44,6*%. 36,.25,.13,.5,.69,.5,6%.42,.29,.15,.95, 
.73,.57,6%.48,.34,.14,.58,.75,.58,6%.49,.35,. 14, .46,.66,.44, 
6#,36,.25,.23,.49,.68,.5,6%.81,.29, 613, 253,671, .56,/6%.47, 
3,5 1 OSS 7 ore TZ6 ON ee els / 

JJ/5, 10,15, 20, 25, 30, 35,40,45,50,55,60,65, 70, 


&75, 80, 85,90,95,100,105,110/ 


DATA 
DATA 
DATA 
DATA 


&! 


NN/8,13,18, 23, 28, 33, 38,43,48,53,58,63,68/ 
FRSTYR/14%18, 21, 26, 31, 36,41,46,51,56,61/ 
LASTYR/20, 25, 30, 35,40,45,50,55,60,65/ 

FILE/ 

RSCUR', "MICUR', "RSPRO', 'MIPRO', "TXFPB', 'MMRSC', THliMIC', 'HURSP', 


& 'MMMIEP', "MMTEP', 'DHRSC', 'DHMIC', 'DHERSP', 'DHMIP', 'DETFP', 'RTRSC', 


oO. 


& 
& 


CALL 


'RIMIC', (RIRSP', 'RTMIP', "RTTFP", 'DBRSC', 'DEMIC', 'DPRSP®, “DENI P 
'DBTRP*/ 


DATA DSDRP/14*1. 002,101. , 1441. 003,1.002,9*1.,14*1. 004, 2*1.003,5 
&. 


,14%*1.005, 2*1.004,1.003, 7*1.,14#1.006,2*1.0055,1.0045,1. 005, b¥1 
16*1.007, 2*1.006,1.004,5*1.,18*1.0085,1.0075,1.005,4*1.,16*1.0 
3*#1,011,1.01,1.007, 3*1.,17*1.015,1.016,2*1.017, 3*1.012,1./ 

RFAC/61.,71.,86.,11%.,115.,167.,148.,232.,189.,509.,; 

PA CNT2 6 S84 205 PELE Ss 50s 8507 19 Ole F 

CLOCK1 


C THE NEXT SECTION IS USED TO INPUT THE VARIOUS POLICY OPTIONS 


Cr TRA 


ARE TO BE TESTED FOR THIS RUN 


19 -TYPE.20 
20 FORMAT(' TYPE IN PLAN TO BE TESTED: ',$) 
ACCEPT 25, FNAME 


DO 21 


IK=1,5 


21 IFCFNAHE.EQ.FILECIK))GO TO 24 


Titre 


23, FNAME 


23 FORMAT(1X,A5,'? TRY AGAIN! '/) 
GO TO 19 

22 FORMAT(1G) 

25 FORMAT(A5) 

24 TYPE 26 

26 FORMAT(' TYPE IN DROPOUT RATE (0-1): ',$) 
ACCEPT22, DRPRT 


Y 2 
07 HORMATC! TYPE IN RETIREMENT PENSION RATE (0-1): ',9) 
ACCEPT22, RETRT 


TYPE 


28 


28 FORMAT(' DO YOU WANT EARNINGS TEST? (YES/NO): °,§$) 
ACCEPT29,ERNTST 
29 FORMAT(A2) 


TYPE 


30 


30 FORMAT(' OUTPUT? (TADLE/DATA): +62 
ACCEPT29, OUT 


CF 4a 62 


THE NEXT SECTION READS IN THE YMPE VALUES, THE POPULATION DATA, 
THE SURVIVAL RATES, AND THE PENSION INDEX TO BE USED DEPENDING 
ON THE POLICY VARIABLES BEING TESTED, 


CALL [FILE (22, FNAME) 
READ (22,1) (YMPECI),1=1,60) 
1 FORMAT(10C) 


CALI 


IEILE (2k. 'ONE') 


P wecgueead) 


i. i 
p . 
v6 “ ; Wel n ; \ i 
r i, 7) 4 f 
‘ 1 I My 7 
® és i ) 
i@ 
i 
i i 
Me 
X 
| 
f 
= 
+4 


io >% 
“7 bi 
Z 


ROM ane 
ie eka ae 


‘it & a} 


thes 


i OB, ' 
t 
Mi 0 


hy re ay , 


De 
ai % eat Ay 

ne Hf y wee en, ry p ly (e 
Nhe ak Ws 

Sais! Ned ae) fy 
et Bis tal) ie # 

fy een. Ate at ry 
Vane ory 

i a wn ie by au G ays 

Se we 

Ak 


Wer 


? 
re" 
j 
of 


Ne 


ea / 
car 


cok = re 


Ki ee 

“ a ra ° 
ose a een i 
aaa haa ria! % = 7 > 


’ 
F-7t tie BE ke. | 5p) Fi _ 
A a ; Bed 
Le Ter . AV ea ay 4 
; : i if : oh? a i, 
v » 200% ih lin a ' 7 
4 Y " in Nn ; | F an ; a 


ie £49 nal 


e « 
. i We 
me) 


7.) 
ae) bo Me 

tai 
4 i 


; 73h 
ay ate arg 

‘an... ) ivi 
Pe > a 

hei nJaRR 
a amen” | 
| er era. ay 
wie. enor 


isin ria? 7b 


Pag | os 
jt gat 


\ + if Ae; 
Lara vi¥) 
ne », vi 


La ? ¢ i niger. : 


ix 


READ(24,1) (((POP(1,J,«),JU=1,60),1=1,19),K=1,2) 
CALL IFILE(25, 'SURVI') 
READ(25,9) (((SURVRT(1,J,K),J=1, 8), 1=1,15),K=1,2) 
9 FORMAT(8G) 
GO, TO: €36, 37, 36,737,358) 1K 
36 CALL IFILE(23,'RSPNX') 
CALL IFILE(26,'RSDRD') 
GO TO 39 
37 CALL IFILE(23, 'MIPNX') 
CALL IFILE(26,'MIDRD") 


GO 70 39 
38 CALL IFILE(23, 'TFPNX') 
Tree. < 


2 FORMAT(C' ECONOMIC ASSUMPTIONS? (RS/MI): ',$) 

ACCEPT 29,ANS 

IFCANS.EQ.'MI')GO TO 37 

GO TO 36 
39 READ(23,1) (PNIND(I),1=1,60) 

READ(26,1) (DREDC(1),1=1,60) 
C THE HEXT SECTION INPUTS THE APPROPRIATE DATA FOR THE MODIFIED 
C MODIFIED AVERAGE EARNINGS AS CALCULATED IN PROGRAM MMAE.3 
3 CALL IFILE(21,FILEC(IK+5)) 

CALL OFILE(27,FILEC(IK+20)) 

CALL OFILE(28,FILE(IK+15)) 

CALL OFILE(29,FILEC(IK+10)) 

READ(21,35) ((C(MMAECI,J,K),J=1,60),1=1,11),K=1,2) 
35 FORMAT (5G) 

1ECCIK/2)*2.EQ. IK.OR. (IK. EQ. 5.AND.ANS.EQ. 'MI'))GO TO 5 

IMl=4 

GO TO 6 
5 IM1=0 
6 JJJ=0 

DO 8 K 

DO 7 | 

pO 7 J 
IF C(K.EQ.1.AND.JU.GE.20) N1=2 
IF (K.EQ.2.AND.U.GE.1) N1=3 
IF (CK. EQ.2.AND.U.GE.10) N1=4 
IFCK.EQ.2.AND.JU.GE.20) N1=5 
IFCK.EQ.2.AND.JU.GE.40) N1=6 


C THE NEXT CALCULATION MULTIPLIES THE MATRIX OF MODIFIED MODIFIED 
C AVERAGE EARNING BY THE APPROPRIATE LABOUR FORCE PARTICIPATION 
-C RATE DEPENDING ON THE SEX,THE YEAR AND THE PLAN BEING TESTED 
7 PMMAE(I,JU,K)=HMMAECI,J,K)*PARTCI,JJJ+N1) 
8 JdJ=JddtIM! 
C THE NEXT SECTION PERFORMS A MATRIX TRANSFORMATION WHEREBY THE 
C MMAE'S FOR ANY AGE GROUP (E.G. AGE -25 IN 1965; AGE 0 IN 
C 1965; ACE 35 IN 1965) ARE TRANSFORMED INTO A ROW OF 
C VALUES. THIS SIMPLIFIES THE DIAGONAL STEP-WISE ADDITION OF 
C MMAE'S FOR ALL AGE GROUPS AND ALSO SIMPLIFIES THE CALCULATION 
C OF EARNINGS DROPPED OUT IN SUBROUTINE MIN. THE TRANSFORMATION 
C 1S DONE FOR THE MATRIX OF MMAE'S AND THE MATRIX OF PARTICIPATI?C 
C MMAE'S (PMMAE). THE TRANSFORMED MATRICES ARE CALLED INDIV AND 
C PINDIV 

DO 70 K=1,2 

DO 70 1=2,10 

DO 70 J=1,4 


INDIV(I,J,K)=MMAE(1,u,K) 
70 PINDIV(I,J,K)=PMMAE(1,u,K) 
DO 80 K=1,2 


DO 80 1=2,9 
DO 80 L=Ii+1,10 
MM=JJCL-1) 
DO 80 Jai, MM+4 
INDIV(IL,J,K)=iMAE(L,J,K) 
80 PINDIV(I,Ju,K)=PHMAE(L,J,K) 


DO 90 K=1,2 
DO 90 1=0,-12,-1 
N2=NN(-(1-1)) 
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1F(N2.GT.60)N2=60 
INDIVCI,1,K)=MMAE(1,N2,K) 
INDIVCI,2,K)=MMAE(1,N2+1,K) 
PINDIV(I,1,K)=PMMAE(1, N2,K) 

90 PINDIV(I,2,K)=PMMAE(1,N2+1,K) 
J=0 
DO 100 K=1,2 
DO 100 I=1,-12,-1 
DO 100 L=2,10 
MM=Ju(L=1) 

DO 100 J=liM, MM+h 
NA=J-((-1*5)+7) 

IJ=uJ 

IF(Id.GT.60)1U=60 
INDIVCL, NA, K)=MMAECL, Iu, K) 

100 PINDIVCI,NA,K)=PMMAE(L, Iu, K) 
DO 105 K=1,2 
INDIV(1,1,K)=MAE(1, 3,K) 
INDIV(1,2,K)=MMAE(1, 4,K) 
PINDIV(1, 1,K)=PMMAE(1, 3,K) 

105 PINDIV(1, 2,K)=PMMAE(1,4,K) 
DO 192 NR#=1,9 
DO 108 K=1,2 
DO 108 I=NR-10,NR+1 

108 SUM(I,K)=0. 

DO 110 |=NR-10,NR+1 
YRSPRT(1)=LASTYR(NR+1)-FRSTYR(1) 
IF CYRSPRT(1).LT.0)YRSPRT(1)=0 

110 EXPOS (1)=DRPRT*YRSPRT(1) 

DO 120 K=1,2 
DO 120 1=NR-10,NR+1 
DO 120 J=1,YRSPRT(1) 

120 SUM(1,K)=SUM(I1,K)+PINDIV(1,u,K) 
CALL MIN 
DO 140 I1=NR-10,NR+1 


Pi=l 
IFCIi1.LE.1)DROP(11,1)=.5*DROP(I1,1)*.65 
[FCIIL.LE.1)DROP(II,2)=.5*DROP(I1,2)*.5 
1FCI1.GT.1)DROP(I1,1)=,.5*DROP(I1,1)*PARTCNR+1, 1) 
IFC1E.GT.1)DROP(11,2)=.5*DROP(I1,2)*PART(NR+1, 3) 
140 CONTINUE 
C CAUTION*****DROP 1S BEING OVERWRITTEN 
C THIS VALUE OF .5 IN THE ABOVE CALCULATION IS TO TAKE ACCOUNT 
© OF THE NIL EARNINGS OF CONTRIBUTORS AND THUS ONLY USE 
C HALF THE MAXIMUM DROPOUT AMOUNT. 
DO 160 K=1,2 
DO 160 !=NR-10,NR+1 
160 TPNER(I1,K)=SUM(1,K)-DROPC(I,K) 
DO 170 1=NR-10,NR+1 
TEM1(1)=YRSPRT(1)*(1.-DRPRT)*(1./RETRT) 
PFECTEM1C1).LT.40. )TEM1C1)=YRSPRTC1)*(1./RETRT) 
170 CONTINUE 
DO 180 K=1,2 
DO 180 I=NR-10,NR+1 
180 RTFAC(I,K)=(TPNER(I,K)/TEM1(1))*DSDRPC1,NR) 
DO 192 K=1,2 
DO 192 1=NR-10,NR+1 
JeJU(NR-1+2) 
192 FACTOR(NR+1,U,K)=RTFAC(I,K) 
DO 195 K=1,2 
DO 195 J=5,60,5 
195 FACTOR(1,U,K)=FACTOR(2,J,K)-(( (FACTOR (2,U,K)+FACTOR(3,U,K))/2.) 
& -FACTOR(2,uU,K)) 
DO 200 K=1,2 
DO 200 N=2,9 
DO 200 J=5,60,5 
= + N+ 
200 FACTOR CH, y/K) ((FACTOR(N, J, K)+FACTOR(N+1,J5,K))/2.) 


DO 203 I=-1,10 
90% TEMCI.KY=TPNFROI.KY/TEMICI) 
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DO 205 K=1,2 
nO 204 I=-1,9 
20h RTFAC(I,K)=((TEMCI,K)+TENC1+1,K)) 
905 RTFAC(10,K)=1EM(10,K)*DSDRP(10, 9) 
DO 208 K=1,2 
DO 208 !=-1,10 
J=Jd(11-1) 
208 FACTOR(10,JU,K)=RTFAC(I,K) 
pO 210 K=1,2 
DO 210 J=5,55,5 
210 FACTOR (11,U+5,K)=FACTOR(10,u,K) 
FACTOR(11,5,1)#172. 
FACTOR(11,10,1)=643. 
FACTOR(11, 5, 2)=30. 
FACTOR(11,10,2)=142. 
FACTOR( 11,15, 1)=. 3273*FACTOR(11,10,1)+.6727*FACTOR(11, 20,1) 
FACTOR(11,15,2)=.2072*FACTOR(11,10,2)+. 7928*FACTOR( 11, 20, 2) 
THESE ADJUSTMENTS ARE NECESSARY TO ACCOUNT FOR THE SHORT 
RANGE PROJECTION OF DEATH BENEFIT FACTORS 
DO 220 K=1,2 
DO 220 N#1,11 
DO 220 J=5,55,5 
DO 220 M=1,4 
220 FACTOR(N, Uti, K)=M*. 2*(FACTOR(N, J+5,K)-FACTOR(N,U,K))+FACTOR(L, 
Racdok) 
DO 230 J=3,60 
230 ESC1(J)=( CC CYMPE(J )+YMPE(JU-1)+YMPE(d-2))/(3./RETRT))/ 
& (PNIND(J)/100.))/1250. ) 
DO 235 J=3,60 
235 ESC2(J)=C((CYMPE (J )+YMPE(J-1)+YMPE(U-2))/(3./RETRT) )/1250. ) 
DO 240 K=1,2 
DO 240 J=11,60 
240 FACTOR(11,uU,K)=FACTOR(11,uU,K)*ESC1(U-5) 
DO 250 K=1,2 
DO 250 N=12,15 
DO 250 J=5,55 
250 FACTOR(N, J+5,K)=FACTOR(N-1,U,K) 
DO 260 J=11,60 
FACTOR(10, JU-5,1)=FACTOR(11,J, 1) 
260 FACTOR(10,uU-5, 2)=FACTOR(11, Jd, 2) 
DO 270 J=56,60 
FACTOR(10,u,1)=FACTOR(10,JU,1)*ESC1(U) 
270 FACTOR(10,u, 2)=FACTOR(10,J,2)*ESC1(J) 
WRITE(29,41) ((CFACTOR(N,J,K), J=1,60),N=1,15),K=1,2) 
DO 310 K=1,2 
DO 310 1=1,15 
DO 280 J=5,9 
280 NMOTH(1,u,K)=(1.-SURVRT(1,2,K))*POP(1 +4, U,K)*FACTOR(1,u,K) 
DO 285 J=10,14 
285 NMDTH(1,Jd,K)=(1.-SURVRT(1,3,K))*POP(I+4,U,K)*FACTOR(1,J,K) 
DO 290 J=15,19 
290 NMDTHC 1, J, K)=(1.-SURVRT(1,4,K))*POP(1+4,d,K)*FACTOR(I,J,K) 
DO 295 J=20, 24 
295 NMDTHCL, JpK)=(1. -SURVRT(1,5,K))*POP(1+4,U,K)*FACTOR(I,d,K) 
DO 300 J=25,29 
300 NMDTH(1, JpK)=(1. -SURVRT(1,6,K))*POP(1 +4, U,K)*FACTOR(1,J,K) 
DO 305 J=30, 34 
305 NMDTH(1,Jd,K)=(1.-SURVRT(1, 7,K))#POP(I+4,d,K)*FACTORCI,J,K) 
DO 310 J=35,60 
310 NMDTH (I,J, K)=(1. -SURVRT(1,8,K))*POP(1+4,U,K)*FACTOR(I,J,K) 
DO 315 K=1,2 
DO 315 J=5,60 
po 315 1=1,9 
315 DTHCS1(U,K)=DTHCS1(J,K)+NMDTH(1,J,K) 
DO 320 K=1,2 
DO 320 J=1,60 
320 DTHCS1(J,K)=DTHCS1(U,K) *ESC2 CU) 
DO 325 K=1,2 
DO 325 J=5,60 
pO 325 1=10,15 
295 NTHCS?(1.K SN THCS2 (41. K+#NMDTHCI.J.K) 
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DO 330 K=1,2 
DO 330 J=1,60 
330 DTHCS2(JU,K)=DTHCS2(JU,K)*(PNIND(U)/100. ) 
DO 335 K=1,2 
DO 335 J=1,60 
335 DTHCS1(J,K)=(DTHCS1(J, K)+DTHCS2(J,K))*.5 
DO 341 J=1,60 
341 DTHCS1(J, 1) =#(1.-DREDC(J)) *DTHCS1(u, 1) 
DO 345 J=1,60 
YRCJ)=1965+d 
345 DTHCS1(J, 3)=DTHCS1(J,1)+DTRCSI1 Cu, 2) 
|F(OUT.EQ.'DA')GO TO 348 
WRITE(27,10) FNAME 
10 FORMAT (! DEATH BENEFIT COSTS UNDER PLAN ',AS, ($000''S)') 
WRITE(27,13) 
WRITE(27,14) 
DO 347 J=5,60 
347 WRITE(27,17) YR(J), (DTHCS1(J,K),K=1, 3) 
GO TO 999 
348 WRITE(27,999) (DTHCS1(J, 3),J=5,60) 
999 FORMAT(5(3PF14.0)) 
DO 1208 K=1,2 
DO 1208 |=-1,10 
J=JJ( 11-1) 
1208 FACTOR(10,J,K)=RTFAC(I,K) 
DO 1210 K=1,2 
DO 1210 J=5,55,5 
1210 FACTOR(11,J+5,K)=FACTOR(10,u,K) 
FACTOR(10,10,1)=241, 
FACTOR(10,10,2)=95. 
FACTOR(11,10,1)=410. 
FACTOR(11,10,2)=85. 
FACTOR(11,15,1)=.25*FACTOR(11,10,1)+. 75*FACTOR(11, 20,1) 
FACTOR(11,15,2)=. 35*FACTOR(11,10,2)+.65*FACTOR(11, 20,2) 
DO 1215 K=1,2 
DO 1215 N=10,11 
DO 1215 J=5,9 
1215 FACTOR(N,J,K)=RFAC(N-9,JU-4,K) 
C THESE ADJUSTMENTS ARE NECESSARY TO ACCOUNT FOR THE SHORT 
C RANGE PROJECTION OF RETIREMENT FACTORS 
DO 1220 K=1,2 
DO 1220 N=10,11 
DO 1220 J=10,55,5 
DO 1220 M=l1, i 
1220 esc J+, K)=M*, 2*(FACTOR(N, J+5,K)-FACTOR(N, J, K))+FACTOR(N, 
& J,K 
DO 1240 K=1,2 
DO 1240 J=11,60 
1240 FACTOR(11,U,K)=FACTOR(11,U,K)*ESC1(JU-5) 
DO 1250 K#1,2 
DO 1250 N=12,15 
DO 1250 J=5,55 
1250 FACTOR(N, J+5,K)=FACTOR(N=1,J,K) 
DO 1260 K=1,2 
DO 1260 J=11,60 
1260 FACTOR(10,JU,K)=FACTOR(10,JU,K)*ESC1(J) 
IFCERNTST.EQ."NO')GO TO 357 
DO 355 J=11,60 
FACTOR(10,J,1)=. 5*FACTOR(10,J,1) 
355 FACTOR(10,J,2)=. 75*FACTOR(10,Ju, 2) 
C THE .5 AND .75 IN THE PREVIOUS STATEMENTS ALLOW FOR THE EARNINGS 
C TEST BETWEEN THE AGES OF 65-69, AND BY OMITTING THIS 
C DO-LOOP, THE EARNINGS TEST IS NEGATED 
357 DO 365 K=1,2 
DO 360 N=10,15 
DO 360 J=1,60 
360 FACTOR(N,J,K)=FACTOR(N, J, K)*POP(N+4,J5,K)*(PNIND(J)/100. ) 
C CAUTI ON****HATRIX FACTOR 1S BEING PARTIALLY OVERWRITTEN 
DO 365 J=1,60 
DO 365 1=10,15 
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365 RTCST(JU,K)=RTCST(U,K)+FACTOR(1,J,K) 
DO 370 J=1,60 
370 RTCST(J, 3)=RTCST(J,1)4RTCSTC(J, 2) 
DO 375 K=1,3 
DO 375 J=1,10 
375 RTCST(J,K)=.98052*RTCST(uU,K) 
C THIS IS AN ADJUSTMENT FOR THE SHORT TERM FOUND ON SHEET F100-1/2 
1\FCOUT.EG. 'DA')GO TO 381 
WRITE(28,12) FNAME 
12 FORMAT(' RETIREMENT BENEFIT COSTS UNDER PLAN ',A5,' ($0900''S)') 
WRITE(28,13) 
13 FORMAT(' YEAR MALES FEMALES 
& TOTAL') 
WRITE(28,14) 
14 FORMATC! mame en 
& -ee-- *) 
DO 380 J=5,60 
380 WRITE(28,17) YR(J), (RTCST(U,K),K=1,3) 
GO TO 998 
381 WRITE(28,998) (RTCST(J,3),JU=5,60) 
998 FORMAT(5(3PF1G.0)) 
17 FORMAT (HX, 14 /8X,F12.1/5X,F12.1,5X,F13.1) 
h1 FORMAT(6(F10.0)) 
T=CLOCK(0,0)*60. 
TYPE 421,FILECIK+10),FILECIK+15),FILECIK+20),T 


421 FORMAT(//," END OF RUN',//," OUTPUT FILES ARE: ',3(A5,'.DAT, 
& ')//," EXECUTION TIME = ',F8.3,' SECONDS',///) 


42 FORMAT(6(F10.4,1X)) 
STOP 
END 
SUDROUTINE MIN 


THIS SUBROUTINE CALCULATES THE EARNINGS DROPPED OUT 
TO FIND TOTAL PENSIONABLE EARHINGS 


C2 <3.) <3 


COMMON /1ODATA/INDIV(-12/11,47, 2), EXPOS (-12/11),DROP(-12/11, 2), 
& NR,DRPRT 
REAL DUP(-12/11, 47, 2) 
NO 195 K=1,2 
DO 195 J=1,47 
pO 195 !=-12,11 
DROP(I,K)=0. 
195 DUP(I,J,K)=INCIV(I,J,K) 
DO 230 K=1,2 
LECK. EQ. 2. AND.DRPRT.LE..25)G0 TO 230 
C THIS CONTROL IS USED TO DROP OUT EARNINGS FOR FEMALES IF THE 
C DROP OUT RATE EXCEEDS 25% 
DO 230 I=NR-10,NR+1 
l=! 
NB=EXPOS(11)/DRPRT 
IECNR.EQ.9.AND.NB.LT.10)GO TO 230 
LL=(1FIXCEXPOS(11)))+1 
DO 220 L=1,LL 
JCNT=1 a. 
Ont epUPtil, CNT, K) 
IE(DUPC(I1,u,K).LT.SMALL)GO TO 200 
Go TO 210 
200 JCNT=u 
SHALL=DUP(I1,J,K) 
210 CONTINUE 7 
DUP(I1,UCNT,K)=DUP(I1,JCNT,K)+5000. 
220 DROP(I1,K)=DROP(I1,K)+#SMALL 
DROP(II,K)=DROPCII,K)= (C1. -CEXPOS(11)=CLFIXCEXPOS(11)))))*SHALL) 
230 CONTINUE 
RETURN 


END 
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VERSIOt! AS OF 12:00 AUGUST 11,1972 

PROGRAIHiE CONT.5.THIS 1S THE FIFTH PROGRAMME IN THE CPP PROCRAITIES 
CONT.5 CALCULATES THE CONTRIBUTIONS TO THE CANADA PENSION PLAI FROM 
196G TO 2025 UNDER VARIOUS USER SPECIFIED OPTIONS. TVESE OPTIONS 
INCLUDE REASONABLE STABILITY AND MODERATE INFLATION UNDER THE CUR- 
RENT AND PROPOSED SCIiEMES. 

PROGRAMME WRITTEN BY T.LETT AND G.JACOBS FOR THE TAXATION AND 
FISCAL POLICY BRANCH MARCH 1972. 


aagagaanangnada 


DOUBLE PRECISION POP(19,60,2),TEMP(11,60,2) 
REAL YMPE(60), YBE(GO0),MMAE(11,60, 2), ETADJ(60), 

+ PART(11,10),ADJ(60),FILE(15),T(11,7), CONRTE(60) 
INTEGER YEAR(60) 


YMPE AND YBE ARE SELF EXPLANATORY. POP CONTAINS THE DATA CEN- 

RATED BY PROGRAMME POP.1.USERS SHOULD BE CAUTIONED THAT THIS ARRAY 
1S OVERWRITTEN SEVERAL TIMES DURING CONTRIBUTION CALCULATIONS AS 

IS THE TEMP ARRAY TO CONSERVE CORE. THE TEMP ARRAY IS FOR TEMPORARY 
COMPUTATIONAL USAGE ONLY. THE MMAE ARRAY IS THE OUTPUT FROM PROG- 
RAMME MMAE.3.THE ARRAY PART CONTAINS TEN SETS OF ELEVEN EACH OF 
PARTICIPATION RATES FOR THE LABOUR FORCE. TWO SETS ARE FOR MALES 

AND THE REMAINING EIGHT ARE FOR FEMALES. THE VECTOR ADJ CONTAINS 
ADJUSTMENT FACTORS(NOT THE SAME AS IN MMAE.3) WHICH ARE FOUND Oli 
WORKSHEET E£108-5 .THE VECTOR FILE CONTAINS THE VARIOUS PROGRAMME OPTIOK 
CODES. TXFPB IS A TAXATION AND FISCAL POLICY OPTION THAT IS SPECIF- 
[ED BEFOREHAND IN OTHER PROGRAMMES. THE ARRAY T CONTAINS THE NEC 
ESSARY COMPONENTS FOR PRINTING OUT TITLES. THE VECTOR ETADJ CONTAINS 
THE EARNINGS TEST ADJUSTMENT FACTORS. THESE ARE FOUND ON 

WORKSHHEET E108-1. 


RAS PENS CO OO? Ce Gr CTO: acs Coe 


DATA ETADJ/7¥*1.,.933,.867,.8,.68,.56, 48, .32, 46%. 2/ 
DATA PART/.65, 3*1., 3*.93, 2*.87,.78,.53,.68, 3*1.,35%. 94,2, 88, 
+779) 15h, .46,.66,.44,6%. 36,.25,.13,.5,.69,.5,6%, 42, .29,,15, 055, 
+ .73,.57,6%.48,.34,.16,.58,.75,.58,6%.49,.35,.14,.46,.66, 44, 


- O64, 356, .25, 213,089, .68,.5,6%, 41, .29,.13,.53). 71, .56,0%. 975 
wl 963,024, 555575, 257, 5%, 88,.58,.14/ 


DATA ADJ/ 
2 58.1 pe 090,45 O92, Se Se 3 DR ORG DC TRE nt ae 
DATAC (TC I,J), J=1, 7), 1, 11)/ 


'REASO!, ARC LER STAR MTL TV"|.* CURRT, TENT.P*, CARS, 
"ODER", "ATE 1, 'NFLAT', "ION C', "URREN', 'T PLA", 'N ey 
'REASO', "'NABLE', ' STAB’, "IETTY', © PROP"? *OSED. °,.* PEAR 
'MODER', "ATE I1', "NFLAT', "ION P", 'ROPOS', "ED PL", ‘AH : 
TUSER to 'SpEC! "FIED t, TOPTIO®, RT, * : 
iwSE ©. VEQUAL®, 9S. FLA), "T RAT‘, (EOF |? 
'yBE ',"propo','"RTION',' OF Y', 'MPE 1! j ; 
TYEAR ', "CONTR', 'IBUTI', ‘ONS ', 'TOTAL', 'aRWTH',' ', 
MALES | 'FEMAL', TES Pelee el wee, 
tae ye ER ie , th 4 4 

4 4 , 4 


ON R', "ATE E', 


7 
tc € FF F 


. fk. 8 & 


t+eeeeeeeeet 


i 
QUALS',' '/ 


Rogar s RUaY, 


DATA FILE/' RSCUR', "MICUR', "RSPRO', "MIPRO', "TXFPB', 'MMRSC', 'MMMIC' 
THERSP', "MMtl P!, "MMTEP?, "RSCCN', 'MICCN', 'RSPCN', 'NIPCN', 'TEPCN'/ 
DATA’ N/1/ 
102 TYPE100 
100 FORMAT(' TYPE It PLAN TO BE TESTED: ',$) 
ACCEPT2, FNAME 
2 RORIARE RS? 
DO 3 IK=1, 
3 J ECENAME. a4 FILECIK)) GO TO 4 
TYPE103, FNAME 
103 FORMAT(1X,A5,'?, TRY AGAIN'/) 
GO TO 102 


: DATA FILE OPSHN CONTAINS THE OPTIONS THAT CAN TESTED. X1YECE 1S THE 


C ELAT RATE OPTION FOR THE YEE.NOTE THAT IT 1S NECESSARY Te LEAVE 
© THIS VALUE AT G00. IF THE USER IS TESTING A PROPORTIONAL VALUE FOP 
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RATE WITH ANY OF THE CURRENT OR PROPOSED SC 
NOTE THAT THE SWITCH MUST CONFORM WITH THE 


BEING USED X2YBE MUST BE .12.CONRTE IS TRE 
[IE .036. 


a 3 CD CY C3 69 S39 Cee Ce 2 


CALL IFILE(20,FILE(IK+5)) 
CALL IFILE(21, FNAME) 
CALL IFILE(22, 'ONE’) 
CALL IFILE(23,'OPSHN') 
READ(20,5)(C((MMAECI,J,K) 
READ(21,6)(YMPE(1),1=1,6 
READ(22,6)(((POP(I,uJ,K), 
READ(23,5) X1YBE,KX1, X2Y 
READ(23,6)(CONRTE(J),J=1,60) 
FORMAT(5G) 

FORMAT(10G) 


2 
6 
7 FORMAT(9G) 
9 


FOR THE YBE.KX1 1S A SWITCH WHICH ALLOWS Th 


,J=1,60),1#1,11),K= 

0) 

J=1,60), 1=1,19),K=1,2) 
BE 

70 


E USER TO USE THE FLAT 
HEMES. THE USER SHOULD 


REFER TO STATEMENT NO. 999 TO SEE THE MECHANICS OF THIS SWITCH. 


STRUCTURE 


OF THE VECTOR FILE.X2YBE IS THE OPTION WHICH ALLOWS THE YBE TO 
EXPRESSED AS A PROPORTION OF THE YMPE.IF THE FLAT RATE OPTION IS 


TOTAL CONTRIBUTION 


RATE(EMPLOYER PLUS EMPLOYEE) AND MUST BE EXPRESSED AS A DECIIAL, 


1,2) 


99 IF (FNAME. EQ.FILE(KX1).0R.FNAME.EQ.FILE(KX1+1)) GO TO 91 


DO 8 1=1,60 
8 YBECI)=X1YBE 
GO TO 10 


91 DO 9 I1=1,60 


9 YBECI)=(1FIXCX2YBE*YMPE(I )+1)/100)*100 
C 
C THIS SECTION FIRST CALCULATES THE PARTICIPATING POPULATION AND 
C THEN COPIES IT INTO TEMP FOR LATER USE.THE PARTICIPATING POPUL~ 
C ATION FIGURES STORED IN POP ARE SUBSEQUENTLY MULTIPLIED BY 
C MODIFIED MODIFIED AVERAGE EARNINGS TO OBTAIN AGE GROUPED EARNINGS 
C FOR THE WORKING POPULATION. 
C 
GOV FOTCS1 57,51, 52, 559,7%K 
51 IMl=4 
GO TO 56 
52 IMI=0 
GO TO 56 
55 TY PES4 
54 FORMAT(1X, ECONOMIC ASSUMPTIONS=-MI/RS?',$) 
ACCEPT55,ANS 


55 FORMAT(A2) 
IE(ANS.EQ.MI) GO TO 52 
GO TO 51 
56 JJ=0 
10 DO 101 K=1,2 
DO 11 1!#1,11 
DO 11 J=1,60 
N=1 
1F(K. EQ. 1.AND.JU.GE.20) N=2 
1F(K.EQ.2.AND.U.GE.1) N=3 
IE (K.EQ.2.AND.JU.GE.10) N= 
iF (K.EQ.2.AND.JU.GE.20) N=5 
1F(K.EQ.2.AND.U.GE.40) N56 
POP(1,U,K)=POP(1+3,u,K)*PART(1,JU+N) 
TEMP(1,u,K)=POP(1,d,K) 
11 POP(1,JU,K)=POP(1,JU,K)*MMAE(I,J,K) 
101 JJ=Ju+IMI 


THE FIRST AGE GROUP IS SCALED BY .59 BEFORE 
BY .& AFTER THAT DATE. THIS OPERATION COMPEN 
THIS AGE GROUP CONTAINS INDIVIDUALS FROM 15 


QadONND 


DO 14 K=1,2 
‘po 14 J=1,60 

1F(u.GE.34) GO TO 120 
POP(1,uU,K)=POP(1,U,K)*. 39 


THE YEAR 2000 AND 
SATES FOR THE FACT THAT 
TO 19 IN ACE WHILE 


CONTRIBUTIONS TO THE CPP DO NOT START UNTIL AGE 18. 
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GO TO 14 | 
120 POP(1,J,K)=POP(1,JU,K)*.4 
14 CONTINUE 


AT THIS POINT THE GROUPED EARNINGS ARE TOTALLED BY YEAR AND 
MULTIPLIED BY AN INDEX TO REFLECT THE GROWTH IN THE YHPE THRU 
TIME.THE TEMP ARRAY IS ALSO TOTALLED. THE CONSTANT 5000 IS THE BASE 
YEAR(1966) YMPE,. 


Cd Cd. CS C3 CD 


DO 15 K=1,2 

DO 15 J=1,60 

POP(11,JU,K)=POP(11,JU,K)*ETADU(J) 

DO’ I> t=2,11 

POP(1,JU,K)=POP(1,JU,K)+(POP(I,JU,K)*(YMPE(J)/5000)) 
15 TEMP(1,JU,K)=TEMP(1,JU,K)+TEMP(1,J,K) 


IN THE NEXT SECTION 1 1S ADDED TO THE ADJUSTMENT FACTOR SO THAT 
IT SCALES THE DATA OVER 100%.THE YBE 1S THEN MULTIPLIED BY THE 

PARTICIPATING POPULATION AND THEN IS SUBTRACTED FROM THE TOTAL 

EARNINGS TO ELIMINATE THE TOTAL EXEMPTIONS FROM EARNINGS. 

TOTAL EARNINGS ARE THEN ADJUSTED AND MULTIPLIED BY THE CONTRB- 

UTION RATE TO GET TOTAL CONTRIBUTIONS. 


CID 3.033 Co. Ges 


DO 60 J=1,60 
60 ADJCJ)=ADUCJ) +1. 
OC Ie ease 
DO 16 J=1,60 
16 POP(1,JU,K)=(POP(1,JU,K)-CYBE( JU) *TEMP(1,J5,K))) *ADJUC JU) *CONRTE CU) 
DO 13 K=1,2 
DO 13 J=1,60 
YEAR(CJ)=1965+d 
13 TEMP(1,J,K)=0. 


C 
C CONTRIBUTIONS FOR MALES AND FEMALES ARE TOTALLED INTO TEMP 
C FOR PRINTING PURPOSES. THE SECOND DIMENSION OF TEMP CONTAINS 
C THE GROWTH RATE OF TOTAL CONTRIBUTIONS. 
C 

DO 18 K=1,2 


DO 18 J=1,60 
18 TEMP(1,J,1)=TEMP(1,J,1)+POP(1,JU,K) 
DO 19 J=1,60 
L=J+1 
IFCL.EQ.61) GO TO 19 
~ TEMP(1,L,2)=¢((TEMPC1,L,1)-TEMPC1,JU,1))/TEMP(1,J,1))*100. 
19 CONTINUE 
IVPES2 
32 FORMAT(1X, "OUTPUT? (TABLE/DATA):',$) 
ACCEPT50,ANS 
IFCANS.EQ.'DA') GO TO 33 


AeEobR ROAM FHE TNENACT SEPALS CORVERS [On THAT FARES BEACE 1 
IT IS LEFT AS A DECIMAL FRACTION. THIS OPERATION IS NECESSARY SO 
THAT THE LOGICAL COMPARISONS WILL WORK, 


C3 co 9 C309 ' CE 


X2YBE=X2YBE*100. 

IFC(X1YBE.EQ.600.) IJ=7 

IF(X2YBE.EQ.12.) IJ=6 

IF (X1YBE.EQ.600. .AND.KX1.EQ.1.AND. X2Y8E.EQ.12.) GO TO (24,24, 


+ 25,25, 26), 1K 
IFC IU-6)22,22, 25 
22 XYBE=X1YBE 
GO TO 28 
23 XYBE=X2YBE/100. 
GO TO 28 
24 XYBE=X2YBE/100. 
| J=7 
GO TO 28 
25 XYBE=X1YBC 
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— 1J=6 


GO TO 28 
296 XYBE=X1YBE 
1J=6 
TYPE27,XYBE 
27 FORMAT(" x*e* WARNING YBE HAS BEEN SET TO',F6.2, '***') 


28 CALL OFILE(24, 'CONT') 
CALL TIME(XTIME) 
CALL TDATE(LD, LM, LY) 
WRITE(24, 91) LDO,LM,LY,XTIME 


21 FORMAT (! 1', 50X, "DATE OF RUN' 1X, RnR yy , "TIME 'AS) 
WRITE(24, 17)(TUIK, 1), f=1,7), (TC1u, | ), | =1,6), XYBE, (T(8,1), 
+ |=1, 4), (T(8,1), 1=2, 6), T(10, 4), (T09,1),1=1,3),T(8, 5), T(10,1), 
+ T(10, 1), T(10, 1), T(10, 2), T(10, 1), TeTG, £iy LGAes Bee T(10, 1), (10, 1) 
17 FORMAT(' ,50X7A5//' ihe 50X6A5, P72, 7/1°0", 4OX,A5, 3X, 


+ 2(3A5, 2X), LXA5, 6X, 2A5/' ’52XA5, 11X2A5, 21XA5/' ,40X, 
+ A5, 3x2 (3A5, 1X), 6XA5, BXA5). 
nO 29 J=#1,60 
WRITEC24, 1)CONRTE(U), YEAR(J), POP(1,uU, 1), POP(1, ud, 2), TEMP(1,u,1) 
+ STEMP(L, Jg2s 
29 CONTINUE 
1 FORMAT('O', 33X, 2PF3.1, ' ¥"3X14, 3X3PF13.0, 2(4X3PF13. 0), 5SXOPF5.2) 
50 FORMAT(A2) 
33 |FCANS.EQ.'TA') STOP 
CALL OFILE(25,FILE(IK+10) ) 
WRITE(25, 31) (TEMPCL, iLL SL; Ge 
31 FORMAT (! , 5(3PF14.0, 1X)) 
TYPE34, Saver ety) 
34 FORMAT(1X, "OUTPUT RYEE<(S “AG ) UAT" 
STOP 
END 
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C VERSION AS OF AUGUST 11,1972 

C PROGRAMME WIDOW.6 

C WRITTEN BY HARRY NEWTON. APRIL 1972. 
COMMON/BLOCK1/ PWEWB (58, 10),MAXPEN(96) 
COMMON /BLOCK2/COEFF(58,16),STOTAL(16) 
REAL HAXPEN 
DIMENSION SURVIV(7,16,2),RMARRY(5,10, 2) 
DIMENSION POP(19,60),SURVM(8,15),FAC68(10) 
DIMENSION YMPE(60), PNIND(60),ESC1(58),DTHFC(15,60),FFAC(58) 
DIMENSION FNAME(5),FOUT(5), DNAME (5) 
DIMENSION LIMITICK), LIMIT2(7), OAS (3),ERPCTUC(3), ERPCTO(3) 
DIMENSION PRPMRD(15, 5), AGEDIS(15, 6), NHAGRP(15), NWAGRP(15) 
DIMENSION ESCDBF (15), DEATHS(15),ZERO( 944 ) 
DIMENSION ERBU65(58), ERBO65(58),FRBU65(58) 
DATA ERBU65, ERBOG5, FRBU65/174*0./, |IKOUNT/5/ 
DATA FAC68/600.,793.,881., 864.,2*855.,832.,785.,699.,405./ 


DATA FFAC/5#0.,23.,32.,61.,64.,88.,111.,155.,158.,181.,204.,22 


* 25154 4190,7 
DATA OAS /300.,5960., 6. /, ERPCTU/. 375, .75,0./, ERPCTO/.6,.75,0./ 
DATA FNAME/"RSCUR', 'MICUR', 'RSPRO', 'MIPRO!, 'TXFPB'/ 
DATA DNAME/'DHRSC', 'DHMIC', 'DHRSP', 'DHMIP', 'DHTFP'/ 
DATA FOUT/'RSCWD', 'MICWD', "RSPWD', "MI PWD', 'TFPWD'/ 
DATA WNHAGRP/3,4,4,5,5,9*6,5/ 
Reh rr tecee ces Wosek ociuacrone wee 1a 
tA ? vA 
DATA LIMI1T2/1970, 1975, 1980,1985,1990,1995,2000/ 


DATA ((PRPMRD(1,J),J=1,5),1=1,15)/5*.3,5%.7,5%. 82, 5%. 86, 9*,. 88, 


eat 


89, 


2*.87,3%.88,.84,2*. 85,2. 86,.82,.83,.84,2%.85,.8,.81,.82,.83, 


Z 
2 38h, 743. 76,.78,.8,.81,.67,.7, 727s 1, < 16,255, .58,55,.02, 259, 
3 


4565548705) . 525,31, 3 3471 3671 505547 


DATA (C(AGEDIS(1,J),U#1,6),1=1,15)/.22,.69,.09, 3*0.,.03,.4,.49, 
36, 
06, 
OF, 
. 36, 


90005 1 O9S4485 34565 095 280. .02; 6.49.38, 42, .08, 0.5503, 13205; 
SOR het OSS 08, 615, 38,2 3%,.075 002,205, . 07, 238, stp 00h 0S; 
PRS K S¥6 5.554 067 5.03)407) 5.1956 365428,...07,.05,.07,. 18, 2547528, 
OBL O7% SOs) Shyh2eG 06,07, +075: 215. 35,.25;,'.05; (08, 208,25, 
O62 O3 7005 1 15. C6yaSSyi07,,02,445,5105,483), 32,422, 85/ 

EQUIVALENCE (COEFF(1,1),ZERO(1)) 

CALL CLOCK1 

34 TYPE 2001 
2001 FORMAT(//,' TYPE IN PLAN TO BE TESTED - ',$) 
ACCEPT 2002, PLAN 
2002 FORMAT(AS) 
DO 40 J=1,5 
IF (PLAN. EQ.FNAME(J)) GO TO 41 
40 CONTINUE 
TYPE 2003, PLAN 
2003 FORMAT(//,' ***ERROR*** NO PLAN CALLED ',A5) 
GO TO 34 
41 KFILE=d 
CALL I FILE(28, PLAN) 
¢ READ IN MAXIMUM PENSIONABLE EARNINGS FOR EACH YEARCYMPE). 
READ( 28,2004) (YMPE(J),J=1,60) 
: GO TO (43,44,43,44,45),KFILE 
45 TYPE 2005 
2005 FORMAT (' He IN GROWTH RATE FOR AVERAGE EARNINGS AFTER 1975', 
* /,3X,'C1.E., FOR 5.5%, TYPE 1.055) === ',$) 
ACCEPT 2004,RATE 
CALL 1FILE(27, 'TEPNX') 
GO TO 46 
43 CALL 1FILE(27,'RSPNX') 
RATE=1.035 
GO TO 46 
bh CALL IFILE(27, "M1 PNX") 


wityrt "200 ) (PNIND(J), 0) 
4 
C *CARCULATE’ MAXIMUM, PENSION OR’ EACH YEAR (IMAXPEN). 
DO 42 J=1,58 
b2 MAXPENCJ )=(YMPE (J) +YMPECU+1)+YMPE(d+2)) #100. /(12. #PNIND(J+2)) 
DO 39 J=59,96 
39 MAXPEN(J)=\MAXPEN(J=1) *RATE 
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xix 


TYPE 1006 
1006 FORNAT(' TYPE IN RETIREMENT PENSION RATE (0-1) -- ',3) 
ACCEPT 2004, RETRT 
CALCULATE ESCALATION COEFFICIENTS. 
NO 47 J=1,58 
47 ESC1(J)=RETRT*(YMPE(J )+YMPE(U+1)+YMPE(U+2))/CPNIND(U+2)*37. 5) 
CALL IFILE(29, 'ONE') 
READ(29,2004) ((POP(I,J),J=1,60),1=1,19) 
CALL IFILE(29, 'SURVM') 
READ(29,2010) ((SURVM( I,J), 11,8), J=1,15) 
2010 FORMAT(8G) 
CALL 1 FILE(29,DNAIHE(KFILE)) 
READ(29,2008) ((DTHFC(K,L),L=1,60),K=1,15) 
2008 FORMAT(6G) 
MALE DEATH BENEFIT FACTORS FOR 1968-9 ARE NOT COMPUTED IN THE 
PROGRAM 'RETBN.4'. HENCE, THE VALUES SUPPLIED BY THE DEPT. OF 
INSURANCE, OTTAWA ARE INSERTED. 
DO 51 K=1,10 
DTHFC(K,3)=FAC68(K) 
S1 DTHFC(K,4)=FAC68(K)+(DTHFC(K, 5)-FACE8(K) )/2. 


IN THIS SECTION, THE FEMALE DEATH BENEFIT FACTORS COMPUTED IN THE 
PROGRAM 'RETBN.4&' ARE ADJUSTED IN THE SHORT TERM FOR THE PURPOSE OF 
COMPUTING REDUCTION FACTORS AT A LATER STAGE IN THIS PROGRAM. 
SINCE THE FEMALE DEATH BENEFIT FACTORS ARE ALMOST IDENTICAL It THE 
SHORT TERM FOR BOTH REASONABLE STABILITY AND MODERATE INFLATION, 
THE PROGRAM DOES NOT DISTINGUISH BETWEEN THESE TWO TYPES OF SCHEMES 
WHEN ADJUSTING THE FACTORS. 

DO 48 Ke1,11 
48 READ(29,2008) (ZERO(L),L=1,19), (FFAC(L),L=18, 58) 
DO 49 K=#8,17 
49 FFAC(K)=FFAC(K)*ESC1(K-7) 
DO 50 K=18,58 
50 FFAC(K)=FFAC(K) *ESC1(K-7)/ESC1(K=-5) 
CALL IFILE(29, 'REMAR') 
READ(29, 2004) (((RMARRY(K,L,M),L=1,10),K=1,5),M=1, 2) 
2004 FORMAT(10G) 
CALL IFILE(29,'SURVP') 
READ(29,1002) (((SURVIV(K,L,M),L=1,16),K=1,7),M=1,2) 
1002 FORMAT(16F4.3) 
CALL IFILE(29,'PWEWB") 
READ(29,1003) ((PWEWB(K,L),L=1,10),K=1,58) 


1003 FORMAT(5X,10F5.4) 


TYPE 1001 

1001 FORMAT(' TYPE IN STARTING YEAR (DEFAULT IS 1968) - ',$) 
ACCEPT 2004, 1YRST 
IFCLYRST.LT.1968.0R.1YRST.GT.2025) IYRST=1968 


TYPE 1004 
1004 FORMAT(' TYPE IN PROPOSED ANNUAL FLAT RATE PAYMENT, AND THE', 
* /,6X,"YEAR IT BEGINS --- ',$) 
ACCEPT 2004,0AS(3),1FR 
| FR=1FR-1967 


OAS (2)=100. *0AS(2)/PNIND(8) 
OAS (3)=100. *OAS(3)/PNIND(IFR#2) 
TYPE 1005 
1005 FORMAT(' TYPE IN PERCENTAGE (BETWEEN 0 AND 1) OF BENEFIT THAT 
* "1S EARNINGS',/,5X,' RELATED (1) FOR WIDOWS UNDER 65 --- ',$) 
ACCEPT 2004, ERPCTU(3) 
TYPE 1008 
1008 FORMAT(14X,'(2) FOR WIDOWS 65 AND OVER --- ',$) 
ACCEPT 2004, ERPCTO(3) 
TYPE 1009 
1009 FORMAT(14X, "BEGINNING IN THE YEAR ',$) 
ACCEPT 2004, 1ER 
[ER=IER-1967 
DO 100 IYR=IYRST, 2025 
DO 1 J=1,4 
IFCIYR.LT-LIMITI(J)) GO TO 2 
1 CONTINUE 
J=5 
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38 
37 
29 

7 


| X=J 
DO 3 1=#1,944 
ZERO(1)=0. 


NX=10+(1YR-1965)/5 

NX=IA1NOCNX,15) 

1YEAR=1YR-1967 

LXO=(1YEAR+7)/5 

DO 38 J=1,7 

IFCIYR.LT.LIMIT2(U)) GO TO 37 

CONTINUE 

J=8 

1 Xl=dJ 

DO 35 11,15 

DEATHS( 1 )=(1.-SURVM(1X1,1))*POP(I1 +4, | YEAR+2) 
DO 7 J=1,9 
ESCDBF(J)=ESC1( 1 YEAR) *DTHFC(J, | YEAR+2) 


DEATH BENEFIT FACTORS FOR AGE GROUPS 65+ HAVE ALREADY SEEN 
ESCALATED IN THE PROGRAM 'RETBN.4' 


12 


wo 


10 


13 
11 


DO 12 J=10,15 
ESCDBF(J)=DTHFC (J, | YEAR+2) 

DO 5 1=1,NX 

NUM=NHAGRPCI ) 

NST=NWAGR PCI) 

NF IN=NST+NUM-1 

DO 6 J=NST,NFIN 

TEMP= =DEATHS(1)*PRPMRD(1, 1X) *AGEDIS(I, J-NST+1)*1000 
COEFF(I YEAR, J)=COEFFCIYEAR, J)+TEMP 
STOTAL(J) =STOTAL(J)+TEMP*ESCDBF (1) 
CONTINUE 

DO 4 1=1,16 

TEMP=COEFF(IYEAR, | ) 

STOTAL(1 )=STOTAL(I)/ANAX1(1., TEMP) 
M=90 

'41=0 

HST=LYR-1966 

[FCNST.GT.58) GO TO 69 

H2=1 

[FC LYR.GT.1969) M2=2 

DO 8 I=NST,58 
IFCL.NE.NST.AND.I.GT.3.AND.MOD(1+2,5).GT.0) GO TO 8 
KYR=LYRtI-NST 

DO 9 J=1,6 

IFCKYR.LE.LIMIT2(J)) GO TO 10 
CONTINUE 

J=7 

| X=J 

LX=1+6 
1F(MOD(LX,5).EQ.4.AND.LX.GT. (28+NST)) M1=M1+1 
LX=LX/5 

IF CLXO.NE.LX) M=M+(LX-LX0) 
UX=14(1+1)/5- (NST#1)/5 
JX=MINOCJX, 5) 

DO 11 J=1,16 

IF((J+i).GT.16) GO TO 11 
JX1=M1N0(10,JU+M1) 

IFCCI-NST).GT.3) GO TO 
COEFF(I,JU)=COEFF(I YEAR, 15 #RMARRY( JX, JX1,M2)#SURVI VIX, J#ti, M2) 
GO TO 11 
COEFF(1,JU)=COEFF(1-5,U)*RMARRY (JX, JX1, M2) *SURVIV(JX,dth,M2) 
CONTINUE 

LXO=LX 

CONTINUE 

NSTI=NST+#1 

NST2=5-MOD(NST+2, 5) 

1F(NST2.EQ.5) GO TO G2 
NST3=NST+#+NST2-1 

STi=1. 

ST2=NST2+1 

NO 60 1=NST,NST3 

DO 61 J=1,16 
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61 COEFF(1,JU)=(COEFF(NST=1, J) *(ST2-ST1)/ST2)+COEFF(NST341, J)*ST1/ST2 

60 ST1=ST1+tl. 

66 NST1L=NST3+2 

62 IF(NST1.GT.58) GO TO 65 
$T1=1. 
DO 63 I=NST1, 58 
IF(MOD(1+2,5).EQ.0) GO TO 64 
N2=34+5%((1+2)/5) 
N1=SN2-5 
DO 67 J=1,16 

57 COEFF(I,Ju)=(COEFF(N1,J)*(5.-ST1)/5. )+COEFF(N2,J)*STI/5. 
GO TO 63 

64 ST1=0. 

63 ST1=ST1+1. 

65 IFCIYR.GT.1968) GO TO 69 
DO 68 J=1,16 

68 COEFF(2,JU)=(COEFF(1,JU)+COEFF(3,U))/2. 

69 CONTINUE 
KX3=10+(1 YEAR+2)/5 
KX3=MAX0(KX3,11) 
KX3=MINO(KX3,16) 
DO 17 K=I YEAR, 58 
JER=1 
IFCKFILE.GT.2.AND.K.GT.5) JER=2 
IFCK.GT. (I ER-1)) JER#=3 
JFR=1 
|FCKFILE.GT.2.AND.K.GT.5) JFR#=2 
IFCK.GT.(CIFR-1)) JFR=3 
H=1. 
IFCK.EQ.1 YEAR) H=.5 
KX=(K+4-1YEAR)/5 
XK=1. -MOD(K=I1 YEAR, 5)/5. 
IF(XK.GT..9) XK=#0. 
KX=10-KX 
KX1=KX+1 
IFCKX) 25,25, 23 

23 DO 16 J#1,KX 

16 ERBU65(K)=ERBUG5(K)+ERPCTU( JER) *STOTAL (J) *COEFF(K,J)*H 
ERBU65(K) =ERBUG5(K)+ERPCTU( JER) *STOTAL(KX1) *COEFF (KAL, J) *KK FE 

25 KX1=MAX0(1,KX1) 
DO 18 J=KX1,KX3 
1YR72=60+1 YEAR-5*J 
I YR72=MINO(58, 1YR72) 
1YR72=MAX0(6, |YR72) 
R=REDFAC(u, 1 YEAR, STOTAL(J), FFAC(IYR72)) 

18 ERBO65(K)=ERBOG5(K)+ERPCTO( JER) *STOTAL (J) *COEFF(K,JU)*R*H 
IFCKX) 17,22, 26 
FFOXKSLT..1)7°>GO°TO 29 

26 1YR72=60+1YEAR-5*KX1 
1YR72=MINO(58, 1YR72) 
1YR72=MAX0(6,1YR72) 
R=REDFAC(KX1, | YEAR, STOTAL(KX1), FFAC(1YR72)) 
ERBO65(K)=ERBOG5(K)-ERPCTO(JER)*STOTAL(KX1) *COEFF (KX1, JU) *R*XK eh 

29 KX4=KX3+1 
IF (KX4.GT.16) GO TO 24 
DO 19 J=KxX4,16 

19 ERBO65(K)=ERBO65(K)+ERPCTO(JER)*STOTAL (J) *COEFF(K,J)*H 

24 DO 20 J=1,KX 
KX2=MAX0(1,JU-6) 

20 ERBUG5(K)=FRBU65(K)+COEFF(K, J) *PWEWB(I YEAR, KX2) *OAS( UFR) #E 

22 KX2=MAX0(1,KX1-6) 
IF(K.EQ.1 YEAR) GO TO 17 
ERBUGS(K)=FRBUG5(K)+COEFF(K,KX1) *PWEWB( 1 YEAR, KX2) #OAS (UFR) *XK*H 

17 CONTINUE 

100 CONTINUE 
DO 21 1=1,58 
ERBU65(1 )=(ERBU65(1)+ERBO65(1))*.009*PNINDC(1 +2) 

91 FRBU65(1)=FRBUG5(1)*. 009*PNIND(1+2) 
CALL OF! LE(23,FOUT(KFILE)) 
WRITE(23, 2010) ERBU65,FRBU65 
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T=CLOCK(0,0)*60. 
TYPE 2009,FOUT(KFILE),T 
2009 FORMAT(//,' END OF RUN',//,' OUTPUT FILE IS AS, ‘DAI’, //, 
* ' EXECUTION TIME = ',F8.3,' SECONDS',///) 
STOP 
END 


FUNCTION REDFAC(AGEGRP, | YR,MFAC, FFAC) 

THIS SUBROUTINE COMPUTES REDUCTION FACTORS. 

WRITTEN BY HARRY NEWTON. APRIL 1972. 
COMMON/BLOCK1/ PWEWB(58, 10), MAXPEN(96) 
INTEGER AGEGRP, YEAR, YR65,AGE 
REAL MAXPEN,MFAC 
DIMENSION PCT(5),DIST(5),BENM(2),BENF(5), PWNRP1(11), PWNRP2(11) 
DATA PCT/.225,.125,.075,.075,0./,DIST/.225,.125,.075,.075,9./ 
DATA PWNRP1/.55,.83,.78,.73,.67,-61,.54, .47,.4,.32,.15/ 
DATA PWNRP2/.85,.5,.75, «69, +63,.57,.5,.43,.36,.25,.11/ 

1 K1=AGEGRP-6 

K1=HAX0(K1, 1) 
X}iXPEN=MAXPENC1I YR) 
AVGM=MFAC/PWEWBCIYR, K1) 
IFCIYR.GT.2) GO TO 6 
PWNRP=PWNRP1(AGEGRP) 
GO TO 7 

6 1OLD=(1YR+4)/5 
1OLD=MINO(IOLD, AGEGRP-1) 
PWNRP=PWNRP2 (AGEGRP=-1 OLD) 

7 AVGF=FFAC/PWNRP 

AGE=10+5*AGEGRP 
IFC(AGE.GT.64) GO TO 25 
[FCAVCM. GT. (.75*XMXPEN)) GO TO 24 
BENM(1)=1. 25*AVGIM 
GO TO 26 

24 BENM(C1)=XMXAPEN 

26 BENIM(2)=2*(AVGM-BENM(1)/2.) 
GO TO 27 

25 BENM(2)=. 75*AVGM 
BENM(1)=2*(AVGM-BENM(2)/2.) 

27 YR65=63+1YR-AGE 


IFCYR65.LT.1). YR65=1 
LESYRER: Gt: 9) TYPE 28, 1YR,AGE,YR65 


98 FORMAT(' YEAR =',15,8X, "AGE =',14,8X,'YR65 =", 15) 
YEPEG6 5=MAXPEN( ¥R65): ; 
IECYR65.LT.9) YMPEGS=YMPE65*(FLOAT(YR65)+1.5)/10. 
BENF (2) =AVGF 
BENF(3)=1. 5*AVGF 
BENF (4 )=YMPE65 
IF (YMPE65.GT. (2. 3*AVGF)) GO TO 31 
BENF(5)=0. 
BTOT=. 25*BENF(2)+.15*BENF(3)+. 15*BENF (4) 
BENF(1)=(AVGF-BTOT)/.45 
GO TO 32 
31 BENF(1)=.4*AVGF 
OTOT=.45*BENF(1)+.25*BENF(2)+. 15*BENF(3)+. OS*BENF (4) 
BENE (5)=10*(AVGF-BTOT) 
32 CONTINUE 
BENM(1)=. 75*BENM(1) 
BENIM(2)=. 75*BENM(2) 
BTOT=0. 
PCT(4)=DIST(4) 
PCT(5)=DIST(5) 
IF(BENF(5).E0.0.) GO TO 10 
PCT(4)=.025 
PCT(5)=.05 
10 YMPEGS=MAXPEN(YR65) 
IECLYR.GT.YR65) YMPE65=XMXPEN 
DO 20 J=1,5 
DO 20 1=1,2 
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IE CCAMAXI(BENM(1),BENF (J) )+. 5*#AMINICBENM (1), BENF(U))). GT. YMPE65) 
GO TO 21 

IF(BENM(1).GT.BENF(J)) GO TO 22 

TEMP=DENMC(I)*.5 


GO TO 23 

TEMP=BENM(1 )-BENF(J)*.5 

GO TO 23 

TEMP=YMPE6 5-AMAX1(BENM(1),BENF(J)) 


TEMP=AMAX1(0., TEMP) 

BTOT=BTOT+TEMP*PCT (uJ) 

CONTINUE 

PROPB P=PWEWB(1YR,K1)*PWNRP 
PROPWP=PWEWB( I YR,K1)*#(1.-PWNRP) 

TEMP=, 75*AVGM 
REDFAC=(PROPBP*BTOT+PROPWP*TEMP) / (TEMP*PWEWECIYR, K1)) 
RETURN 

END 
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VERSION AS OF 12:00 AUGUST 17,1972 

PROGRAIIiE DISBN.7.THIS PROGRAMME CALCULATES DISABILITY BENEFITS 
FOR THE CANADA PENSION PLAN. IT ALSO COMPUTES WIVES' BENEFITS AS A 
SPECIFIED OPTION. 


€> > C2 €3 C2 


REAL POP(19,60,2),DTHFAC(15, 60,2), PRVRTE(9,6,2),DISINS(4), 
+ PROP(9,10),FILE(20),DISFR(9,60,2), YMPE(60),PNIND(60), 
+ WIFEBN(19,60, 2),0AS(3) 


REAL DISERN(9,60, 2), PHWHWC(10), FACTOR(9),FMTH65 (4,10) 

INTEGER YEAR(60),YR 

DATA (CC PRVRTEC(I,u,K),1=1,9),J=1,6),K=1,2) 
/.0k,.1,.24,.4,.66,.96,1.78,2.88,5.,.05,.13,.35,.6, 
-96,1.39,2.61,4.27, 7.64, .05,.13,.37,.67,1.07,1.55, 2.95,4. 78, 
8.66,.05,.13,.36,.69,1.12,1.63,3.07,4.97,9.01,.05,.15,.56,.68, 
1,.13,1.67, 3.15, 5.08,9.21,.05, #13,.37,+468, 1.12, 1.68, 3.16,5.12, 
931, 305, 12, o Spo ys Beg de dls Code Sods. 0s Lye Voge ad tee a0e 
.62,1.,1.45,2.73,4.46,8.09,.05,.14,.38,.68,1.08,1.57,2.96,4. 82, 
8. 76,.05,.14,.38,.7,1.13, 1.64, 3.08,4.98,9.02,.05,.14,.58,.7, 
1.15,1.67,3.15,5.08,9.21,.05,.18,.38,.7,1.15,1.69,3.18,5.13,9.3/ 
. FITHE5/.49,1.365,3.605,7.35,.731,2.421,6.532,14.112, .963, 
3.257,9.124,19. 846,1. 04, 3.655,10. 393, 23.189,1. 087,53. 844, 
11.064, 24. 881,1.109,3.942,11.3587, 25. 79, 1.12, 5.991, 11.552, 

OG. 223- 1.612, &, 211.593, 26. 36,1.12,4.,11.6, 26.392, 1.22,4%., 
11.6,26.4/ 


+ 


t+eetetet t+ 


DAT 


> 


++ + + 


DATA ((PROP(I,J),1=#1,9),J=1,10) 
+ /.4,3%.95, 3%*.9,2%.85,.4,3%.95, 3%.91,2*.86,.2,.45,.4, 
+ 5%.35,.3,.21,-52, -47,5%.42,.36,.23,.6,.54, 5%. 49,.43,.24,.61, 
$755.58. 5, Hh, 2p ob 5p oh, SR oy aap a clg ete nhl se OMe the wlesang 
Foe eae Oe Gs bas toe £07 + ep OR Noe Foe 
DATA PHWHWC/.576,.738,.942,.947,.894,.843,.741,.558, .356,.114/ 
DATA FACTOR/386.4,683.5, 791.7, 828. 2, 834. 7,835.6, 813. 2,777. 8, 
+ 696.6/ 
DATA OAS/300.,960.,0./ 
DATA DISINS/.95,.7,.75,.8/ 
DATA FILE/'RSCUR', 'MICUR', 'RSPRO', 'MIPRO', 'TXFPB', 'DHRSC', 
+ "DHMIC', 'DHRSP', 'DHMIP', 'DHTFP', 'DSRSC', 'DSMIC', 'DSRSP!, 
+ 'DSMIP', 'DSTFP', 'RSCWFE', "MICWF', "RSPWF', 'MI PWE', 'TFPWF'/ 
EQUI VALENCE( POP, WI FEBN) 
101 TYPE1 
1 FORMAT(' TYPE IN PLAN TO BE TESTED: '$) 
ACCEPT2, FNAHE 
2 FORMATC(AS) 
DO 3 IK=1,5 
3 |FCFNAME.EQ.FILE(IK)) GO TO & 
TYPE100, FNAME 
100 FORIMAT(1X,A5,' ?, TRY AGAIN'/) 
GO TO 101 
bh CALL IFILE(20,FILECIK+5)) 
CALL 1 FILE(21, FNAME) 
- CALL IFILE(22, 'ONE') 
CALL IFILE(23, 'MI PNX') 
IFCIK/2*2.NE.1K) CALL PEI LE Ce epee 
K.EQ. 1FILE(23, 'TFEPNX 
tL ETERS MCE LEC tes109) 
READ (20, 20)(((DTHFAC(I,u,K),J=1,60),1=1,15),K=1,2) 
READ(21,5)(YNPE(J),J#1,60) 
READ(22,5)(((POP(I,d,K),J=1, 60), 1=1,19),K=1,2) 
READ(23,5)(PNIND(J),J=1,60) 
FORMAT (10G) 


THIS SECTION SCALES THE FIRST AGE GROUP FOR THOSE 22-24 YEARS OF 

AGE. THE FLAT RATE CALCULATION BEGINS WITH MULTIPLYING THE 

DISABILITY PREVALENCE RATES BY THE CORRESPONDING POPULATION STATISTICS. 
THE CALCULTION OF THE EARNINGS RELATED PORTION IS SIMILAR BUT IT ALSO 
INCLUDES THE DEATH RATE OF THE CORESPONDING POPULATION AGE GROUPS. 
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SCALE=1. 

IFCI.EQ.1) SCALE=.6 

DO 7 J=5,60 

N=(J-5)/5 

N=HAxX0(1,N) 

N=MINOC6,N) 

YEAR(J)=YEAR(J)+1965 
DISFR(I,u,K)=POP(1+4,U,K)*PRVRTECI,N, 2) *SCALE 


7 DISERN(1,Jd,K)=POP(1+4,JU,K)*PRVRTE(I,N,1)*SCALE*DTHFAC(1,d,K) 
IFC IK/2*2.EQ.1K) GO TO 8 
1il=4 
GO TO 9 
8 IMI=0 
9 JJ=0 
IFCIK.EQ.5) GO TO 10 
GO» 10.45 
C 
C IFCIK.EQ.5), THEN THE TXFPB OPTION HAS BEEN SPECIFIED AND THE 
C USER MMUST DECIDE WHICH PENSION INDEX HE WILL USE.THE INPUT IS 
C THROUGH THE TELETYPE. 
C 
10 TYPE11 
11 FORMAT(' ECONOMIC ASSUMPTIONS - (RS/MI)"',$) 
ACCEPT12,ANS 
12 FORMAT(A2) 


coe) ©). Cec) fC) oD 


IFCANS.NE.'MIE') IMI =4 


THIS SECTION MULTIPLIES THE APPROPRIATE PROPORTION INSURED 

FOR DISABILITY BY THE DISABLED CONTRIBUTOR POPULATION, THE 
RESULTS ARE PUT INTO THE MATRIX WIFEBN WHICH IS THE POPULATION 
ARRAY (NO LONGER USED AND THUS EQUIVALENCED TO WIFEBN) BECAUSE 
THE CALCULATION OF WIVES' BENEFITS IS IDENTICAL TO THAT OF THE 
FLAT RATE TO THIS POINT. 


Ls DO. 15. K=l, 2 


14 
15 


CC) «3 C3 3.6) Co 


[op] 


16 
31 
32 


DO 14 1=1,9 
DO 14 J=5,60 
N=1 


[FCK.EQ.1.AND.U.GE.25) N=2 
IF(K.EQ.2.AND.J.GE. 1) N&3 
IFC(K.EQ.2.AND.JU.GE.15) N=4 
IFCK.EQ.2.AND.U.GE.25) N=5 
IFC(K.EQ.2.AND.U.GE.45) N=6 
DISFR(I,u,K)=DISFR(1,JU,K)* 

J,1 


PROP(1,JUJU+N) 
) 


WIFEBN(I1,JU,1)=DISFR(I,u, 

Ud=Ju+ IMI 

THIS SECTION TOTALS THE EARNINGS RELATED AND FLAT RATE BENEFITS 
THE COMPUTED GO TO SETS THE APPROPRIATE FLAT RATE BENEFIT TO 

BE USED. THE USER CAN SPECIFY A NEW FLAT RATE TO BE USED IN THE 
FUTURE VIA THE TELETYPE.FOR THE YEARS PRECEDING THE NEW FLAT 
RATE THE PROGRAM DEFAULTS TO 300 OR 960 AMOUNTS DEPENDING ON THE 
BEING TESTED.IF THE USER IS TESTING A TXFPB OPTION THE PROGRAM WILL 
ASK WHICH INITIAL FLAT RATE HE WISHES TO USE, 

TYPEG 

FORMAT(' ANNUAL FLAT RATE BENEFIT AND YEAR IT BEGINS? '/) 
ACCEPTS, OAS(3),YR 

DO 16 K=1,2 

DO 16 J=5,60 

DO 16 1=2,9 


DISFR(1,U,K)=DISFR(1,JU,K)+DISFR(1,d,K) 
DISERN(1,u,K)=DISERN(1,JU,K)+DISERNCI, 
GO TO (31,31, 32, 32,33),1K 
| J=1 
GO TO 34 
| J=2 
GO TO 34 


J,K) 


33 TYPE35 
35 FORMAT(' FLAT RATE - (CURRENT/PROPOSED)',$) 
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ACCEPT12,FR 
IF(FR.EQ.'CU') GO TO 31 
GO TO 32 
34 TYPE29 
29 FORMAT(' WIVES BENEFIT - (YES/NO):',$) 
ACCEPT12,ANS 
1 L=2 
IFCANS.EQ."NO') IL=1 


THE EARNINGS RELATED PORTION IS MULTIPLIED BY THE RATE OF GROWTE 
OF RETIREMENT PENSIONS AND THE PROPORTION OF CONTRIBUTORS INSURED 
FOR DISABILITY. 


qooQqonan 


DO 17 K=1,2 
DO 17 J=5,60 
N=1 
JJ=1 
IFC IlJ.EQ.2.AND.JU.GE.8) JJ=2 
Ll=JJ**35 
IF(K.EQ.2.AND.JU.GE. 1) N=2 
IFC(K.EQ.2.AND.JU.GE.15) N=3 
1FCK.EQ.2.AND.JU.GE.25) N=4 
DISERN(1,d,K)=(CCYMPE(JU-2)+YMPE(U-1)+YMPE(JU))/15000)* 
+ (DISINS(N)*DISERN(1,J,K))) 
IF(J.GE. CYR-1966+1)) JuJ=3 
IF (J.GE. (YR-1966+1)) I|1=YR-1966+1 
17 DISFR(1,U,K)=( (OAS (Jd) *(PNIND(J)/PNINDC IT) ))*DISFR(1,u,K)) 
DO 18 J=5,60 
DISERN(1,J,1)=DISERN(1,J,1)+DISERN(1, J, 2) 
18 DISFR(1,JU,1)=DISFR(1,U,1)+DISFR(1,d, 2) 
TYPE180 
180 FORMAT(' RATIO OF IMPUTED RETIREMENT PENSION AND'/ 
+ ' YEAR IT BEGINS: ',$) 
ACCEPT5,Bu, YR1 
IFCIK.EQ.5) GO TO 36 
GO FTO" SS 7 
ao 1 ¥PES50 
350 FORMAT(' EARNINGS RELATED - (CURRENT/PROPOSED)',$) 
ACCEPT12,FR 
| J=2 
IFCFR.EQ.'CU') [Jel 


EARNINGS RELATED DISABILITY PENSIONS ARE SCALED BY 
THE RATIO FOR THE IMPUTED RETIREMENT PENSION. 


LT AP ere TF, 


7,59 3¢1 J=5,60 


IFClJ.EQ.2.AND.U.GE.8) AJ=l1. 
IF(U.GE. (YR1-1965+1)) Aud=Bu 
181 DISERN(1,U,1)=DISERN(1,J,1)*Ad 
WRITE(24, 19) (DISERN(1, J, 1),J=5,60) 
-WRITE(24, 19)(DISFR(1, 3, J=5, 60) 
13 FORMAT (' ,10(1PF12.0)) 
20 FORMAT(6C). 
Cc 
C THIS SECTION CALCULATES THE WIVES’ BENEFITS. 


CALL OFILE(25,FILECIK+15)) 
GO TO (30,22),1L 


C 
C THIS SECTION CALCULATES THE NUMBER OF WIVES ENTITLED TO BENEFITS. 
C 


22 DO 23 1=1,9 
DO 23 J=5,60 
93 WIFEBN(I,J,1)=WIFEBN(I,JU,1)*FACTOR(1) 
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C THIS CALCULATES THE NUMBER OF WIVES OVER 65 WHO ARE 
C ENTITLED TO BENEFITS UNDER THE SCHEME. 
c 

NO 24 1=1,4 

DO 24 J=5,60 

N=(J-5)/5 


N=MAX0(1,N) 
N=MINO(10,N) 
WIFEBNC I,J, 2)=WIFEBN(1+9,J, 2)*FIITH65(1,N) 
24 WIFEBN(10,JU,1)=0. 
C 
C THE WIVES OVER 65 WHO ARE ENTITLED TO BENEFITS ARE ADDED YO THOSE 
C UNDER 65 WHO RECEIVE BENEFITS. 
C 
DO 25 1=1,4 
DO 25 J=5,60 
25 WIFEBN(10,J,1)=WIFEBN(10,U,1)+WIFEBN(I,u, 2) 
DO 250 J=5,60 
250 WIFEBN(10,JU,1)=WIFEBN(10,U,1)*100 
C 
C THIS CALCULATES THE NUMBER OF HUSBANDS AND WIVES AT HOME 
C WITH CHILDREN. 
C 
DO 26 1=1,10 
DO 26 J=5,60 
26 WIFEBN(I,J,1)=WIFEBN(I1,JU,1)*PHWHWC(1I ) 
C 
C THE NUMBER OF WIVES ENTITLED TO BENEFITS IS TOTALLED AND 
C COSTED AND MULTIPLIED BY THE APPROPRIATE PENSION INDEX. 
C 
DO 27 1=2,10 
DO 27 J=5,60 
27 WIFEBN(1,J,1)=WIFEBN(1,JU,1)+WIFEBN(I,u,1) 
DO 270 J=5,60 
JJ=2 
11=8 
[F(UJ.GE.YR-1966+1) JJ=3 
[F(J.GE.YR-1966+1) | 1=YR-1966+1 
270 WIFEBN(1,J,1)=WIFEBN(1,J,1)*(OAS (Jd) *(PNIND(J)/PNINDCI1))) 
WRITE(25,190)(WIFEBN(1,J,1),J=5,60) 
190 FORNAT(' ',10(-2PF12.0)) 
_ GO TO 1000 
30 DO 134 J=1,60 
134 WIFEBN(1,J,1)=0.0 


WRITE(25,190) (WIFEBN(1,U,1),J=1,60) 
1000 TYPE4O,FILE(IK+10),FILECIK+15) 
bO FORMAT(' OUTPUT FILES ARE ',A5,'.DAT AND ',A5,'.DAT'/) 


STOP 
END 
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VERSION AS OF 12:00 AUGUST 11,1972 
PROGRAMME CHILD.8.THIS PROGRAMME CALCULATES BENEFITS 
FOR ORPHANS AND CHILDREN OF DISABLED CONTRIBUTORS. 


DIMENSION DUMMY(10) 
REAL POP(4,60, 2), PROB(4,13),ADU(5), PNIND(60),SDFIFD(4,6), 

- FILE(10), PROB1(4, 58), ORPHBN(4, 58, 2), CHLDBN(4, 58, 2), PROB2(4, 58) 
INTEGER YEAR(60),YR 
DATA ((PROB(I,J),J#1,13),1=1,4) 

+ /.105,.425,.528,.517,.486,.468, .437,.423,.403,.4, 
.399,.395,.401,.15499,.67,1.71,1.773,1. 711, 1.648, 1. 574, 1.502, 
1.424, 1. 38, 1.378, 1. 376,1. 382,.234, .984, 2.589, 3.563, 3.571, 3.427, 
3,296, 3.156, 3.023, 2. 888, 2. 796, 2.805, 2. 801, . 324, 1. 383, 3.577, 
4.947, 5.717, 56522, 56297, 5.097, 4. 878,4.672, 4.466, 4. 576, 4. 567/ 
DATA ((SDFIFD(I,J),J=1,6), 1=1,4) 

+ /,399,.488,.525,.541,.546,.548,.655,.811,.877,.909, 

+ .923,.929,.986, 1.232, 1. 338,1.392,1.42,1.432,1.301,1.636,1.778, 

+ 1.853,1.895,1.916/ 

DATA AGG / 175,68, 5857695 695/ 
DATA FILE/'RSCUR', 'MICUR', 'RSPRO', 'MIPRO', 'TXFPB', 'RSCCH', 

+ "MICCH', 'RSPCH', 'MIPCH', 'TFPCH'/ 

TYPE1 

FORMAT(' TYPE IN PLAN TO BE TESTED: ',$) 
ACCEPT2, FNAME 

FORMATCAS) 

DO 3 IK=1,5 

IECFEWAME.EQ.FILECIK)) GO TO (31,32,31,32,33),1K 
TYPES, FNAME 

EFORMAT(1X,A5,'?, TRY AGAIN'/) 

GO TO 6 


= 
+ 
+ 
+ 


31 CALL IFILE(20, "RSPNX') 


GO TO 34 


32 CALL IFILE(20, 'MIPNX') 


GO TO 34 


33 CALL IFILE(20, 'TFPNX') 
34 CALL IFILE(21, 'ONE') 


co 


Aaa o wo 


CALL OFILE(22,FILE(IK+5)) 
READ(20,9)(PNIND(J),J=1, 60) 

READ (21,9)((POP(1,U,1),J=1,60),1=1,4) 
DO 8 J=1,90 

READ(21,9) DUHNY 

READ(21,9)((POP(I,u, 2),JU=1,60),1=1,4) 


- FORMAT(10G) 


THIS SECTION DOES A LINEAR INTERPOLATION FOR THE 
INTER YEAR PROBABILITIES OF BEING ORPHANED FOR CHILDREN 
FROM 0 TO 17 YEARS OF AGE. 


DO 11 I=1,4 

JN=0 

DO 11 JK=1,12 

N=5 

IF(JK.EQ.1) N=2 
DELTA=PROB( 1, JK+1)-PROB(I, JK) 
XINC=DELTA/N 

DO 11 JL=1,N 

JN=JN+1 

IF(JL.GT.1) GO TO 10 
PROD1(1,JUN)=PROBC 1, JK) 
GO TO ll 


10 PROB1(1,JN)=PROB1(1,JUN-1)+XI NC 


IF(JN.EQ.57) PROB1(1,58)=PROB(1,13) 


11 CONTINUE 

C 

C THE FIRST AGE GROUP IS SCALED BY .6 BECAUSE IT CORTAINS 

© PERSONS AGED 15 TO 19 WHILE WE ARE ONLY CONCERNED WITH THOSE 
C WHO ARE 15 TO 17. 

C 


pO 12 J=1,58 


12 PROB1(4,U)=PROB1(4,U)*.6 
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C 
C THE PROBABILITIES OF BEING ORPHANED AE;NOW MULTIPLIED 
C BY THE CORRESPONDING POPULATION. 
C 
DO =; K=1,2 
DO 13 /=1,4 


DO 13 J=1,58 
13 ORPHBN(I,JU,K)=PROB1(1,U)*POP(1,J+2,K) 


C 
C THE NUMBERS OF ORPIIANS ARE TOTALLED FOR EACH YEAR. 
C 
DO 14 K=1,2 
DO 14 J=1,58 
DO 14 1=2,4 
14 ORPHB(1, JK) =ORPHBN(1,JU,K)+ORPHBN(I,JK) 
4¥ PEELS 
15 FORMAT(' ENTER ANNUAL VALUE OF ORPHANS BENEFIT 
+ "AND YEAR IT BEGINS ',$) 
ACCEPTS, PAYMANT, YR 
1 FC PAYMNT.EQ. 348. )G0 TO 180 
1FCPAYMNT.NE.300.) GO TO 160 
GO TO 222 


C 
C THE NEXT SECTION MULTIPLIES THE ANNUAL ORPHANS BENEFIT BY 
C THE NUMBER OF ORPHANS AND THEN ESCALATES THE RESULTS BY 

C THE APPROPRIATE YEAR'S PENSION INDEX. 

C 

1 


80 DO 181 K=1,2 
DO 182 J=1,6 
182 ORPHBN(1, J, K)=ORPHBN(1,J,K)*300. *(PNIND(J+2)/100.) 
DO 181 J=7,58 
181 ORPHBN(1,U,K)=ORPHBN(1,U,K)*PAYMNT 
GO TO 170 
222 DO 16 K=1,2 
DO 16 J=1,58 
YEAR(J)=1967+d 
16 ORPHBN(1, J, K)=ORPHBN(1,JU,K)*300. *(PNIND(J+2)/100) 
GO TO 179 
160 LOOP=YR-1968 
DO 162 K=1,2 
DO 161 J=1, LOOP 
YEAR(J)=1967+d 
161 ORPHBN(1, JU, K)=ORPHBN(1,J,K)*300. *C(PNIND(J+2)/100) 
DO 162 J=LOOP+1,58 
YEAR(J)=1967+d 
162 ORPHBN(1,U,K)=ORPHBN(1, JU, K)*PAYMNT*( PNIND(J+2)/PNINDCLOOP+3)) 


C 
C THE BENEFITS ARE THEN TOTALLED BY SEX AND STORED FOR PRINTING. 
C 
170 DO 17 J=1,58 
ORPHBN(2,u,1)=ORPHBN(1,J,1)+ORPHBN(1,d, 2) 
IE(J.GE.3.AND. J. LE. 7)ORPHBN(2,J, 1)=ORPHBN(2,J,1)*ADJU(J-2) 
17 COKTINUE 


THE NEXT SECTION DOES A LINEAR INTERPOLATION TO CALCULATE 
THE INTER YEAR PERCENTAGES OF SURVIVING PISABLED FATHERS 
INSURED FOR DISABILITY.THE FIRST DATA YEAR 1S 1975. 


3). ta t3.¢c9 


DO 19 I=1,4 
JN=0 
DO 19 JK=1,5 
DELTA=SDFIFD(1,JUK+1)-SDFIFD(1, JK) 
XINC=DELTA/5 
DO 19 JL=1,5 
JN=JN+1 
1F(JL.GT.1) GO TO 18 
PROB2(1,uUN)=SDFIFD(1,JK) 
GO TO 19 
18 ppo32(1,UN)=PROB2(1,UN-1)+XI NC 
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19 CONTINUE 


C 
© FOR THE YEARS BEYOND 2000 WE USE THE VALUES FOR 2000. 
C 
DO 20 1=1,4 
PROB2(1,26)=SDFIFD(1,6) 
DO 20 J=27,51 
20 PROD2(1,J)=PROB2(1, 26) 


C 
C MULTIPLY THE PERCENTAGES BY THE CORRESPONDING 
C POPULATION TO CALCULATE THE NUMBER OF CHILDREN OF 
C DISABLED FATHERS. 
Cc 
DO 21 K=1,2 
SCALE=1. 
DO 21 1=1,4 


DO 21 J=8,58 
IFC1.EQ.4) SCALE=.6 
21 CHLDBN(1,u,K)=PROB2(1,U-7)*POP(1,J+2,K)*SCALE 
TYPEZ2Z 
22 FORMAT(' ENTER ANNUAL VALUE OF CHILDRENS BENEFIT a 
+ "AND YEAR IT BEGINS ',$) 
ACCEPTS, PAYMNT, YR 


C 
C THE NEXT SECTION FIRST CALCULTES THE 1970 VALUE OF 
C ELIGIBLE CHILDREN BY TAKING .75 OF THE 1975 VALUE. 
C THE VALUES IN BETWEEN THESE YEARS ARE THEN INTERPOLATED 
C USING A TWO POINT LAGRANGEAN PROCEDURE. 
C 

DO 23 K=1,2 

DO 23 1=1,4 
93 CHLDBN(I,3,K)=CHLDBN(I, 8,K)*. 75 

DO 24 K=1,2 

DO 24 11,4 

DO 24 J=4,7 

JK=J-3 

AJ=((5.-JK)/5.) 

BueJK/5. 


“ CHLDBN(I,u,K)=AU*CHLDBN(1, 3,K)+BU*CHLDBN(1, 8, K) 


2 

C 

C THE NEXT SECTION TOTALS THE NUMBER OF CHILDREN ELIGIBLE 
C TO RECEIVE BENEFITS AND THEN CALCULATES THE ANNUAL 

C- COST OF THESE BENEFITS AND ESCALATES TEE RESULTS BY THE 
C APPROPRIATE YEAR'S PENSION INDEX. 

C 


DO 25 K=1,2 
DO 25 J=1,58 
DO 25 1=2,4 
fie CHLDBN(1, JU, K)=CHLDBN(1, JU, K)+CHLDEN(1,u,K) 
}ECPAYMNT.EQ. 348. )GO TO 1180 
4 ] . 
LE CPAYINT NE. 300. ) GO TO 260 


1180 DO 1181 K=1,2 
po 1182 v=1,6 
1182 CHLDBN(1, J, K)=CHLDBN(1,Jd,K) #300. *(PNIND(J+2)/100. ) 
DO 1181 J=7,58 
1181 CHLDBN(1,J,K)=CHLDEN(1,J,K)*PAYMNT 
GO TO 1270 
1222 DO 26 K=1,2 
DO 26 J=1,58 
26 CHLDBN(1,u,K)=CHLDBN(1, Ju, K)#300. *(PNIND(J+2)/100) 
GO TO 1270 
260 LOOP=YR-1968 
DO 262 K=1,2 
DO 261 J=1,LO0P 
061 CHLDBN(1,u,K)=CHLDBN(1,J,K)#300. *(PNIND(J+2)/100) 
DO 262 J=LOOP+1,58 
962 CHLDBN(1, J, K)=CHLDBN(1, J, K)*PAYMIUT#( PNIND (u+2)/PHIND(LOOP+3)) 
1270 DO 27 J=1,58 
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C VERSION AS OF AUGUST 28,1972 


C PROGRAMME FUND.9 
CANADA PENSION PLAN -- FUND CALCULATIOR 
WRITTEN BY HARRY NEWTON, 


C 
C 


* 


* 


1 


1, 


2 
3 
4 


12 
1001 


2001 


10 

1002 

id 

15 

1003 
* 


2002 
HF 


MAY 1972. 


REAL MATURE, NETREV, NETRV1, NRGROW 
INTEGER RECALL 


DIMENSION STABLE(2,4),REGION(12,4),CURPRO(2, 2) 
DIMENSION TOTAL(56),XINFO(7), 


PROGRAM 


INFOC(7), ANNINT(60), AILESPY(56) 


DIMENSION FNAME(5), RTNAM(5), DBNANI(S), WDNAN(S5), DSA (5), WEHAMNCS), 
CHNAM(5), CNNAM(5), FNOUT(5), TOTINV(60), MATURE( 56), CRATE(56) 


DIMENSION RTBEN(56), DHBEN(56), ERWID (5G), FRWID(56), ERDSB(56), 
FRDSB(56), WEDSB(56), ORPHA .N(56), CHILD(56), CPPCN(5 6), RATE(56) 
CNNAN/'RSCCN", 
FNAME/'RSCUR', 
RTNAM/RTRSC', 
DBNAM/'DBRSC', 
WDNAM/'RSCWD', 
DSNAM/'psRsc', 
WFNAM/'RSCWF', 
CHNAM/' RSCCH', 
FNOUT/'RSCFD! 
1 TAB, | PAGE, KURPRO/1, 0, 1/, CNRATE/3. 6/, 1XAUST/9/ 
TOTINV/383. 3,764.5, 718. 7,810. 
ANNINT/20.9, oe: 6, 712, 
CNST69/690600000, /, RATE6S/ 1. 0737/, TOT6/0./, TOT7/2684500000./ 
TRGROW, NRGROW, PROV/2*0., 
e122), | 
J), J=1,4),1 
'NFLAT! 
"h 


DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


TICK, 
2*! 
unett 


CALL CLOCK1 


(CCURPRO CI 
((STABLECI, 
"MODER', 
DATA (CREGIONCI, 
VE EDW', 'ARD 
, 'QUEBE', 
'TCHEW', 


, "SASKA', 


TAY, 


TYPE 1001 
FORMAT(/,' 
ACCEPT 2001, PLAN 
FORMAT(A5) 


DO 10 


1=1,5 


pea 


tAbE les 


'CANAD!, 


MHiCcHt, 
'MICUR", 


"RINIC!, 
'pBMICc', 
"M1CWD', 
"psiic', 
"MICWF', 
"MICCHS, 


, ‘HICED*, 


Jd, Us 1,4), 
»'SLAND', 


on! 
TAN!, 
tAY, 


Ct, 


2x! 


URSPCNT, 
'RSPRO!, 


'eTRSPL, 
'pBRSP!, 
™RSPWD!, 
'DSRSP!, 
"RSPWE', 
ML PCH, 


"RSPCH! 
'RSPFD!, 


"HEPCN', 
"MI PRO!, 
"RTMIP!', 
'DBMIP!, 
"MI PWwD!', 
'DSMIP!, 
"M1 PWE', 


"MIPFD’, 


56*0./ 


17h 3, 56%0. / 


Lag 
7. 2)/" 


=1,2)/' REASO', 


CUR’, 


,12)/"NEWFO", 


'ONTAR 
4 ALBER', 
» "CANAD! 4 


F 

=1 
OVA ' 
v 


IFC PLAN. EQ. FNAME(I)) GO TO 11 


CONTINUE 


TYPE 1002, PLAN 
FORMAT(/, ! ERROR - NO PLAN CALLED ',A5) 


GO TO 12 


NFILE=] 


GO.T0-(13,14,13, 14,15), NFILE 
TYPE 1003 


ACCEPT 2002,ANSWER 
FORMAT(A1) 
IFCANSWER.EQ.'Y') GO TO 14 


KSTAB=1 


CALL IFILE(29,"RSINT') 


GO 


14 KSTAB=2 


16 


100% FORMAT(/,* TYPE IN POPULATION SEGMENT USING 
'] = NFLD; 
'7 - MAN; 


1 
2 


CALL 


TO 16 


[FILE(29, 


TYPE 1004 


GX, 
6X, 


2°=.PEls 
8 - SASK; 


THE INT') 


3 - 


QO «= 
y 


1S TXFPB PLAN MODERATE 
yey ate 1S REASONABLE STABILITY) --- 


NS; 
ALTA; 


, ‘SCOTT: 


hk =- NB; 
10 - BC; 


ie tiot, 


IN PLAN TO BE TESTED - ',$ 


(TAT, 
CALLES*, 


‘TEPCN'/ 
VTXF PB? / 
"RITEP!/ 
"DBTFP'/ 
'TFPWD'/ 
"NSTFP'/ 
"TF PW '/ 


'TFPCH'/ 


'TFPFD'/ 


"RENT', 
"NABLE', 


, UND LA', 
Tar, 


Qe 


) 


get 


'PROP', 


'ND', 


"OSED'/ 


" S1AS > 


5° PR 


SLATS 


Pies 


, NEW AO ae 


: , "MANIT!, TOBA', 
, *BRETE, "SH CO", 


"S» QUES, 


INFLATION? CY OR WD 


",$) 


5 - QUE; 
ll - 


36X,'12 - CANADA LESS QUE (DEFAULT !S 12) <== 


ACCEPT 2003, 
2003 FORIMAT(104) 
1F(1POP.LT.1.OR. 


| POP 


1F (1 POP. EQ.6)PROV=. 55 


DO 1 


l=1,h 


|POP.GT.12) 


TOTINV V(I)= =TOTINV( 1 )*PROV*10**6 


| POP=12 


'seEC'/ 


; 


THESE CODES:' 

6 - ONT;' 
CANADA; 
",9) 


PY Pet ge 
ak. 
Py ee 


7 


ee 
ee 
peg ' 


* 


on 2 J 
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4 ANNINTCL)=ANNINTC1L) ®PROV#*®10**6 


oO QOooaOMoo 0 


CNST69=CNST69*PROV 
TOT7=TOT7*PROV 
TYPE 1017 
1017 FORMAT(/,' IF EXPENSES EXCEED REVENUE, WHICH LOANS WILL BE RECAL', 
1 'LED?',/,' (TYPE -1 FOR MOST RECENT; 0 FOR NONE; 1 FOR Ae 
2 ‘EARLIEST): ‘',$) 
ACCEPT 2003, RECALL 
TYPE 1016 
1016 FORMAT(/,' IN THE FOLLOWING QUESTIONS, TYPE IN ALL RATIOS AND', 
* /," PERCENTAGES AS NUMBERS BETWEEN 0 AND 1.') 
TYPE 1005 
1005 FORMAT(/,' TYPE IN YEAR''S BASIC EXEMPTION - ',$) 
ACCEPT 2003, YBE 
IFCYBE.LT.1.) YBE=YBE*100 
CALL IFILE(29, 'OPSHN') 
READ(29,1006) D1,D2,D3,D4,CRATE 
1006 FORMAT(6(/,10G)) 
TYPE 1007 
1007 FORMAT(/," RETIREMENT BENEFITS',/,1H ,19('="),//,1H , "TYPE IN RET 
LIREMENT BENEFIT RATIO, AND PERCENTAGE OF',/,1H , ‘LOWEST AER''S TO 
2DROP OUT --- ',$) 
ACCEPT 2003, BENRAT, AEROUT 
AEROUT=AEROUT#*100 
TYPE 1015 
1015 FORMAT(/,' IS THERE AN EARNINGS TEST? (Y OR N) - ',$) 
ACCEPT 2002, ERTEST 
TYPE 1008 
1008 FORMAT(/,' DISABILITY BENEFITS',/,1H ,19(¢'-")) 
TYPE 1009 
1009 FORMAT(/," TYPE IN ANNUAL FLAT RATE (IN DOLLARS), AND YEAR IT , 
* 'BEGINS',/,11X,$) 
ACCEPT 2003, XINFO(1),INFO(1) 


TYPE 1010 
1010 FORMAT(/," TYPE IN PERCENTAGE OF BENEFIT THAT IS EARNINGS RELATED, 
«'/,1H , "AND YEAR IT BEGINS --- ',$) 


ACCEPT 2003, XINFO(2), INFO(2) 
XINFO(2)=X1INFO(2)*100 


TYPE 1011 

1011 FORMAT(/,' CHILDREN''S BENEFITS',/,1H ,19('-")) 
TYPE 1009 
ACCEPT 2003,XINFO(3),1NFOC3) 

- TYPE 1012 

1012 FORMAT(/,' WIVES'' BENEFITS',/,1H ,15('=')) 
TYPE 1009 
ACCEPT 2003, XINFO(4), INFOC4) 
TYPE 1013 

1013 FORMAT(/,' WIDOWS'' BENEFITS',/,1H ,16('=')) 
TYPE 1009 
ACCEPT 2003,XINFO(5),INFO(5) 
TYPE 1010 


ACCEPT 2003,XINFO(6),1NFO(6) 

XINFO(6)=XINFO(6)*100 

TYPE 1014 | 
1014 FORMAT(/,' ORPHANS'' BENEFITS’, /,2H »17C'="2) 


TYPE 1009 
ACCEPT 2003, XINFO(7), INFOC(7) 


READ-IN SECTION. 


D1,D2,D3,D4 ARE DUMMY VARIABLES USED TO STORE DATA FOR YEARS 
PRIOR TO 1970. 


READ [N EFFECTIVE ANNUAL RATE OF INTEREST ON NEW INVESTMENTS. 
READ(29,2003) RATE 
READ IN TOTAL CONTRIBUTIONS. 
CALL 1FILE(29, CNNAMCNFILE)) 
READ(29, 2004) D1,02,03,04&,CPPCN 
2004 FORMAT(5G@) 
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C READ IN RETIREMENT BENEFITS. 

CALL IFILE(29,RTNAIMCNFILE)) 
READ(29,2004) RTGSEN 

C READ IN DEATH BENEFITS. 

CALL I1FILE(29,DBNAM(CNFILE)) 
READ(29, 2004) DHBEN 

Cc READ IN WIDOWS' BENEFITS -- EARNINGS RELATED AND FLAT RATE. 
CALL IFILE(29,WDNAMCNFILE)) 

READ(29, 2005) D1,D2,ERVID,D3,D4,FRWID 
2005 FORMAT(8G) 

C READ IN DISABILITY BENEFITS -- EARNINGS RELATED AND FLAT RATC. 
CALL IFILE(29,0SNAMC(NFILE)) 
READ(29,2003) ERDSB 
READ(29, 2003) FRDSB 

C READ IN WIVES' BENEFITS. 

CALL IFILE(29,WFNAMCNFILE)) 
READ(29,2003) WFDSB 

C READ IN CHILDREN'S AND ORPHANS’ BENEFITS. 
CALL IFILE(29,CHNAMC(NFILE)) 
READ(29,2003) D1,D2,O0RPHAN 
READ(29,2003) D1,D2,CHILD 


C OUTPUT SECTION. 
C 


300 CALL OFILE(28,FNOUT(NFILE)) 
CALL TIMECTCLOCK) 
CALL TDATE(LDAY, LIIONTH, LYEAR) 
|FCNFILE.GT.2) KURPRO=2 
WRITE(28,4016) LDAY, LMONTH, LYEAR, TCLOCK 
WRITE(28, 4009) (STABLE(KSTAB, J), J=1,4), (CURPRO(KURPRO, J),J=1, 2), 
* (REGION(1 POP, J),J=1,4) 
IFCNFILE.EQ.5) WRITE(28, 4020) 
IFCYBE.LT.100.) GO TO 301 
WRITE(28,4010) YBE 
GO TO 302 
301 WRITE(28,4011) YBE 
302 WRITE(28,4012) CNRATE, ISTCNR 
WRITE(28, 4023) 
WRITE(28,4013) BENRAT, AEROUT 
ERTST='YES' 
IFCERTEST.EQ.'N') ERTST="NO' 
WRITE(28,4018) ERTST 
WRITE(28,1008) 
WRITE(28,4014) XINFO(1),!NFO(1) 
WRITE(28,4015) XINFO(2),1NFO(2) 
WRITE(28, 1011) 
WRITE(28,4014) XINFO(3), 1 NFOC3) 
WRITE(28,1012) 


WRITES, 40243 XINFO(4), 1 NFOC4) 
WRITE(28, 1013 


WRITE(28,4014) XINFO(5), INFO(5) 

WRITE(28, 4015) XINFO(6),1INFO(6) 

WRITE(28, 1014) 

WRITE(28,4014) XINFO(7), INFO(7) 

WRITE(28, 4019) 

DO 400 1=1,56 

| YEAR=1+1969 

1F(1.EQ.1.0R.1.£Q.29) GO TO 401 

GO TO 402 

401 1 PAGE=I PAGE+1 

WRITE(28, 4017) 

WRITE(28, 4005) 1 TAB, LDAY, LMONTH, LYEAR, TCLOCK, | PAGE 
WRITE(28, 4004) 

RITEC28,4003) (STABLE(KSTAB, J), J=1,4), (CURPROCKURPRO, J), J=1, 2), 
x (REGION(! POP, J), J=1, 4) 

iECNEILE.EQ.5) WRITE(28, 4020) 

WRITE(28, 4000) 

402 EXPENS=,001*CPPCN(I)/CRATE(I) 

TOTALC| JERTBEN(! )+FRDSB(1)+ERDSBE( 1 )+CHILD C1 )+NFDSD C1 )+FRWIDCI D+ 
é ERWID(1 )+ORPHAN(I )+DHBEN( 1 )+EXPENS 
FECMOD(1,4).EQ.1) WRITE(28,4002) 
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KOO WRITE(28,4001) | YEAR, RTBEN(1),FRDSB(1),ERDSB(1), CHILD(1),WFDSACI), 

* FRWID(1), ERWID(1),ORPHAN( I), DRBENCI), EXPENS, TOTALC! ) 

| TAB=2 

| PAGE=0 

DO 500 |1=1,56 

1 YEAR=1+1969 

1FC1.EQ.1.0R.1.EQ.29) GO TO 501 

GO TO 502 

501 | PAGE=! PAGE+1 

WRITE(28,4017) 

WRITE(28,4005) 1 TAB, LDAY, LHONTH, LYEAR, TCLOCK, | PAGE 

WRITE(28, 4008) 

WRITE(28, 4003) (STABLE(KSTAB, J), J=1,4), (CURPRO(KURPRO,J),J=1,2), 
* (REGION(1 POP, J), J=1, 4) 

IFCNFILE.EQ.5) WRITE(28,4020) 

WRITE(28, 4006) 

502 |RECAL=1+4 
IF(RECALL.EQ.1) IRECAL=1 
ALESPY( 1 )=ANNINTC1+3) 
1F(1.GT.17) ALESPY(1)=AIESPY(I )-ANNINT(I-17) 

MATURECI )=0. 
1F(1.GT.16) MATURE(I! )=TOTINV( I-16) 

504 TOTINV(I1+4)=CPPCN(I )+Al ESPYCI )+MATURE( 1 )-TOTALCI D 
IFC(TOTINV(I+4).GT.0..OR.RECALL.EQ.0) GO TO 503 
IRECAL=IRECAL+RECALL 
IFC IRECAL.GT. (1+3).0R. IRECAL.LT.1) GO TO 503 
MATURE( 1) =MATURE(I )+#TOTINVC(IRECAL) 

TOTINVCIRECAL)=0. 
GO TO 504 
503 TOT6=TOT6+MATURE(CI ) 
TOT7=TOT7+TOTINV(I +4) 
XINT=TOTINV(1+4)*RATECI)/100. 
ANNINTC 144) =ANNINTCI4+3)4+X0NT 
COL12=TOT7-TOT6 
ADJ13=ANNINTCI+4) 
1FC(1.GT.16) ADJ13=ADJ13-ANNINTC I-16) 
COL13=COL12+.25*ADJ13 
COL14=CNST69 *RATE69** 1 
COL15=.125*(COL14+TOTAL(I )-CPPCN(I)) 
FUND=COL13+COL15 
IE CFUND. LT.0..AND. 1 XAUST.EQ.0) WRITE(28,4022) IYEAR 
[FCFUND.LT.O.) | XAUST=1 
IF CFUND.LT.O.) FUND=0. 
IF(1.GT.16) TOTINV(I-16 


)=0 
IF(MOD(1,4).EQ.1) WRITE(28,4002) 
500 WRITE(28,4007) IYEAR,CPPCN(1),TOTAL(CI), FUND 
| TAB=3 j 
| PAGE=0 


DO 600 1=1,56 
| YEAR=1+1969 
1FC(1.EQ.1.0R.1.EQ.29) GO TO 601 
GO TO 602 
601 | PACE=I PAGE+1 
WRITE(28, 4017) 
WRITE(28, 4005) | TAB, LDAY, L!ONTH, LYEAR, TCLOCK, I PAGE 
WRITE(28,4024) 
WRITE(28, 4003) (STABLE(KSTAB, J), J=1,4), (CURPRO(KURPRO, J),J=1, 2), 
* (REGION(1 POP, J),J=1,4) 
IECNFILE.EQ.5) WRITE(28,4020) 
WRI TE(28,4025) 
602 TOTREV=CPPCN( 1 )+ALESPY(1 )+MATURECI ) 
NETREV=TOTREV-TOTAL(I ) 
IF(1.EQ.1) GO TO 603 
TRGROW= (TOTREV-TOTRV1) #100. /ABS(TOTRV1) 
NRGROW=(NETREV-NETRV1) #100. /ABS(NETRV1) 
603 TOTRV1=TOTREV 
NETRVI=NETREV 
1F(MOD(1,4&).EQ.1) WRITE(28, 4002) 
600 WRITE(28, 4026) | YEAR, CPPCN(1),AIESPY(1 ),MAATURECI ), TOTREV, TROP OW, 
* TOTAL( I), NETREV, NRGROW 
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WRITE(28,4017) 
T=CLOCK(0,0)*60. 
TYPE 4021,FNOUT(NF!ILE),T 

4000 FORMAT( /,25X,'------- DISABILITY PENSIONS---------- WIDOWS'' PE 
INSIONS', 25X, 'EXPENSES',/,' CALENDAR RETIREMENT',61X, "ORPHANS"! 

2 DEATH', 8X, 'OF',/,1H ,2X, 'YEAR', 6X, 'PENSIONS FLAT EARNINGS °', 
3" CHILDREN''S wWIVES'' FLAT EARNINGS BENEFITS  BENEFITS' 
4, 3X, "ADMINIS=",5X,*TOTAL',/,1H ,24X,"RATE RELATED BENEFITS' 

5, 3X, 'PENSIONS RATE RELATED', 26X, 'TRATION') 

4001 FORMATC1H ,2X,14,6X, -6PF7.1,2(3X,-6PF6.1),2(5X,-6PF6.1), 3X, -EPFE.1 
* ,2X,-6PF8.1,5X,3(-GPF7.1,4X),-6PF8.1) 

4002 FORMAT(2X) 

4003 FORMAT(/,1X,4A5,1X, 2A4,' PLAN',T85, "REGION: ',4A5) 

4004 FORMAT( /,41X,'CANADA PENSION PLAN',/,1H ,30X, 'BENEFITS AND EXPEN 
peeent ADMINISTRATION',/,1H ,32X, "CALL FIGURES IN MILLIONS OF DOLL 
2ARS 

4005 FORMAT(1H1,'"TABLE NO.',12,20X,*RUN DATED',13,'/',12,'/',12, 3X, 

* "TIME: ',A5,20X,'PAGE', 12," OF 2') 
4006 FORMAT(/,' CALENDAR', 7X, ‘CONTRI BUTIONS', 8X, 'BENEFITS &',9X, 
* "FUND AT',/,' YEAR', 31X, 'EXPENSES', 9X, 'YEAR END") 

4007 FORMAT(1H ,2X,14,9X, -GPF11.1, 8X, -6PF11.1,6X,-6PF11. 1) 

4008 FORMAT(/,41X, "CANADA PENSION PLAN',/,1H ,41X,'FUND PROJECTIONS’, 
* /,1H ,32X,'(ALL FIGURES IN MILLIONS OF DOLLARS)') 

4009 FORMAT(/, ' CANADA PENSION PLAN‘,/,1H , "SCHEME: ',4A5,1X, 
IZAR oAinPEAN: pfsah » "REGION: ‘',4A5,//,1H , "ASSUMPTIONS FOR 
2 THIS SCHEME') 

4010 FORMAT(/,' YEAR''S BASIC EXEMPTION IS_',F6.0," DOLLARS') 

4011 FORMAT(/,' YEAR''S BASIC EXEMPTION IS',F6.1," PER CENT') 

4012 FORMAT (1H , "TOTAL CONTRIBUTION RATE IS ',F6.2,' PER CENT, BEGINNING 
*G IN',15) 

4013 FORMAT( /," RETIREMENT BENEFIT RATIO IS ',F6.3,/,1H , "PERCENTAGE 0 
*F LOWEST AER''S TO DROP OUT IS ',F5.1,' PER CENT') 

4014 FORMAT(/,' ANNUAL FLAT RATE IS ',F6.0," DOLLARS, BEGINNING IN', 

* 

4015 FORMAT(1H ,'PERCENTAGE OF BENEFIT THAT IS EARNINGS RELATED Is L 
* F5.1,' PER CENT, BEGINNING IN',15) 

4016 FORMAT (1H1, ‘RUN DATED?» V3 60s a 12e Ty ee She TIME: "AS, SUK, PAGES 
» OF 1") 

4017 FORMAT( /,' TABLES PREPARED BY TAXATION AND FISCAL POLICY BRAK 


SALULPREY. RE vERAAEHEY or “SAREL G9 3° / BH 29% aD INTERGOVERKN ENT 


4018 FORMAT(1H , "EARNINGS TEST? - ',A3) 

kO19 FORMAT(//,2X) 

4020 FORMATC1H+,35X,'- TXFPB PROPOSALS") 

4021 FORMAT(//," END OF RUN',//,' OUTPUT FILE IS ',A5,'.DAT',//, 
* ' EXECUTION TIME = ',F8.3,' SECONDS',///) 


‘4 


AL 


4022 FORMAT(/,1X,16,15X, '**** CPP FUND IS EXHAUSTED xexxe!) 

b023 FORMAT(/,' RETIREMENT BENEFITS',/,1H ,19('-')) 

4024 FORMAT(/,41X, "CANADA PENSION PLAN',/,1H ,40X, ‘REVENUE PROJECTIONS' 
* ,/,1H ,32X,'CALL FIGURES IN MILLIONS OF DOLLARS) ') 

4025 FORMAT(/,' CALENDAR CONTRIBUTIONS | NTEREST LOAN', 6X, 

~ 1 "TOTAL REVENUE GROWTH IN BENEFITS & NET REVENUE',4&X, 

9 ‘GROWTH IN',/, 3X, 'YEAR', 32X, "REPAYMENTS ', 18X, ‘TOTAL REVENUE', 
3 4X," EXPENSES', 18X, "NET REVENUE’, /, 72X, °(%)", 39X, °€39") 

402G FORMATC1H ,16,5X,-6PF11. 1, 2(3X,-6PF9.1),5X,-6PF11.1,6X, OPF7.1, 
* 2X,2(3X,-6PF11.1),6X, OPF7.1) 
STOP 

END 
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C VERSION AS OF AUGUST 17,1972 
C PROGRAMME CPPBN.10 


REAL YHPE(60), RETSN(60),WIDOW(60),WIDFLT(60),WIDERN(E9) 
REAL LMPSH(60),DISABL(60),DISWIF(60),MATUR(60),DISFLT(69) 
REAL DISERN(60), CHLDBN(60), YR(60), CPI (120), 0AS(60),GIS(60) 
DATA MATUR/O0'0> ol, «Zoe Be ey 057 Oe ely Oss 9, 9081. / 


DATA CPI/6*1. 0,1. 036, 1.034, 1.032, 51*1. 03,6*1.,1.028,1.022,1.01t 
&,51*1.01/ 


150 
101 


CALL OFILE(29, 'BENFT') 

TYPE 100 

EORMAT(' SPECIFY PLAN TO BE TESTED: ',$) 

ACCEPT 150,FNANE 

FORMAT(AS) 

TYPE 101 

EORMAT(' SPECIFY PENSIONER''S LIFETIME AVERAGE EARNINGS RATIO: : 


&, 3) 


103 


104 


&$ ) 


108 


ACCEPT 151, AER 

FORMAT(G) 

TYPE...1G2 

FORMAT(' SPECIFY CEILING ON PENSION INDEX IN 1973 ET SEC: pet 
ACCEPT 151,CEI LNG 

TYPE 105 

FORMAT(' SPECIFY RETIREMENT RATE: ',$) 

ACCEPT 151,RETRT 

wire 104 

FORMAT(' SPECIFY WIDOW''S MONTHLY FLAT RATE IN 1973: "D9 
ACCEPT 151,WIDFLT(8) 

TAPE. 105 

FORMAT(' SPECIFY WIDOW''S EARNINGS RELATED RATIO: ar 
ACCEPT 151,WIDRT 

TYPE 106 

EQORMAT(' SPECIFY ORPHAN''S MONTHLY FLAT RATE IN 1973: pee 
ACCEPT 151, CHLDBN(8) 

wire 10s 

FORMAT(' SPECIFY DISABILITY MONTHLY FLAT RATE IN 1973: aP 


ACCEPT 151,DISFLT(8) 
TYPE 108 
FORMAT(' SPECIFY DISABILITY EARNINGS RELATED RATIO IN 1973 ET ', 


&'SEQ: ' 


109 


&,$ 


110 


&2. 


210 


,$) 
ACCEPT 151,DISRT 
TYPE 105 
FORMAT(' SPECIFY MONTHLY BENEFIT FOR WIFE OF DISABLED IN 1973: ' 
) 
ACCEPT 151,DISWIF(8) 
CALL 1FILE(20, FNAME) 
READ(20,110) (YMPE(J),J=1,60) 
FORMAT (10G) 
DO 200 J=8,60 
RETBN(d ) =RETRT*AER*MATUR (J) #( (YMPE (J) +YMPE(U-1)+YMPE(J-2))/3. 0/4 


WIDERN( J) =WIDRT*®(RETBN(U)*(1. /MATURCU)) 
DISERN(J)=DISRT*(RETBN(J)*(1. /MATURCU) ) 
LIUPSM(J)=. 1¥YMPECJ) 

DO 210 J=9,60 
WIDFELT(J)=WIDFLT(J-1)*CEI LNG 
DISFLT(J)=DISFLT(JU-1)*CEI LNG 

CHLDBN(J )=CHLDBN(J-1)*CEI LNG 
DISWIF(J)=DISWIF(J-1)*CEI LNG 


7 iY 


C DO 215 J=9,60 
C215 CHLDBN(J)=29. 


220 


230 
JJ=0 


DO 220 J=8,60 
WI DOW(J)=W1 DERN(J)+HIDFLT (J) 
DISABL(J)=DISERN(J)+DISFLT(J) 
po 230 J=1,60 
YR(J)=1965+d 


IECENAME. EC. 'RSCUR'. OR, FNAME, EQ. 'RSPRO' )JJ=60 
0AS(7)=82. 88 

G1S(7)=67.12 

DO 240 J=8,60 
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OAS (J )=OAS(J-1)*CPI (J-1+Ju) 
240 GIS(J)=GIS(U-1)*CPI (J-1+Jd) 
WRITE(29, 298) FNAME, AER, CEILNG, RETRT 
298 FORMAT(' PLAN: ',A5,5X,! AER: 
&" CEILING: ',F4.2,5X,' RETIREMENT RATE: 
WRITE(29,299) WIDRT,DISRT 
299 FORMAT(' WIDOW''S EARNINGS RELATED RATIO: 
&"ABILITY EARNINGS RELATED RATIO: 
WRITE(29, 300) 
300 FORMAT(' 
&iM 
WRITE(29, 301) 


",F4.2,//) 


MAXI IAUM 


MAXIMUM MAXIMUM MONTHLY 


301 FORMAT(! YEAR''S MAXIMUM MONTHLY 
&'y MAXIMUM MONTHLY PENSION TO 
WRI TE(29, 302) 
302 FORMAT(' YEAR OF PENSIONABLE RETIREMENT 
&ltts LUMP SUM DISABILITY WIFE OF 
&' MONTHLY ') 
WRITE(29, 303) 
303 FORMAT(' PENSION EARNINGS PENSION 
&"N PAYMENT PENSION DISABLED 


a! OF A6s. Gs Sold 
WRITE(29, 304) 
lt alacant 
& le 


DO 350 J=8,60 


WRI TE(29, 310) YR(J), YMPE(U), RETBN(J), WIDOW(CJ), 


&, DISABL(J),DISWIF (J), CHLDBN(J), OAS (J), GIS CJ) 


ee ee i} PENS | ON INDEX', 
Fas ty sf 


'F5.3,5X,' DIS’, 


MAX TIMUM MAXIMU', 
MAXIMUM MONTHLY‘) 
MONTHLY MONTHL', 
PENSION TO') 
WiDOW''S ORPHAN', 
CHILD OF', 
PENSION PENSIO', 
DISABLED', 


CHLDBN(J), LMPSMCd ) 


310 FORMAT(2X,F 5.0, 5X,F 7.0, 7X,F 7.2, 5X,F7.2, 4X, FG. 2, 5X, F 6.0, 5X, F 7.2, 
&8X,F6. 2, 12X,F6. 2, 7X,F6.2, 2X, F6. 2) 
350 CONTINUE 
STOP 


END 
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